Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,285 G→C 7.8% 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 7.8% 338.7 / 15.5 103A124G (GCT→GGT) sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (47/48);  new base C (5/3);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 106 total aligned reads displayed.

CTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGG  >  NC_000913/1724137‑1724434
                                                                                                                                                    |                                                                                                                                                     
cttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGc                                                                                                                                                      >  1:397045/1‑150 (MQ=255)
cttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGc                                                                                                                                                      >  2:1289617/1‑150 (MQ=255)
cttaattacACCACACGCATAGCGTTCACCGCCACCGCCCAGTGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGc                                                                                                                                                      <  2:87299/150‑1 (MQ=255)
 ttaattacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCt                                                                                                                                                     <  1:1161920/150‑1 (MQ=255)
      tacACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGcc                                                                                                                                                >  1:637042/1‑150 (MQ=255)
       acacCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCg                                                                                                                                               >  2:44654/1‑150 (MQ=255)
           cacaGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTGATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCat                                                                                                                                           >  1:829995/1‑150 (MQ=255)
             caGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAAAGCTACCGATGCCGTCatcg                                                                                                                                         <  2:1552860/150‑3 (MQ=255)
               ggCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAAAGCTACCGATGCCGTCatcgcg                                                                                                                                       <  1:1230422/150‑5 (MQ=255)
                gCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTg                                                                                                                                      <  2:442629/150‑1 (MQ=255)
                 cATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCAAAGCTACCGATGCCGTCatcgcgcc                                                                                                                                     <  2:1003934/150‑7 (MQ=255)
                   aaGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgacg                                                                                                                                   <  2:1318273/149‑1 (MQ=255)
                     gCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgacgac                                                                                                                                 >  1:965222/1‑150 (MQ=255)
                     gCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATAGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgacgac                                                                                                                                 >  1:11863/1‑150 (MQ=255)
                      cGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGAcc                                                                                                                                <  1:1637757/150‑1 (MQ=255)
                        ttCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAg                                                                                                                              <  1:1526393/150‑1 (MQ=255)
                          caccgccaccgcCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTg                                                                                                                            >  1:626149/1‑150 (MQ=255)
                          caccgccaccgcCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTg                                                                                                                            >  2:1090859/1‑150 (MQ=255)
                           accgccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgc                                                                                                                           >  2:1499092/1‑150 (MQ=255)
                               ccaccgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggc                                                                                                                       >  2:1224140/1‑150 (MQ=255)
                                caacgccCAGCGGTTTAGGTTGATCGTACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCTGTGATTTCAGACGAGGAGCGATGACGGAATCGGTAGCTTTGCCGTCATTATTGACGAAGAGTgacggca                                                                                                                      >  2:1106560/4‑150 (MQ=255)
                                   cgccCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGAt                                                                                                                   <  2:818543/150‑1 (MQ=255)
                                     cccAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCg                                                                                                                 >  2:1719154/1‑150 (MQ=255)
                                       cAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGcc                                                                                                               <  1:790139/150‑1 (MQ=255)
                                           ggTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTaaa                                                                                                           <  2:1286961/150‑1 (MQ=255)
                                           ggTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTaaa                                                                                                           >  1:1429536/1‑150 (MQ=255)
                                                 ggTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCg                                                                                                     >  2:917551/1‑150 (MQ=255)
                                                 ggTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCg                                                                                                     >  1:1114281/1‑150 (MQ=255)
                                                 ggTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCg                                                                                                     >  1:1206210/1‑150 (MQ=255)
                                                   ttGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCTTCCAGTGATTTCAGACGAGGCGCGATGACGTCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGc                                                                                                   <  2:733438/150‑1 (MQ=255)
                                                          gACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCt                                                                                            >  2:1038100/1‑150 (MQ=255)
                                                          gACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCGTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCt                                                                                            <  1:688764/150‑1 (MQ=255)
                                                            cATATTATCGGCGACAACGTGGACCATCAGCGCTATGTCTTTGATTTCCTCCAGTGAATTCAGACGAGGCGCGATGACGGCATTGGTAGCTTTGCCGTCATTATTGACGACCAGTCCAGGCAGATCGCCTAAATGCCCGGCACCTTCTgg                                                                                          <  1:1120482/150‑1 (MQ=255)
                                                            cATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCgg                                                                                          <  2:455829/150‑1 (MQ=255)
                                                               attatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGccc                                                                                       <  2:1472309/150‑1 (MQ=255)
                                                                 tatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCtt                                                                                     <  2:618251/150‑1 (MQ=255)
                                                                  atCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGAGCGCCTAAATGCCCGGCACCTTCTGGCCCTTc                                                                                    <  1:352643/150‑1 (MQ=255)
                                                                     gccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATg                                                                                 <  1:1262572/150‑1 (MQ=255)
                                                                     gccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGACGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGACCGGCACCTTCTGGCCCTTCATg                                                                                 >  2:811247/1‑150 (MQ=255)
                                                                      ccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGt                                                                                <  1:730395/150‑1 (MQ=255)
                                                                      ccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGAGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGt                                                                                <  2:18054/150‑1 (MQ=255)
                                                                      ccgccAACGGGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGt                                                                                <  1:1161572/150‑1 (MQ=255)
                                                                       cgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCAAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGtt                                                                               >  1:802568/1‑150 (MQ=255)
                                                                        gccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGttt                                                                              <  1:975451/150‑1 (MQ=255)
                                                                         ccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTa                                                                             >  1:1209415/1‑150 (MQ=255)
                                                                         ccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTa                                                                             >  2:521396/1‑150 (MQ=255)
                                                                          cAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTAATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTAc                                                                            >  1:899625/1‑150 (MQ=255)
                                                                           aaCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCTCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTAcc                                                                           <  2:1167249/150‑1 (MQ=255)
                                                                            aCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCg                                                                          >  1:835853/1‑150 (MQ=255)
                                                                             cGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCgg                                                                         <  2:1425697/150‑1 (MQ=255)
                                                                              gTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGt                                                                        <  1:61388/150‑1 (MQ=255)
                                                                                   ccATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACATTCTGGCCCTTCATGTTTACCGGTAtttt                                                                   >  1:254508/1‑150 (MQ=255)
                                                                                    cATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTtg                                                                  <  1:1289617/150‑1 (MQ=255)
                                                                                     aTCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTAACGGTATTTtgt                                                                 >  1:935105/1‑150 (MQ=255)
                                                                                         gcgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCAGGCCCTTCATGTTTACAGGTATATTGTGGAt                                                             >  2:288555/1‑150 (MQ=255)
                                                                                           gcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCTTTATTGACGCCCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCTTGTTTACCGGTATTTTGTGGATCa                                                           >  1:91545/1‑150 (MQ=255)
                                                                                           gcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCa                                                           >  2:1707953/1‑150 (MQ=255)
                                                                                            cTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCaa                                                          >  2:463158/1‑150 (MQ=255)
                                                                                             tttGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAg                                                         <  1:612406/150‑1 (MQ=255)
                                                                                              ttGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGa                                                        >  2:774925/1‑150 (MQ=255)
                                                                                                gTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCCCCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGCTCAAGATg                                                      >  2:292456/1‑150 (MQ=255)
                                                                                                       aTTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCt                                               >  1:749719/1‑150 (MQ=255)
                                                                                                          tCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCg                                            <  2:1612453/150‑1 (MQ=255)
                                                                                                           cATCCAGTGATTTCAGCCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGGGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCgg                                           <  1:1171494/150‑1 (MQ=255)
                                                                                                           cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGa                                                        <  2:88132/137‑1 (MQ=255)
                                                                                                           cATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGa                                                        >  1:88132/1‑137 (MQ=255)
                                                                                                            aTCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCGTGCATGTTTAGCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGa                                          >  1:1057489/1‑150 (MQ=255)
                                                                                                             tCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTTTTGACGACCAGTGCAGGCAGATCGCCTAAATGGCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGAt                                         <  2:220564/150‑1 (MQ=255)
                                                                                                              ccAGTGATTTCAGACGAGGCGCGATGACGGCAACGGTAGCTTTGCCTTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCACGGCACCTTCAGGCCCTTCATGATTACCGGTATATTTTGGATCAAGATGCCCGCCTGCGGAtt                                        >  1:312546/1‑150 (MQ=255)
                                                                                                               cAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTc                                       >  1:1541764/1‑150 (MQ=255)
                                                                                                                aGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTcc                                      >  1:1652706/1‑150 (MQ=255)
                                                                                                                aGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTcc                                      <  2:626149/150‑1 (MQ=255)
                                                                                                                        cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCt                              <  1:1090859/150‑1 (MQ=255)
                                                                                                                        cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCt                              >  2:567558/1‑150 (MQ=255)
                                                                                                                        cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCt                              >  1:1723905/1‑150 (MQ=255)
                                                                                                                        cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCGCGGAATCcgcaggcggg                             <  1:994686/150‑8 (MQ=255)
                                                                                                                          gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTgg                            >  2:529876/1‑150 (MQ=255)
                                                                                                                               ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTg                       <  2:1711782/150‑1 (MQ=255)
                                                                                                                               ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTg                       >  1:985945/1‑150 (MQ=255)
                                                                                                                                gcgcGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                      <  1:1294413/150‑1 (MQ=255)
                                                                                                                                gcgcGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                      <  1:1470282/150‑1 (MQ=255)
                                                                                                                                gcgcGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCGTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                      <  1:142609/150‑1 (MQ=255)
                                                                                                                                 cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTGTGcc                     >  1:1430541/1‑150 (MQ=255)
                                                                                                                                   cgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCAt                   <  1:1018117/150‑1 (MQ=255)
                                                                                                                                    aaTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTAGTTTGTGGATCAAGATGCCCGCCTGCGGAGTCCGCGGCGCTGGCTTTGCCAtg                  >  1:1424823/2‑149 (MQ=255)
                                                                                                                                     aTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCt                 >  2:860295/1‑150 (MQ=255)
                                                                                                                                     aTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCt                 >  1:1634116/1‑150 (MQ=255)
                                                                                                                                     aTGACGGCAAAGCTACCGATGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCt                 >  1:622666/1‑150 (MQ=255)
                                                                                                                                      tGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCtt                <  1:620086/150‑1 (MQ=255)
                                                                                                                                      tGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCtt                <  1:1269564/150‑1 (MQ=255)
                                                                                                                                             aTCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggc         <  2:637042/150‑1 (MQ=255)
                                                                                                                                              tCGGTAGCTTTGCCGTCATTATTGACTACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggct        <  1:951339/150‑1 (MQ=255)
                                                                                                                                               gggTAGCTTTGCCGTCATTATTGACGACCAGTGCATGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctg       <  2:1373747/149‑1 (MQ=255)
                                                                                                                                               cGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctg       <  2:125553/150‑1 (MQ=255)
                                                                                                                                                ggTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCGATCTTTGGtggctgg      >  1:1669915/1‑150 (MQ=255)
                                                                                                                                                ggTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctgg      >  1:398898/1‑150 (MQ=255)
                                                                                                                                                 gTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggc     <  1:1547558/150‑1 (MQ=255)
                                                                                                                                                    gCTTTGCCGTCATTATTGACGACCTGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctgg  <  2:899625/150‑1 (MQ=255)
                                                                                                                                                    gCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCTCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctgg  <  2:1620089/150‑1 (MQ=255)
                                                                                                                                                    gCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCGTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGCTTCCGCGGCGCTGGCTTTGCGATCTTTGGtggctggctgg  >  2:72048/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGG  >  NC_000913/1724137‑1724434

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: