Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,288 T→A 6.7% K123I (AAA→ATA) ‡ 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,2880TA 6.7% 302.7 / 9.1 105K123I (AAA→ATA) ‡sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (45/53);  new base A (4/3);  total (49/56)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 106 total aligned reads displayed.

AATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGC  >  NC_000913/1724140‑1724435
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                           >  2:243193/1‑150 (MQ=255)
                         ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                           >  2:383792/1‑150 (MQ=255)
                         ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                           >  1:1008694/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)
                                    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: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)
                                                                                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                                                               >  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:350160/150‑1 (MQ=255)
                                                                                                        cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                            <  1:1378970/150‑1 (MQ=255)
                                                                                                        cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                            <  2:1420399/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                                   <  1:1453556/150‑1 (MQ=255)
                                                                                                                 aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                                   >  1:1641073/1‑150 (MQ=255)
                                                                                                                 aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                                   <  2:17920/150‑1 (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)
                                                                                                                            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: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)
                                                                                                                                      cGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtgg              <  2:1578202/150‑1 (MQ=255)
                                                                                                                                       ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCCTGTTTACCGGTATTTTGTGGTTCAAGCTGCCCGCCTGCGGAGTCCGCGGCGCTGGCTTTGCCATCTTtggt             >  1:1273128/1‑150 (MQ=255)
                                                                                                                                       ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggt             >  1:1565913/1‑150 (MQ=255)
                                                                                                                                       ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggt             <  1:565203/150‑1 (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)
                                                                                                                                                aGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctg    >  2:571846/1‑150 (MQ=255)
                                                                                                                                                  cTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTAGTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctggc  <  2:916458/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGC  >  NC_000913/1724140‑1724435

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: