Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,218,461 G→T 5.5% intergenic (+48/‑135) 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,4610GT 5.5% 749.8 / 20.7 218intergenic (+48/‑135)aceA/aceKisocitrate lyase/isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (107/99);  new base T (2/10);  total (109/109)
Fisher's exact test for biased strand distribution p-value = 3.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

Only 100 of 232 total aligned reads displayed.

GAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTG  >  NC_000913/4218312‑4218610
                                                                                                                                                     |                                                                                                                                                     
gAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                       <  1:658907/150‑1 (MQ=255)
gAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                       <  1:1328563/150‑1 (MQ=255)
gAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                       <  2:1081895/150‑1 (MQ=255)
     gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                  <  1:1799134/150‑1 (MQ=255)
     gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                  >  1:296167/1‑150 (MQ=255)
     gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                  >  2:525546/1‑150 (MQ=255)
     gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                  >  2:1492853/1‑150 (MQ=255)
     gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                  >  2:1315919/1‑150 (MQ=255)
       gTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                <  1:1376048/150‑1 (MQ=255)
       gTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                                >  2:531825/1‑150 (MQ=255)
        tACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                               >  2:1502834/1‑150 (MQ=255)
        tACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                               <  2:184081/150‑1 (MQ=255)
          cAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                             <  1:1101756/150‑3 (MQ=255)
            ggTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATc                                                                                                                                           >  1:1970059/1‑150 (MQ=255)
                 cTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGc                                                                                                                                      >  1:1886479/1‑150 (MQ=255)
                     gATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGa                                                                                                                                  >  1:512606/1‑150 (MQ=255)
                       tAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGaag                                                                                                                                <  2:862966/150‑3 (MQ=255)
                         aaGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGcc                                                                                                                              <  2:23469/150‑1 (MQ=255)
                             gacgacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGAt                                                                                                                          <  1:133624/150‑1 (MQ=255)
                                gacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCg                                                                                                                       <  1:1036881/150‑1 (MQ=255)
                                      tattCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                             <  2:1995916/122‑3 (MQ=38)
                                      tattCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                             >  1:1995916/1‑120 (MQ=38)
                                       attCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGTCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTg                                                                                                                >  2:1257005/1‑150 (MQ=255)
                                       attCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTg                                                                                                                >  1:1700182/1‑150 (MQ=255)
                                          cAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgcgc                                                                                                             >  1:1898996/1‑150 (MQ=255)
                                          cAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgcgc                                                                                                             >  2:1152180/1‑150 (MQ=255)
                                              gcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCt                                                                                                         <  2:1928970/150‑1 (MQ=255)
                                                     gTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATg                                                                                                  >  1:1191031/1‑150 (MQ=255)
                                                       cTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCgg                                                                                                                                         >  2:1958590/1‑109 (MQ=255)
                                                       cTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGcc                                                                                                <  1:171535/150‑1 (MQ=255)
                                                         tCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                                     >  1:563074/1‑95 (MQ=255)
                                                         tCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                                              <  1:1315919/150‑1 (MQ=255)
                                                          cAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTAc                                                                                             <  1:531825/150‑1 (MQ=255)
                                                           aGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACa                                                                                            <  1:1493778/150‑1 (MQ=255)
                                                            gTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTg                                                                                                                            >  2:1599125/1‑117 (MQ=255)
                                                            gTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTg                                                                                                                            <  1:1599125/117‑1 (MQ=255)
                                                            gTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAg                                                                                           <  1:1149828/150‑1 (MQ=255)
                                                               aCCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCg                                                                                                                       <  1:361633/119‑1 (MQ=255)
                                                               aCCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCg                                                                                        <  1:1007927/150‑1 (MQ=255)
                                                               aCCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTGGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCg                                                                                        <  1:1859265/150‑1 (MQ=255)
                                                                 cgcgCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGtt                                                                                      >  1:233124/1‑150 (MQ=255)
                                                                 cgcgCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTATCATGCCTACAGCCGTTg                                                                                     <  2:733947/150‑1 (MQ=255)
                                                                  gcgcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTg                                                                                     <  1:740937/150‑1 (MQ=255)
                                                                   cgcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCt                                                                                               <  1:913610/139‑1 (MQ=255)
                                                                    gcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGcc                                                                                   <  1:324983/150‑1 (MQ=255)
                                                                    gcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGcc                                                                                   <  1:1593390/150‑1 (MQ=255)
                                                                    gcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGcc                                                                                   <  2:1940860/150‑1 (MQ=255)
                                                                     cTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGccc                                                                                  >  2:1684567/1‑149 (MQ=255)
                                                                       gACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGaa                                                                                >  2:1202315/1‑150 (MQ=255)
                                                                        aCCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAAc                                                                               <  2:1700182/150‑1 (MQ=255)
                                                                         ccGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACg                                                                              >  1:16490/1‑150 (MQ=255)
                                                                          cGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGt                                                                             >  1:1804333/1‑150 (MQ=255)
                                                                           ggCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTa                                                                            <  1:1316410/150‑1 (MQ=255)
                                                                               ccACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCt                                                                        >  2:1797971/1‑150 (MQ=255)
                                                                                cagtgatgaATGGCTGTTCTATGCAACAACAACCGTTGCTGACTGTAGGGCGGATAAGGCGTTCTCGCCGCATCCGGCATTCGGTGCAGGATGCCTGATGCGACTCTTGGGCGTCTTATCATGCCTTCAGCCGTTGCCGATCGTAGGGTg                                                                       <  2:1636940/147‑1 (MQ=25)
                                                                                cACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                                       >  1:1275971/1‑150 (MQ=255)
                                                                                   tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt                                                                    >  1:69286/1‑150 (MQ=255)
                                                                                    gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATa                                                                   <  1:1929226/150‑1 (MQ=255)
                                                                                    gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATa                                                                   >  2:271440/1‑150 (MQ=255)
                                                                                      agaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAg                                                                 >  1:712282/1‑150 (MQ=255)
                                                                                        aaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGc                                                               <  2:407054/150‑1 (MQ=255)
                                                                                         aTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCg                                                              >  2:1343240/1‑150 (MQ=255)
                                                                                         aTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCg                                                              >  1:877812/1‑150 (MQ=255)
                                                                                            gCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGttt                                                           >  2:176291/1‑150 (MQ=255)
                                                                                             cAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTa                                                          >  2:1037231/1‑150 (MQ=255)
                                                                                               gTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcg                                                        >  2:1777908/1‑150 (MQ=255)
                                                                                                 tCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgcc                                                      >  1:353854/1‑150 (MQ=255)
                                                                                                 tCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgcc                                                      <  2:937332/150‑1 (MQ=255)
                                                                                                  cTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgccg                                                     <  2:954528/150‑1 (MQ=255)
                                                                                                   tAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                                                                                            >  2:1445751/1‑94 (MQ=255)
                                                                                                   tAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                                                                                            <  1:1445751/94‑1 (MQ=255)
                                                                                                    aaGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCa                                                   >  1:458342/1‑150 (MQ=255)
                                                                                                     aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCAt                                                  >  1:837436/1‑150 (MQ=255)
                                                                                                        aacaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCg                                               <  2:1088779/150‑1 (MQ=255)
                                                                                                           aacaacCGTTTCTGACTGTAGGCCTGATAAGGAGTTCACGACGCATACGTCAATCGGTGAACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCa                                            >  2:961301/1‑150 (MQ=25)
                                                                                                            acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCaa                                           >  1:1671953/1‑150 (MQ=255)
                                                                                                             caacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAAt                                          >  2:1571055/1‑150 (MQ=255)
                                                                                                                  gTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTctct                                     <  2:256262/150‑1 (MQ=255)
                                                                                                                     gCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCt                                  >  2:696957/1‑150 (MQ=255)
                                                                                                                      cTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTc                                 >  1:1020774/1‑150 (MQ=255)
                                                                                                                        gACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCt                               >  1:175872/1‑150 (MQ=255)
                                                                                                                          cTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtga                             >  1:1826227/1‑150 (MQ=255)
                                                                                                                            gTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtgatg                           <  1:1944491/150‑1 (MQ=255)
                                                                                                                            gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCCCCtgatg                           <  2:1275971/150‑1 (MQ=255)
                                                                                                                             tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtgatga                          >  1:205785/1‑150 (MQ=255)
                                                                                                                              aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTccctcatca                         <  2:293451/150‑8 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAgg                        >  1:1466878/1‑150 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAgg                        >  1:545324/1‑150 (MQ=255)
                                                                                                                                 ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGgc                      <  2:1671953/150‑1 (MQ=255)
                                                                                                                                   ggATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGgcgc                    >  2:1372642/1‑150 (MQ=255)
                                                                                                                                     aTAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTa                  <  1:99546/150‑1 (MQ=255)
                                                                                                                                         ggCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAgg                        <  2:1873490/140‑1 (MQ=255)
                                                                                                                                         ggCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAgg                        >  1:1873490/1‑140 (MQ=255)
                                                                                                                                         ggCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGc             >  1:1170048/1‑150 (MQ=255)
                                                                                                                                              tCACGCCGCATCAGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCcgcg         >  2:1221172/1‑150 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATGCGGCGTAAACGCCTTATcca                                               >  2:1317289/1‑108 (MQ=21)
                                                                                                                                                 cgccgcATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTgg      >  2:298515/1‑150 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTgg      <  1:271440/150‑1 (MQ=255)
                                                                                                                                                   ccgcATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATTCCGCGTGGcc    <  2:1719901/150‑1 (MQ=255)
                                                                                                                                                   ccgcATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGcc    >  2:1879959/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTG  >  NC_000913/4218312‑4218610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: