Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,218,447 A→G 6.0% intergenic (+34/‑149) 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,4470AG 6.0% 723.7 / 24.2 232intergenic (+34/‑149)aceA/aceKisocitrate lyase/isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (112/106);  new base G (3/11);  total (115/117)
Fisher's exact test for biased strand distribution p-value = 5.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 235 total aligned reads displayed.

TATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAA  >  NC_000913/4218298‑4218595
                                                                                                                                                     |                                                                                                                                                    
tATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATa                                                                                                                                                      >  1:256262/1‑150 (MQ=255)
   ctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAgg                                                                                                                                                   >  1:1649094/1‑150 (MQ=255)
   ctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAgg                                                                                                                                                   >  2:133624/1‑150 (MQ=255)
    tcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGc                                                                                                                                                  >  1:862966/1‑150 (MQ=255)
    tcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGc                                                                                                                                                  >  1:905149/1‑150 (MQ=255)
     cACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCg                                                                                                                                                 >  2:1929226/1‑150 (MQ=255)
     cACCAGCAGGAAGTGGGTACAGGTTACTTCGAAAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGATACACTGAAGAATCGCAGTTCTAAGCAACAAAAACCAATCCTGACTGTAGGCCGGATAAGGCg                                                                                                                                                 >  2:851387/1‑150 (MQ=255)
      aCCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGt                                                                                                                                                <  1:1976417/150‑1 (MQ=255)
          gcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                            <  1:1270141/150‑1 (MQ=255)
             ggAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                         >  1:887732/1‑150 (MQ=255)
             ggAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                         >  2:355186/1‑150 (MQ=255)
              gAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                        <  1:658907/150‑1 (MQ=255)
                  tGGGTACTGGTTACTTGGATAAAGTGAGTACAATTATTCAGGGGGGCACGTCTTCCGCCACCGCGTGCACGGGCTCCACTTATCAATGGCAGTTCTAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                   <  2:1513959/150‑1 (MQ=255)
                   gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                   <  1:1799134/150‑1 (MQ=255)
                   gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                   <  2:1227820/150‑1 (MQ=255)
                   gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                   >  1:296167/1‑150 (MQ=255)
                   gggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                   >  2:525546/1‑150 (MQ=255)
                     gTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                 >  2:531825/1‑150 (MQ=255)
                      tACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                >  1:779948/1‑150 (MQ=255)
                       aCAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                               <  2:257174/150‑2 (MQ=255)
                        cAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                              <  1:1101756/150‑3 (MQ=255)
                          ggTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATc                                                                                                                            <  1:1609717/150‑1 (MQ=255)
                           gTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATaa                                                                                                                                                     >  2:1300949/1‑124 (MQ=255)
                           gTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTGCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctaca                                                                                                                           <  2:1371045/150‑6 (MQ=255)
                               cTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGc                                                                                                                       <  2:1097123/150‑1 (MQ=255)
                                 tCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATc                                                                                                                            <  1:32245/143‑1 (MQ=255)
                                   gATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGa                                                                                                                   >  1:512606/1‑150 (MQ=255)
                                        aGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCa                                                                                                                      >  1:1178272/1‑142 (MQ=255)
                                             cgacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGc                                                                                                         <  1:1333188/150‑1 (MQ=255)
                                              gacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCg                                                                                                        <  1:1036881/150‑1 (MQ=255)
                                                cTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGAc                                                                                                      <  2:1027715/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)
                                                     attCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTg                                                                                                 >  1:1341888/1‑150 (MQ=255)
                                                       tCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgcg                                                                                               <  1:192589/150‑1 (MQ=255)
                                                        cAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgcgc                                                                                              >  2:1780840/1‑150 (MQ=255)
                                                         aGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                                                                             <  1:980739/150‑1 (MQ=255)
                                                             cggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCtt                                                                                         <  2:1373754/150‑1 (MQ=255)
                                                              ggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTa                                                                                        <  1:1400663/150‑1 (MQ=255)
                                                                  cGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCAt                                                                                    >  2:171529/1‑150 (MQ=255)
                                                                     cTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGcc                                                                                 >  2:1960965/1‑150 (MQ=255)
                                                                     cTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGcc                                                                                 <  1:171535/150‑1 (MQ=255)
                                                                      ttCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCt                                                                                <  2:240381/150‑1 (MQ=255)
                                                                       tCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                      <  2:563074/95‑1 (MQ=255)
                                                                       tCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                               <  1:1315919/150‑1 (MQ=255)
                                                                       tCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                               <  2:13161/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)
                                                                          gTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAg                                                                            <  1:1149828/150‑1 (MQ=255)
                                                                            cACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGcc                                                                          <  1:1540199/150‑1 (MQ=255)
                                                                             aCCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCg                                                                                                        <  1:361633/119‑1 (MQ=255)
                                                                             aCCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCg                                                                         >  2:1069219/1‑150 (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                                                                                >  2:913610/1‑139 (MQ=255)
                                                                                  gcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGcc                                                                    <  1:1593390/150‑1 (MQ=255)
                                                                                    tGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGa                                                                  >  1:1937595/1‑150 (MQ=255)
                                                                                     gACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGaa                                                                 >  2:1202315/1‑150 (MQ=255)
                                                                                      aCCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAAc                                                                <  2:1700182/150‑1 (MQ=255)
                                                                                        cGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGt                                                              >  1:1804333/1‑150 (MQ=255)
                                                                                             ccACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCt                                                         >  1:1897987/1‑150 (MQ=255)
                                                                                              cagtgatgaATGGCTGTTCTATGCAACAACAACCGTTGCTGACTGTAGGGCGGATAAGGCGTTCTCGCCGCATCCGGCATTCGGTGCAGGATGCCTGATGCGACTCTTGGGCGTCTTATCATGCCTTCAGCCGTTGCCGATCGTAGGGTg                                                        <  2:1636940/147‑1 (MQ=25)
                                                                                              cACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                        <  2:296167/150‑1 (MQ=255)
                                                                                              cACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                        >  1:1275971/1‑150 (MQ=255)
                                                                                              cACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                        >  1:67428/1‑150 (MQ=255)
                                                                                                 tGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt                                                     >  1:69286/1‑150 (MQ=255)
                                                                                                  gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATa                                                    >  2:271440/1‑150 (MQ=255)
                                                                                                    agaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAg                                                  <  1:22294/150‑1 (MQ=255)
                                                                                                      aaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGc                                                <  2:407054/150‑1 (MQ=255)
                                                                                                       aTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCg                                               <  2:1898996/150‑1 (MQ=255)
                                                                                                          gCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGttt                                            >  2:176291/1‑150 (MQ=255)
                                                                                                               tCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgcc                                       >  2:108524/1‑150 (MQ=255)
                                                                                                               tCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgcc                                       >  1:353854/1‑150 (MQ=255)
                                                                                                                cTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                        <  2:1878379/132‑1 (MQ=255)
                                                                                                                cTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                        >  1:1878379/1‑132 (MQ=255)
                                                                                                                cTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTg                                                        <  2:1803928/132‑1 (MQ=255)
                                                                                                                cTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgccg                                      <  2:954528/150‑1 (MQ=255)
                                                                                                                 tAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                                                                             <  1:1445751/94‑1 (MQ=255)
                                                                                                                 tAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                                                                             >  2:1445751/1‑94 (MQ=255)
                                                                                                                 tAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgccgc                                     >  2:1944491/1‑150 (MQ=255)
                                                                                                                  aaGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCa                                    >  1:458342/1‑150 (MQ=255)
                                                                                                                     caacaacaaCCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATGCGGCGTAAACGCCTTATcc                                 <  1:724611/150‑1 (MQ=25)
                                                                                                                        caacaccCGTTGCTGACTGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGc                              <  1:1959699/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)
                                                                                                                                    cTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTc                  >  1:1020774/1‑150 (MQ=255)
                                                                                                                                      gACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCt                >  1:175872/1‑150 (MQ=255)
                                                                                                                                        cTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtga              >  1:1826227/1‑150 (MQ=255)
                                                                                                                                          gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtgatg            >  1:132385/1‑150 (MQ=255)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAgg         >  1:545324/1‑150 (MQ=255)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAgg         >  1:1466878/1‑150 (MQ=255)
                                                                                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGgc       <  2:1671953/150‑1 (MQ=255)
                                                                                                                                                    tAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTaa  >  2:1500222/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAA  >  NC_000913/4218298‑4218595

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