Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,285 G→C 5.5% A124G (GCT→GGT)  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,2850GC 5.5% 385.9 / 8.4 118A124G (GCT→GGT) sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (61/49);  new base C (4/3);  total (66/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

ACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTG  >  NC_000913/1724136‑1724433
                                                                                                                                                     |                                                                                                                                                    
acttaattacCCGACTGTCATTGCGTTCTCCGCAACCGGCCAGCGGTTTAGGTTGATCGGACAATTTCTCGCCGCCAACGTGGCCCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAg                                                                                                                                                      <  1:157672/150‑1 (MQ=255)
  ttaattacACCACAGGCATAGCGTTCCCCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTGTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCt                                                                                                                                                    >  1:217124/1‑150 (MQ=255)
   taattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCtt                                                                                                                                                   >  1:1919796/1‑150 (MQ=255)
    aattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGCTTTCAGACGAGGCGCGATGACGGCATCGGTAGCttt                                                                                                                                                  >  2:248932/1‑150 (MQ=255)
      ttacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCCTCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGc                                                                                                                                                >  2:821811/1‑150 (MQ=255)
      ttacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGc                                                                                                                                                <  2:1817837/150‑1 (MQ=255)
       tacTCCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGcc                                                                                                                                               <  2:1301508/150‑1 (MQ=255)
       tacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGcc                                                                                                                                               >  2:1925526/1‑150 (MQ=255)
        acacCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCAATGTCTTTGATTTCATCCAGTGATTTCAGAAGAGGCGCGATGACGGCATCGGTATCTTTGccc                                                                                                                                              >  2:721919/1‑149 (MQ=255)
           ccacAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCCTCCAGTGATTTCAGCCGAGGCGCGATGACGGCATCGGGAGCTTTGCCGTCa                                                                                                                                           >  1:1789441/1‑150 (MQ=255)
           ccacAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCa                                                                                                                                           <  1:702817/150‑1 (MQ=255)
                ggCATAGCGTTCACCGCCACCGCCCAGCGGTGGAGGTGGCTCGGACATTTTATCGCCGCCAACGTGGAGCATCAGCGCATTGGCGTTGATTACATCCAGGGATTTCATCCGAGGCGCGATGAAGGCATCGGGAGCATTGCCCTCagtatt                                                                                                                                      >  1:304853/1‑150 (MQ=255)
                 gCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTg                                                                                                                                     <  2:405984/150‑1 (MQ=255)
                   aTAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAAAGCTACCGATGCCGTCatcgcgcct                                                                                                                                   <  2:1175625/150‑8 (MQ=255)
                      gCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGACGTCATTATTgacgac                                                                                                                                >  1:1433740/1‑150 (MQ=255)
                       cGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGAcc                                                                                                                               >  1:1280198/1‑150 (MQ=255)
                        gTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCa                                                                                                                              >  1:873695/1‑150 (MQ=255)
                         ttCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAg                                                                                                                             >  1:1988329/1‑150 (MQ=255)
                         ttCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAg                                                                                                                             <  1:242732/150‑1 (MQ=255)
                            accgccaccgccCAGCGGTTTAGGTTGCTCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgc                                                                                                                          <  2:1417045/150‑1 (MQ=255)
                             ccgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGCCGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                         >  2:1307738/1‑150 (MQ=255)
                             ccgccaccgccCAGCGGTTTAGGTTGATAGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATGTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                         >  1:1433303/1‑150 (MQ=255)
                             ccgccaccgccCAGCGGGTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                                                                                                         <  1:1785307/150‑1 (MQ=255)
                              cgccaccgccCAGCGGTTTAGGTTGCTCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCAGTCATTATTGCCGACCAGTgcag                                                                                                                        >  2:1785252/1‑150 (MQ=255)
                                ccaccgccCAGCGGTTTAGGTTTTTCGGACATATTATCGCCGCCAACGTTGACCATCAGCGCTTTGTATTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggc                                                                                                                      <  2:2017855/150‑1 (MQ=255)
                                ccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggc                                                                                                                      <  1:1428456/150‑1 (MQ=255)
                                     gccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATc                                                                                                                 >  2:102526/1‑150 (MQ=255)
                                        cAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGcc                                                                                                              <  2:240353/150‑1 (MQ=255)
                                               ttAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGc                                                                                                       <  1:905860/150‑1 (MQ=255)
                                                tAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGcc                                                                                                      >  1:1488374/1‑150 (MQ=255)
                                                  ggTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCg                                                                                                    >  2:445501/1‑150 (MQ=255)
                                                   gTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCgg                                                                                                   >  2:1393571/1‑150 (MQ=255)
                                                     tGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCa                                                                                                 >  1:2002102/1‑150 (MQ=255)
                                                      gATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCAc                                                                                                >  2:1435821/1‑150 (MQ=255)
                                                         cGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCtt                                                                                             >  1:1061349/1‑150 (MQ=255)
                                                         cGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCtt                                                                                             >  1:452994/1‑150 (MQ=255)
                                                          ggACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCTCCTTc                                                                                            >  1:342174/1‑150 (MQ=255)
                                                           gACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCt                                                                                           >  2:723250/1‑150 (MQ=255)
                                                           gACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCt                                                                                           <  1:39950/150‑1 (MQ=255)
                                                           gACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCt                                                                                           <  1:719896/150‑1 (MQ=255)
                                                               tattatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGcc                                                                                       >  1:266657/1‑150 (MQ=255)
                                                                 ttatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCCGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATTCCCGGCACCTTCTGGccca                                                                                     >  2:808363/1‑149 (MQ=255)
                                                                    tCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCa                                                                                  <  2:1302938/150‑1 (MQ=255)
                                                                     cgccgcCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCAt                                                                                 >  1:1595795/1‑150 (MQ=255)
                                                                      gccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATg                                                                                >  1:1234870/1‑150 (MQ=255)
                                                                      gccgccAACGTGGACCAGCAGCGCTTTGGCTATGATTTCATCCAGTGATTTCAGACGCGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATg                                                                                <  1:1229731/150‑1 (MQ=255)
                                                                       ccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACTAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGt                                                                               >  1:1797433/1‑150 (MQ=255)
                                                                       ccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGt                                                                               >  1:1918150/1‑150 (MQ=255)
                                                                       ccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGt                                                                               >  1:1480931/1‑150 (MQ=255)
                                                                       ccgccAACATGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGt                                                                               <  2:1161990/150‑1 (MQ=255)
                                                                         gccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGCCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATAGCATAAATGCCCGGCACCTTCGGGCCCTTGATGttt                                                                             >  2:1881508/1‑150 (MQ=255)
                                                                         gccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCACCTAAATGCCCGGCACCTTCTGGCCCTTCATGttt                                                                             <  2:839885/150‑1 (MQ=255)
                                                                           cAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCATCAGGTTTAc                                                                           >  1:381024/1‑150 (MQ=255)
                                                                            aaCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTAcc                                                                          <  1:229628/150‑1 (MQ=255)
                                                                            aaCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTAcc                                                                          >  1:1409265/1‑150 (MQ=255)
                                                                                 ggACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTAt                                                                     >  1:1038757/1‑150 (MQ=255)
                                                                                 ggACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTACATGTTTACCGGTAt                                                                     >  2:855585/1‑150 (MQ=255)
                                                                                  gACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTAtt                                                                    <  1:132691/150‑1 (MQ=255)
                                                                                      aTCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCCCCTTCTGGCCCTTCATGTTTACCGGTATTTtgt                                                                >  1:1967492/1‑150 (MQ=255)
                                                                                       tCAGCGCTTTGTCATTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGACCGGCACCATAAGGCCCTTCATGATTACCGTTATTTtgag                                                               >  1:1196853/1‑148 (MQ=255)
                                                                                         aGCGCTTTGTCTTTGATGTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCACTTCATGTTTACCGGTATTTTGTGGa                                                             <  1:81670/150‑1 (MQ=255)
                                                                                          gcgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGAt                                                            >  1:1993498/1‑150 (MQ=255)
                                                                                           cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATc                                                           <  2:738993/150‑1 (MQ=255)
                                                                                           cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACAGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATc                                                           <  2:1641762/150‑1 (MQ=255)
                                                                                                tGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGAt                                                      >  1:605886/1‑150 (MQ=255)
                                                                                                 gTCTTTGATTTCATCCAGTGAATTCAGACGAGGCGCGAAGACGGCATCGGTACCTTAGCCGTCATTAATGACGACCAGTGCAGGCAGCTCGCCTAAATGCCCGGCAACTTCTGGACATTCAAGTTTACAGGTATTATTTCGATCAAGaaa                                                     >  2:1261356/1‑148 (MQ=255)
                                                                                                   cTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCGGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGcc                                                   >  2:1338431/1‑150 (MQ=255)
                                                                                                     ttGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCg                                                 >  1:1138295/1‑150 (MQ=255)
                                                                                                     ttGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCg                                                 <  2:148383/150‑1 (MQ=255)
                                                                                                     ttGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCAATATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTAATGTTTACCGGTATTTTGTGGATCAAGATGCCCg                                                 >  1:1366564/1‑150 (MQ=255)
                                                                                                       taTTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGcc                                               <  1:448110/149‑1 (MQ=255)
                                                                                                           tCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCg                                           <  1:248932/150‑1 (MQ=255)
                                                                                                            cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                          >  1:1612613/1‑150 (MQ=255)
                                                                                                            cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                          >  1:971596/1‑150 (MQ=255)
                                                                                                              tCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGAt                                        >  1:624463/1‑150 (MQ=255)
                                                                                                               ccAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCAGGCACCTTCTGGCCCTTCAGGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGAGGAtt                                       >  2:1906064/1‑150 (MQ=255)
                                                                                                                 aGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTcc                                     >  2:1827100/1‑150 (MQ=255)
                                                                                                                    gATTTCCGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCcgcg                                  <  1:584683/150‑1 (MQ=255)
                                                                                                                      tttCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcggc                                >  1:345171/1‑150 (MQ=255)
                                                                                                                        tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGgcgc                              >  2:1467476/1‑150 (MQ=255)
                                                                                                                        tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGAATCcgcaggcg                              <  1:763864/150‑6 (MQ=255)
                                                                                                                          tgACGAGGCTCGATGACGGGATCGTTAGCTTTGCCGGCATTATTGACGACCAGAGGAGGCAGATCGGCTAAATGCCCTGCACCTTCTGGCCCTTCATGTTTAGCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGTTTCCGCGGCGCTg                            <  1:57337/149‑1 (MQ=255)
                                                                                                                               aGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCttt                       <  2:1793359/150‑1 (MQ=255)
                                                                                                                                 gcgcGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                     <  2:1061348/150‑1 (MQ=255)
                                                                                                                                  cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGcc                    <  1:1996870/150‑1 (MQ=255)
                                                                                                                                  cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGcc                    <  2:2002101/150‑1 (MQ=255)
                                                                                                                                  cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGcc                    <  1:821812/150‑1 (MQ=255)
                                                                                                                                   ataaTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGCTCAAGATGCCCGCCTGCGGCTTCCGCGGCGCTGGCTTTGCCa                   >  1:702839/4‑150 (MQ=255)
                                                                                                                                     gATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATc                 <  1:1006541/150‑1 (MQ=255)
                                                                                                                                     aaTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCCCCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATc                 >  2:24735/2‑150 (MQ=255)
                                                                                                                                      aTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCt                >  1:436910/1‑150 (MQ=255)
                                                                                                                                         aCGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGAGCGCGTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtg             <  2:1636410/150‑1 (MQ=255)
                                                                                                                                          cGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtgg            >  1:1146013/1‑150 (MQ=255)
                                                                                                                                          cGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTGtgg            >  1:1292217/1‑150 (MQ=255)
                                                                                                                                          cGGCATCGGTAGCTTAGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtgg            <  2:673560/150‑1 (MQ=255)
                                                                                                                                           ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggt           <  1:1273582/150‑1 (MQ=255)
                                                                                                                                            gCATCGGTAGCTTTGCCGTCATTATTGACGACCGGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtg          <  1:1338432/150‑1 (MQ=255)
                                                                                                                                              aTCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggc        <  2:1766751/150‑1 (MQ=255)
                                                                                                                                                    cgCTTTGCCGTCATTATTGTGGACCTGAGCAGGCAGATCGGCTTAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTCTTTTGTGGATCTAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctg  <  1:1319746/149‑1 (MQ=255)
                                                                                                                                                    aGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctg  >  1:140406/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTG  >  NC_000913/1724136‑1724433

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