Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,279 T→A 6.6% D126V (GAT→GTT) ‡ sodC ← superoxide dismutase, Cu, Zn, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,724,2790TA 6.6% 281.8 / 9.1 107D126V (GAT→GTT) ‡sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (47/51);  new base A (4/3);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

Only 100 of 107 total aligned reads displayed.

TAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTG  >  NC_000913/1724131‑1724425
                                                                                                                                                    |                                                                                                                                                  
taattacttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATc                                                                                                                                                   <  2:1509555/150‑1 (MQ=255)
 aattacttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCACCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCg                                                                                                                                                  >  2:1770921/1‑150 (MQ=255)
  attacttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCgg                                                                                                                                                 <  1:1481034/150‑1 (MQ=255)
  attacttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCgg                                                                                                                                                 <  1:1467634/150‑1 (MQ=255)
     acttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAg                                                                                                                                              <  2:242600/150‑1 (MQ=255)
     acttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAg                                                                                                                                              >  1:1322296/1‑150 (MQ=255)
         aattacACCACAGGCGTAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCttt                                                                                                                                          >  2:1092843/1‑150 (MQ=255)
          attacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTg                                                                                                                                         >  2:1643191/1‑150 (MQ=255)
           ttacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGc                                                                                                                                        <  2:1703741/150‑1 (MQ=255)
             acacCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCg                                                                                                                                      <  2:1240532/150‑1 (MQ=255)
                ccacAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCa                                                                                                                                   >  2:1243295/1‑150 (MQ=255)
                ccacAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCa                                                                                                                                   >  1:1213204/1‑150 (MQ=255)
                ccacAGGCATAGCGTGCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTGCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCa                                                                                                                                   <  1:915901/150‑1 (MQ=255)
                  acaGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGCCCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCatt                                                                                                                                 >  2:1453556/1‑150 (MQ=255)
                  acaGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCatt                                                                                                                                 >  1:1824935/1‑150 (MQ=255)
                  acaGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAAAGCTACCGATGCCGTCatc                                                                                                                                 <  1:1277392/150‑2 (MQ=255)
                    aGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAAAGCTACCGATGCCGTCatcgc                                                                                                                               <  2:482666/150‑4 (MQ=255)
                    aGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAAAGCTACCGATGCCGTCatcgc                                                                                                                               <  1:689868/150‑4 (MQ=255)
                     ggCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATTTTATCGCCGCCAACGTGGACGATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGTTAGCTTTGCCGTCattatt                                                                                                                              <  1:1198929/150‑1 (MQ=255)
                     ggCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGATTTGCCGTCattatt                                                                                                                              >  1:1058590/1‑150 (MQ=255)
                       cATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCTTCGGTAGCTTTGCCGTCATTATTga                                                                                                                            >  1:916458/1‑150 (MQ=255)
                        aTAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCGACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                                                                                                           >  1:622830/1‑150 (MQ=255)
                        aTAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCGGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                                                                                                           >  1:1578202/1‑150 (MQ=255)
                         tAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgacg                                                                                                                          >  1:1525155/1‑150 (MQ=255)
                               tCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGt                                                                                                                    <  1:337384/150‑1 (MQ=255)
                                caccgccaccgcCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTg                                                                                                                   >  2:1310991/1‑150 (MQ=255)
                                  ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                 >  1:1008694/1‑150 (MQ=255)
                                  ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                 >  2:243193/1‑150 (MQ=255)
                                  ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                 >  1:701721/1‑150 (MQ=255)
                                   cgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCCGCGCTTTGTCTTTGATTTCATCGAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcag                                                                                                                >  1:1175859/1‑150 (MQ=255)
                                      caccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggca                                                                                                             <  1:758478/150‑1 (MQ=255)
                                          gccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGCCCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATc                                                                                                         >  2:461294/1‑150 (MQ=255)
                                             cAGCGGTTTAGGTTGATCGGACATTTTATCGGCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGcc                                                                                                      <  1:100131/150‑1 (MQ=255)
                                               gCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                                                                                                    <  2:883309/150‑1 (MQ=255)
                                                       ggTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCg                                                                                            <  1:1774601/150‑1 (MQ=255)
                                                        gTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCgg                                                                                           <  1:461632/150‑1 (MQ=255)
                                                         ttGATCGGACATATTATCGCCGCCAACGGGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAGCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGc                                                                                          <  2:1363210/150‑1 (MQ=255)
                                                             tCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCt                                                                                      <  1:1044096/150‑1 (MQ=255)
                                                              cGGACATATTATCGCCGCCACCGTGGCCCATCAGCGCTTTGTCTTTGATTTCATCAAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGCTCGCCTAAATGCCCGGCACCtt                                                                                     >  1:337930/1‑150 (MQ=255)
                                                                gACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCt                                                                                   >  2:833899/1‑150 (MQ=255)
                                                                  cATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTgg                                                                                 <  1:1439763/150‑1 (MQ=255)
                                                                  cATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTgg                                                                                 <  1:583951/150‑1 (MQ=255)
                                                                     attatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGccc                                                                              <  2:1797102/150‑1 (MQ=255)
                                                                      ttatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCt                                                                             <  2:1805073/150‑1 (MQ=255)
                                                                       tatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCtt                                                                            >  1:1281052/1‑150 (MQ=255)
                                                                        atCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTc                                                                           <  1:350468/150‑1 (MQ=255)
                                                                         acgccgcCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGAATTCAGACGAGGCTCGATGACGGAATCCGTATCTTTGCCGTCATTATAGACGACCAGTGCAGGCAGCTCGCCTAAATGCCCGGCTACTTCTGGCCCATCa                                                                          >  2:285200/2‑150 (MQ=255)
                                                                           gccgccAACGTGGACCATCAGCGCTTTGTCTTTGAATTCAACCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTAATGACGACCAGAGCAGTCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTGCATg                                                                        >  2:48162/1‑150 (MQ=255)
                                                                                       gACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTAtt                                                            >  2:23927/1‑150 (MQ=255)
                                                                                          cATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTtg                                                         >  1:1310588/1‑150 (MQ=255)
                                                                                           aTCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGACCGGCAGCTTCTGGCCCTTCATGTTTACCGGTATTTtgt                                                        >  2:905403/1‑150 (MQ=255)
                                                                                            tCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTtgtg                                                       <  1:1741856/150‑1 (MQ=255)
                                                                                             cAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTgg                                                      >  2:571686/1‑150 (MQ=255)
                                                                                              aGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGa                                                     >  1:35580/1‑150 (MQ=255)
                                                                                              aGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGa                                                     <  2:1322296/150‑1 (MQ=255)
                                                                                              aGCGCTTTGGCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGCAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGa                                                     <  1:508817/150‑1 (MQ=255)
                                                                                                cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATc                                                   <  1:383792/150‑1 (MQ=255)
                                                                                                  cTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCaa                                                 <  1:1002277/150‑1 (MQ=255)
                                                                                                  cTTTGTCCTTGATTTCCACCAGTGATTTCAGCCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGAACGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCaa                                                 <  2:1203755/150‑1 (MQ=255)
                                                                                                   tttGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGa                                               <  2:1824935/150‑1 (MQ=255)
                                                                                                    ttGTCTTTGATTTCATGCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGTTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCATAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGa                                               <  1:1310991/150‑1 (MQ=255)
                                                                                                        cTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGCCCAGTGCAGGCAGATCGCCTAAATGCCCGGCCCCTTCTGGCCCTTCATGTTTACAGGTATTTTGTGGATCAAGATGcc                                           >  2:767979/1‑150 (MQ=255)
                                                                                                          ttGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCg                                         <  2:803132/150‑1 (MQ=255)
                                                                                                            ggtttAATACAGTGTTTTGAGAAGAGGCGCGATGACGGCATCGGTTGCGTTGCCGGGATTATTTACGACCCGTGCAGGCAGATAGCCAAAAATCCCCGCCCCGTGGGGGCCATCATGTTTAACGGGTTTTTTTTGTTAAAGAAGcccccc                                       >  2:672291/3‑150 (MQ=255)
                                                                                                             aTTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCt                                      >  1:60912/1‑150 (MQ=255)
                                                                                                                tCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCg                                   <  2:604474/150‑1 (MQ=255)
                                                                                                                 cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                  <  1:1378970/150‑1 (MQ=255)
                                                                                                                 cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                  <  2:1420399/150‑1 (MQ=255)
                                                                                                                 cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                  <  1:350160/150‑1 (MQ=255)
                                                                                                                  aTCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGa                                 >  2:1725195/1‑150 (MQ=255)
                                                                                                                     cAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTc                              <  1:196236/150‑1 (MQ=255)
                                                                                                                      aGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTcc                             <  1:1562167/150‑1 (MQ=255)
                                                                                                                       gTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCcg                            >  1:725443/1‑150 (MQ=255)
                                                                                                                          aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                         <  2:17920/150‑1 (MQ=255)
                                                                                                                          aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                         <  1:1453556/150‑1 (MQ=255)
                                                                                                                          aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                         >  1:1641073/1‑150 (MQ=255)
                                                                                                                           tttCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcggc                        >  2:1393779/1‑150 (MQ=255)
                                                                                                                               aGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTg                    >  1:381400/1‑150 (MQ=255)
                                                                                                                                gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTgg                   >  2:525496/1‑150 (MQ=255)
                                                                                                                                 ccttggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGc                  <  1:194067/146‑1 (MQ=255)
                                                                                                                                  cGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGCCCAGTGCAGGCAGCTCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCt                 >  1:51882/1‑150 (MQ=255)
                                                                                                                                  cGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCt                 >  2:1564256/1‑150 (MQ=255)
                                                                                                                                  cGACGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGGAGATCGCCTAAATGCCCGGCAGCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCt                 <  2:1552859/150‑1 (MQ=255)
                                                                                                                                     ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTg              <  2:381400/150‑1 (MQ=255)
                                                                                                                                     ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTg              >  1:907128/1‑150 (MQ=255)
                                                                                                                                      caataaTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACATTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc             >  2:73782/6‑150 (MQ=255)
                                                                                                                                       cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGcc            <  2:1525155/150‑1 (MQ=255)
                                                                                                                                       cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGcc            <  2:1213204/150‑1 (MQ=255)
                                                                                                                                        gcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCa           <  1:243193/150‑1 (MQ=255)
                                                                                                                                        ataaTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCa           >  1:1539214/4‑150 (MQ=255)
                                                                                                                                           aTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCTAGATGCCCGCCTGCGGATTCCGCGGCGCTTGCttttcaatgt        >  1:159702/1‑143 (MQ=255)
                                                                                                                                            tGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCtt       <  2:581025/150‑1 (MQ=255)
                                                                                                                                             gACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCttt      <  1:551884/150‑1 (MQ=255)
                                                                                                                                              aCGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtg     <  1:41542/150‑1 (MQ=255)
                                                                                                                                              aCGGCAGCGGTCGCGTTGCCGTCATTAGTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGTATTCCGCGGCGCTGGCTTTGCCATCTTtg     <  1:345831/150‑1 (MQ=255)
                                                                                                                                               cGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtgg    <  2:1578202/150‑1 (MQ=255)
                                                                                                                                                ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCCTGTTTACCGGTATTTTGTGGTTCAAGCTGCCCGCCTGCGGAGTCCGCGGCGCTGGCTTTGCCATCTTtggt   >  1:1273128/1‑150 (MQ=255)
                                                                                                                                                ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggt   <  1:565203/150‑1 (MQ=255)
                                                                                                                                                ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggt   >  1:1565913/1‑150 (MQ=255)
                                                                                                                                                 gCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtg  <  1:1243295/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTG  >  NC_000913/1724131‑1724425

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