Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,282 G→C 6.8% T125S (ACC→AGC)  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,2820GC 6.8% 334.8 / 13.0 103T125S (ACC→AGC) sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (46/49);  new base C (4/3);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 109 total aligned reads displayed.

ATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGC  >  NC_000913/1724133‑1724431
                                                                                                                                                     |                                                                                                                                                     
attacttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCgg                                                                                                                                                       <  1:1467634/150‑1 (MQ=255)
attacttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCgg                                                                                                                                                       <  1:1481034/150‑1 (MQ=255)
  tacttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTATGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTTATTTAATACAGTGATTTCAGACGAGTCGCGATGACGGCATCGGtg                                                                                                                                                     >  1:1781861/1‑149 (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)
                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                                                                                                                       >  2:243193/1‑150 (MQ=255)
                                ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                       >  1:701721/1‑150 (MQ=255)
                                ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                       >  1:1008694/1‑150 (MQ=255)
                                ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                       >  2:383792/1‑150 (MQ=255)
                                 cgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCCGCGCTTTGTCTTTGATTTCATCGAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcag                                                                                                                      >  1:1175859/1‑150 (MQ=255)
                                    caccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggca                                                                                                                   <  1:758478/150‑1 (MQ=255)
                                       cgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGaa                                                                                                                >  2:1576369/1‑149 (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)
                                                       ttGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGc                                                                                                <  1:1336748/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)
                                                                      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)
                                                                                       ccATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTAtttt                                                                >  1:471299/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                                                           <  2:1322296/150‑1 (MQ=255)
                                                                                            aGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGa                                                           >  1:35580/1‑150 (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                                        <  1:350160/150‑1 (MQ=255)
                                                                                                                aTCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGa                                       >  2:1725195/1‑150 (MQ=255)
                                                                                                                   cAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTc                                    <  1:196236/150‑1 (MQ=255)
                                                                                                                     gTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCcg                                  >  1:725443/1‑150 (MQ=255)
                                                                                                                        aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                               <  2:17920/150‑1 (MQ=255)
                                                                                                                        aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                               >  1:1641073/1‑150 (MQ=255)
                                                                                                                        aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                               <  1:1453556/150‑1 (MQ=255)
                                                                                                                         tttCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcggc                              >  2:1393779/1‑150 (MQ=255)
                                                                                                                             aGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTg                          >  1:381400/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)
                                                                                                                                   tcaataaTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCACTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTg                    >  2:1120049/7‑150 (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)
                                                                                                                                      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)
                                                                                                                                                    ggTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctgg   <  2:337930/150‑1 (MQ=255)
                                                                                                                                                    ggTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTGCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGCTTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctgg   >  2:75106/1‑150 (MQ=255)
                                                                                                                                                     gTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggc  <  1:1092843/150‑1 (MQ=255)
                                                                                                                                                     gTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATACCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggc  <  1:1770921/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGC  >  NC_000913/1724133‑1724431

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: