Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,218,448 A→C 7.9% intergenic (+35/‑148) aceA → / → aceK isocitrate lyase/isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,218,4480AC 7.9% 519.4 / 26.5 147intergenic (+35/‑148)aceA/aceKisocitrate lyase/isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (67/68);  new base C (5/7);  total (72/75)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

CACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAA  >  NC_000913/4218303‑4218596
                                                                                                                                                 |                                                                                                                                                    
cACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATaa                                                                                                                                                      >  1:368401/1‑146 (MQ=255)
cACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATaa                                                                                                                                                      <  2:368401/146‑1 (MQ=255)
 aCCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCg                                                                                                                                                  <  4:1806971/149‑1 (MQ=255)
 aCCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCg                                                                                                                                                  >  2:198142/1‑149 (MQ=255)
            gTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAg                                                                                                                                                     <  1:652327/135‑1 (MQ=255)
             tGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                      >  3:1392762/1‑149 (MQ=255)
              gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                     <  1:1487252/149‑6 (MQ=255)
                        tACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCgg                                                                                                                           >  2:246374/1‑149 (MQ=255)
                           ttCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGAGCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                 >  1:1640816/1‑140 (MQ=255)
                           ttCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                     <  2:1621096/136‑6 (MQ=255)
                           ttCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                     >  1:1621096/1‑131 (MQ=255)
                           ttCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGc                                                                                                                        >  3:1533621/1‑149 (MQ=255)
                              gATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                     <  3:758931/133‑6 (MQ=255)
                                  aaGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGc                                                                                                                 >  4:1544077/1‑149 (MQ=255)
                                    gTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                   >  3:839789/1‑129 (MQ=255)
                                        cgacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAg                                                                                                                                                     <  3:1429451/107‑1 (MQ=255)
                                        cgacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAg                                                                                                                                                     >  4:1429451/1‑107 (MQ=255)
                                         gacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCt                                                                                                               <  2:588688/144‑1 (MQ=255)
                                         gacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCt                                                                                                               >  1:588688/1‑144 (MQ=255)
                                            tattatTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGAc                                                                                                       >  4:720722/1‑149 (MQ=255)
                                             attattCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACg                                                                                                      >  4:543708/1‑149 (MQ=255)
                                              ttattCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGc                                                                                                     >  4:802980/1‑149 (MQ=255)
                                                 ttCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                 <  2:39343/118‑1 (MQ=38)
                                                 ttCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                 >  1:39343/1‑118 (MQ=38)
                                                 ttCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGa                                                                                                             <  4:967186/138‑1 (MQ=255)
                                                 ttCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGa                                                                                                             >  3:967186/1‑138 (MQ=255)
                                                 ttCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTg                                                                                                  >  4:850948/1‑149 (MQ=255)
                                                 ttCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTg                                                                                                  >  1:679219/1‑149 (MQ=255)
                                                    aGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGc                                                                                                          >  2:642236/1‑138 (MQ=255)
                                                    aGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgcgc                                                                                               >  4:1872684/1‑149 (MQ=255)
                                                      ggcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGt                                                                                             <  2:679219/149‑1 (MQ=255)
                                                       gcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTc                                                                                            <  3:413323/149‑1 (MQ=255)
                                                        cggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCt                                                                                           <  3:850948/149‑1 (MQ=255)
                                                        cggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCt                                                                                           <  3:1812940/149‑1 (MQ=255)
                                                           cACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTAt                                                                                        <  1:246374/149‑1 (MQ=255)
                                                            aCGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                        <  2:1764553/100‑1 (MQ=255)
                                                            aCGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                        >  1:1764553/1‑100 (MQ=255)
                                                            aCGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGcc                                                                                                                >  4:1858340/1‑124 (MQ=255)
                                                              gTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCAt                                                                                     >  2:1529899/1‑149 (MQ=255)
                                                              gTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCAt                                                                                     >  4:58855/1‑149 (MQ=255)
                                                                   cAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGc                                                                                                                                                   <  4:1243369/82‑1 (MQ=255)
                                                                   cAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGc                                                                                                                                                   >  3:1243369/1‑82 (MQ=255)
                                                                     gTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACa                                                                              <  1:1782287/149‑1 (MQ=255)
                                                                      tCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAg                                                                             >  3:1531869/1‑149 (MQ=255)
                                                                             gcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGc                                                                      <  2:192222/149‑1 (MQ=255)
                                                                                gACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGa                                                                   >  2:1127983/1‑149 (MQ=255)
                                                                                   cGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                     <  4:180228/80‑1 (MQ=255)
                                                                                     gCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTa                                                              <  3:67425/149‑1 (MQ=255)
                                                                                      cTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAg                                                             <  4:1531869/149‑1 (MQ=255)
                                                                                       tCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCt                                                                                                               <  3:228928/98‑1 (MQ=255)
                                                                                       tCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGc                                                                      <  2:294870/139‑1 (MQ=255)
                                                                                          aCTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                         <  3:1870406/149‑1 (MQ=255)
                                                                                            tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTAt                                                                                        >  1:421705/1‑116 (MQ=255)
                                                                                            tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTa                                                              >  1:558816/1‑142 (MQ=255)
                                                                                            tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTa                                                              <  3:1449386/142‑1 (MQ=255)
                                                                                            tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTa                                                              <  2:558816/142‑1 (MQ=255)
                                                                                            tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa                                                       <  1:1127983/149‑1 (MQ=255)
                                                                                            tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa                                                       >  4:9229/1‑149 (MQ=255)
                                                                                             gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt                                                      <  1:1529899/149‑1 (MQ=255)
                                                                                             gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt                                                      <  3:58855/149‑1 (MQ=255)
                                                                                               agaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATaa                                                    >  4:1614305/1‑149 (MQ=255)
                                                                                               agaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATaa                                                    <  3:1614305/149‑1 (MQ=255)
                                                                                               agaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATaa                                                    >  4:745016/1‑149 (MQ=255)
                                                                                                gaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAg                                                   <  3:1872684/149‑1 (MQ=255)
                                                                                                gaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAg                                                   >  3:1855044/1‑149 (MQ=255)
                                                                                                 aaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAgg                                                  <  1:198142/149‑1 (MQ=255)
                                                                                                   tCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGt                                                                                             <  3:1585208/104‑1 (MQ=255)
                                                                                                   tCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGt                                                                                             >  4:1585208/1‑104 (MQ=255)
                                                                                                     gCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGtt                                              <  1:782475/149‑1 (MQ=255)
                                                                                                       aGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTa                                            <  4:1855044/149‑1 (MQ=255)
                                                                                                              aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCt                                                                                                               <  2:549442/75‑1 (MQ=255)
                                                                                                              aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCAt                                                                                     <  3:1782902/101‑1 (MQ=255)
                                                                                                              aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCAt                                                                                     >  2:1488965/1‑101 (MQ=255)
                                                                                                              aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCAt                                                                                     <  1:1488965/101‑1 (MQ=255)
                                                                                                                 aaccacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGt                                                                                                                          <  4:1395179/61‑1 (MQ=255)
                                                                                                                 aacaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGt                                                                                                                          >  3:1395179/1‑61 (MQ=255)
                                                                                                                 aacaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATcc                                  >  3:1157054/1‑149 (MQ=255)
                                                                                                                     acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCa                              >  1:1514759/1‑149 (MQ=255)
                                                                                                                     acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCa                              >  3:232532/1‑149 (MQ=255)
                                                                                                                      caacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAAc                                                                 >  2:121618/1‑113 (MQ=255)
                                                                                                                       aacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGc                                                 <  1:876409/128‑1 (MQ=255)
                                                                                                                       aacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGc                                                 >  2:876409/1‑128 (MQ=255)
                                                                                                                        acCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGc                                                                      >  4:1277410/1‑106 (MQ=255)
                                                                                                                         ccGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtc                          >  2:1721215/1‑149 (MQ=255)
                                                                                                                         ccGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtc                          >  3:1751754/1‑149 (MQ=255)
                                                                                                                         ccGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtc                          >  4:324187/1‑149 (MQ=255)
                                                                                                                          cGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtct                         <  3:324187/149‑1 (MQ=255)
                                                                                                                           gTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtctc                        <  2:720311/149‑1 (MQ=255)
                                                                                                                           gTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtctc                        <  4:1508146/149‑1 (MQ=255)
                                                                                                                           gTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtctc                        <  4:1751754/149‑1 (MQ=255)
                                                                                                                           gTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATtctc                        <  4:232532/149‑1 (MQ=255)
                                                                                                                             tGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTg                      <  3:994353/149‑1 (MQ=255)
                                                                                                                                tGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAc                                           >  3:34318/1‑125 (MQ=255)
                                                                                                                                tGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTc                   >  3:473830/1‑149 (MQ=255)
                                                                                                                                    tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtga               <  1:1202097/149‑1 (MQ=255)
                                                                                                                                    tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtga               <  4:207401/149‑1 (MQ=255)
                                                                                                                                       aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtgatga            <  3:745016/149‑1 (MQ=255)
                                                                                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGaa                                                                  >  1:337822/1‑94 (MQ=255)
                                                                                                                                             gATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGgcgc      >  1:945935/1‑149 (MQ=255)
                                                                                                                                               tAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTcc                  <  2:928301/135‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAA  >  NC_000913/4218303‑4218596

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: