Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,952,587 T→A 5.8% W108R (TGG→AGG)  yecN → MAPEG family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,952,5870TA 5.8% 538.6 / 14.3 152W108R (TGG→AGG) yecNMAPEG family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (67/76);  new base A (4/5);  total (71/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

TGCTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTA  >  NC_000913/1952441‑1952735
                                                                                                                                                   |                                                                                                                                                    
tgCTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTgg                                                                                                                                                    <  1:601928/149‑1 (MQ=255)
  cTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGgtg                                                                                                                                                  >  3:114503/1‑149 (MQ=255)
   tGATGCTGATTATGGAAATGAATGATGTAGAAGCGTGGATGGTGCATATTTGCGGCATCGTTTTGCTAACTGGTTGACTGATGCATTATTCCGGTTTGCCTCACCGTCTGTTCCGCTTGCGGCGTTACGGCAGGAGCGC‑CTGCTGgtgt                                                                                                                                                 <  4:1382593/149‑1 (MQ=255)
    gATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGtgtg                                                                                                                                                >  1:118029/1‑149 (MQ=255)
     aTGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTgc                                                                                                                                               <  1:1437680/149‑1 (MQ=255)
       gCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              <  4:1226513/149‑5 (MQ=255)
        cTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCt                                                                                                                                            >  4:506800/1‑149 (MQ=255)
           tttATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGtt                                                                                                                                         <  1:885540/149‑1 (MQ=255)
             tATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGGTG‑GCGCtcatgc                                                                                                                                       <  4:120805/149‑6 (MQ=255)
             tATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGGTG‑GCGCtcatgc                                                                                                                                       <  4:233678/149‑6 (MQ=255)
                 gAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              <  3:1443175/139‑5 (MQ=255)
                 gAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              >  4:1443175/1‑135 (MQ=255)
                  aaaTGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATg                                                                                                                                  <  2:213514/149‑1 (MQ=255)
                  aaaTGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTGTTGCTTGCTGGTCGTCTGAGGCATTATTACGGTTTTGATCACCGTCTGTTCCGCTGGCGACGGTCCGGCATGAGCGC‑CACCTGGCGTGCGCTGTTGCTGATg                                                                                                                                  >  4:473026/1‑149 (MQ=255)
                   aatgaatgGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATgg                                                                                                                                 <  2:887829/149‑1 (MQ=255)
                     agaatgGCGCAGAAACCCGGATGGTACATATTGGCGGCATCGTGTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTg                                                                                                                               >  4:32755/2‑149 (MQ=255)
                          ggCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGc                                                                                                                          <  1:342647/149‑1 (MQ=255)
                             gcAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGaa                                                                                                                       >  4:479099/1‑149 (MQ=255)
                              cAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAAt                                                                                                                      <  2:304019/149‑1 (MQ=255)
                               aGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATc                                                                                                                     >  2:239264/1‑149 (MQ=255)
                                gAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCt                                                                                                                    >  3:496971/1‑149 (MQ=255)
                                 aaaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCtt                                                                                                                   <  1:66753/149‑1 (MQ=255)
                                  aaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCttt                                                                                                                  <  1:5902/149‑1 (MQ=255)
                                  aaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCttt                                                                                                                  <  2:1296146/149‑1 (MQ=255)
                                  aaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCttt                                                                                                                  <  3:25195/149‑1 (MQ=255)
                                   aCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTg                                                                                                                 >  3:1055447/1‑149 (MQ=255)
                                   aCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTg                                                                                                                 <  1:990165/149‑1 (MQ=255)
                                    ccTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTgg                                                                                                                <  3:1152921/149‑1 (MQ=255)
                                           ggTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGc                                                                                                         <  1:770895/149‑1 (MQ=255)
                                             tGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGccc                                                                                                       <  4:1205590/149‑1 (MQ=255)
                                              gCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGa                                                                                                                                    >  2:1415408/1‑119 (MQ=255)
                                              gCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCt                                                                                                      >  3:1084142/1‑149 (MQ=255)
                                               cATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTg                                                                                                     <  3:180819/149‑1 (MQ=255)
                                                atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTgg                                                                                                    <  4:68257/149‑1 (MQ=255)
                                                atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTgg                                                                                                    >  1:279183/1‑149 (MQ=255)
                                                 tatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTggg                                                                                                   >  2:1014971/1‑149 (MQ=255)
                                                  atTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGa                                                                                                  <  1:974722/149‑1 (MQ=255)
                                                  atTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGa                                                                                                  <  1:231440/149‑1 (MQ=255)
                                                   tttGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAg                                                                                                 >  1:895194/1‑149 (MQ=255)
                                                   tttGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAg                                                                                                 >  4:53608/1‑149 (MQ=255)
                                                     tGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGtt                                                                                               >  2:124127/1‑149 (MQ=255)
                                                      gCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTg                                                                                              <  1:1014971/149‑1 (MQ=255)
                                                      gCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTCCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTg                                                                                              <  2:619559/149‑1 (MQ=255)
                                                        ggCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              >  2:531125/1‑96 (MQ=255)
                                                        ggCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            <  1:10324/149‑1 (MQ=255)
                                                         gCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGtt                                                                                           <  1:1444747/149‑1 (MQ=255)
                                                         gCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGtt                                                                                           >  1:374992/1‑149 (MQ=255)
                                                            tCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTgc                                                                                                                                               >  2:614935/1‑94 (MQ=255)
                                                             cGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTtct                                                                                       <  3:32755/149‑1 (MQ=255)
                                                              gTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTtctc                                                                                      >  4:616641/1‑149 (MQ=255)
                                                               ttttGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTcc                                                                                     <  3:224184/149‑1 (MQ=255)
                                                                   gcttgctGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGc                                                                                 <  3:566795/149‑1 (MQ=255)
                                                                     ttgctGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGt                                                                               >  2:406727/1‑149 (MQ=255)
                                                                         tGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAgc                                                                           >  1:1258094/1‑149 (MQ=255)
                                                                         tGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAgc                                                                           >  1:334325/1‑149 (MQ=255)
                                                                           gtcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            <  2:468352/130‑1 (MQ=255)
                                                                           gtcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            >  1:468352/1‑130 (MQ=255)
                                                                           gtcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAgcgc                                                                         >  1:395756/1‑149 (MQ=255)
                                                                            tcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGca                                                                        <  2:334325/149‑1 (MQ=255)
                                                                                    tGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCTTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGc                                                                <  3:66526/149‑1 (MQ=255)
                                                                                    tGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCAGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGGCCTGGGAGCTGGTTTTCTCCCTGCGTTAGCGCACAATACGc                                                                >  4:1191181/1‑149 (MQ=255)
                                                                                      cATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCa                                                              <  1:29337/149‑1 (MQ=255)
                                                                                       attattACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTg                                                                                                                 <  3:337115/97‑1 (MQ=255)
                                                                                        ttattaCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCAct                                                            >  1:1372842/1‑149 (MQ=255)
                                                                                            taCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttct                                                        >  3:125781/1‑149 (MQ=255)
                                                                                            taCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttct                                                        >  4:377089/1‑149 (MQ=255)
                                                                                              cGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttcttt                                                      <  1:1062043/149‑1 (MQ=255)
                                                                                               ggTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            >  2:494911/1‑110 (MQ=255)
                                                                                               ggTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            <  1:494911/110‑1 (MQ=255)
                                                                                               ggTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCtttt                                                     <  3:1263452/149‑1 (MQ=255)
                                                                                                gTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCttttt                                                    <  4:835751/149‑1 (MQ=255)
                                                                                                 ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTc                                                   >  3:475868/1‑149 (MQ=255)
                                                                                                 ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTc                                                   >  3:1385066/1‑149 (MQ=255)
                                                                                                 ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTc                                                   >  2:336511/1‑149 (MQ=255)
                                                                                                      atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGa                                              <  2:395756/149‑1 (MQ=255)
                                                                                                      atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGa                                              <  2:485185/149‑1 (MQ=255)
                                                                                                       tcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGAt                                             >  4:1338250/1‑149 (MQ=255)
                                                                                                        caCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGAGGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGAtt                                            >  4:1127728/1‑149 (MQ=255)
                                                                                                             tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTAc                                       <  1:239264/149‑1 (MQ=255)
                                                                                                             tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTAc                                       >  3:576979/1‑149 (MQ=255)
                                                                                                                gTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGtt                                    <  3:1338250/149‑1 (MQ=255)
                                                                                                                 ttCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTGTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGACTTTTACGTTa                                   >  3:576621/1‑149 (MQ=255)
                                                                                                                 ttCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTa                                   <  1:435889/149‑1 (MQ=255)
                                                                                                                 ttCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTa                                   >  4:459823/1‑149 (MQ=255)
                                                                                                                  tCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTAt                                  <  4:1084142/149‑1 (MQ=255)
                                                                                                                   ccGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATg                                 <  4:571916/149‑1 (MQ=255)
                                                                                                                           gACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGGGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATg                                 <  2:813979/141‑1 (MQ=255)
                                                                                                                              gTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGa                      <  2:279183/149‑1 (MQ=255)
                                                                                                                               ttCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                     >  1:65091/1‑149 (MQ=255)
                                                                                                                                 ccGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGATTCTTTGGTATATTCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGTACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGacac                   >  4:1069364/1‑149 (MQ=255)
                                                                                                                                   ggCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                 >  2:365662/1‑149 (MQ=255)
                                                                                                                                    gCATGAGCGC‑CACCTGGTGTGCACTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTTTGTCTCACCGCGACACGCt                >  2:1127546/1‑149 (MQ=255)
                                                                                                                                       tGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTAtt             >  4:1028214/1‑149 (MQ=255)
                                                                                                                                        cagcgcacACCAGGTG‑GCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGTCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTAttt            >  3:777538/8‑149 (MQ=255)
                                                                                                                                           cgc‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCt         <  2:271019/149‑1 (MQ=255)
                                                                                                                                             c‑cACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGc       >  3:945018/1‑149 (MQ=255)
                                                                                                                                             c‑cACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGc       >  1:768771/1‑149 (MQ=255)
                                                                                                                                               cACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcc      <  4:907166/149‑1 (MQ=255)
                                                                                                                                                  cTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTAc                                       >  2:760024/1‑113 (MQ=255)
                                                                                                                                                   tGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTa  >  4:1107632/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGC‑CACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTA  >  NC_000913/1952441‑1952735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑