Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,015,693 G→A 8.3% E272K (GAG→AAG)  fliG → flagellar motor switch protein FliG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,015,6930GA 8.3% 213.0 / 7.0 72E272K (GAG→AAG) fliGflagellar motor switch protein FliG
Reads supporting (aligned to +/- strand):  ref base G (29/37);  new base A (2/4);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GGCGAGCTGGCGCAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGA  >  NC_000913/2015546‑2015840
                                                                                                                                                   |                                                                                                                                                   
ggCGAGCTGGCGCAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCgag                                                                                                                                                   <  1:1287673/150‑1 (MQ=255)
        ggCGCAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGaa                                                                                                                                                 >  2:1457567/1‑144 (MQ=255)
        ggCGCAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGaa                                                                                                                                                 <  1:1457567/144‑1 (MQ=255)
         tcccAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTg                                                                                                                                          <  1:565962/147‑1 (MQ=255)
            cAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTgcgc                                                                                                                                       <  1:833336/150‑1 (MQ=255)
            cAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTgcgc                                                                                                                                       <  1:490353/150‑1 (MQ=255)
            cAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCAAGAATTTCTCgcgc                                                                                                                                       <  1:22396/150‑1 (MQ=255)
              gAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCAAGAATTTCTCGCGCag                                                                                                                                     >  1:1635404/1‑149 (MQ=255)
                 aaTCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCAAGAATTTCTCGCGCagtgg                                                                                                                                  <  2:1299557/150‑5 (MQ=255)
                 aaTCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCAAGAATTTCTCGCGCagtgg                                                                                                                                  <  1:1646137/150‑5 (MQ=255)
                  atcatcGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCAAGAATTTCTCGCGCagtggc                                                                                                                                 <  2:1635404/150‑6 (MQ=255)
                      tcGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGc                                                                                                                             >  2:876860/1‑150 (MQ=255)
                      tcGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGc                                                                                                                             <  2:1176939/150‑1 (MQ=255)
                      tcGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGc                                                                                                                             <  1:1647332/150‑1 (MQ=255)
                          cgaGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCg                                                                                                                         <  2:911055/150‑1 (MQ=255)
                          cgaGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCg                                                                                                                         >  1:641509/1‑150 (MQ=255)
                            agaTGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTg                                                                                                                       >  1:1963983/1‑150 (MQ=255)
                            agaTGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTg                                                                                                                       <  2:1237529/150‑1 (MQ=255)
                                      gTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGatat                                                                                                             >  2:1687860/1‑150 (MQ=255)
                                        tCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTc                                                                                                           >  2:1532389/1‑150 (MQ=255)
                                        tCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTc                                                                                                           >  1:1349138/1‑150 (MQ=255)
                                           agaATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTgc                                                                                                        >  1:353982/1‑150 (MQ=255)
                                                cTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGcgac                                                                                                   >  2:735132/1‑150 (MQ=255)
                                                     ggATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGAtct                                                                                              <  1:398447/150‑1 (MQ=255)
                                                     ggATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGAtct                                                                                              <  1:1076202/150‑1 (MQ=255)
                                                        tGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGc                                                                                           >  1:282185/1‑150 (MQ=255)
                                                               gaTCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGt                                                                                    <  2:751111/150‑1 (MQ=255)
                                                                   gcagcaTTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTc                                                                                >  1:1404127/1‑150 (MQ=255)
                                                                      gcaTTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCgg                                                                             <  2:399366/150‑1 (MQ=255)
                                                                      gcaTTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCgg                                                                             <  2:641509/150‑1 (MQ=255)
                                                                        aTTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTg                                                                           >  1:58353/1‑150 (MQ=255)
                                                                          tCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCg                                                                         <  2:749210/150‑1 (MQ=255)
                                                                           cAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGt                                                                        >  2:160957/1‑150 (MQ=255)
                                                                            aGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTc                                                                       <  1:977821/150‑1 (MQ=255)
                                                                             gCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCt                                                                      >  2:1635655/1‑150 (MQ=255)
                                                                             ccGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCt                                                                      >  1:1420531/2‑150 (MQ=255)
                                                                               gTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGt                                                                    >  2:1019183/1‑150 (MQ=255)
                                                                                tCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTc                                                                   <  1:990880/150‑1 (MQ=255)
                                                                                tCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTc                                                                   <  2:283039/150‑1 (MQ=255)
                                                                                  tGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGc                                                                 <  2:266245/150‑1 (MQ=255)
                                                                                   gTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCa                                                                >  2:1980137/1‑150 (MQ=255)
                                                                                      gCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGt                                                             >  2:1684608/1‑150 (MQ=255)
                                                                                      gCAGGAAGTGGATTCCGAATCGCTGTTGATATCGCTGAAAGGAGCAGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTCGTCCGTTGCGTCTGACGCAGGt                                                             >  2:1871488/1‑150 (MQ=255)
                                                                                             gTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAAc                                                      <  1:1113808/150‑1 (MQ=255)
                                                                                              tGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACg                                                     <  1:1033958/150‑1 (MQ=255)
                                                                                               ggATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGa                                                    >  2:25366/1‑150 (MQ=255)
                                                                                               ggATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGa                                                    <  1:876860/150‑1 (MQ=255)
                                                                                                  ttCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGATAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATTTTCTGCGCGACTATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACa                                                 <  2:1464511/150‑1 (MQ=255)
                                                                                                        aTCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGc                                           <  2:1963983/150‑1 (MQ=255)
                                                                                                          cGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGa                                         <  2:972961/150‑1 (MQ=255)
                                                                                                            cTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGAtt                                       >  1:240271/1‑150 (MQ=255)
                                                                                                              gTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTct                                     >  2:1175727/1‑150 (MQ=255)
                                                                                                                 gATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTTGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTctgct                                  <  2:1420531/150‑1 (MQ=255)
                                                                                                                 gATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTctgct                                  <  1:1056124/150‑1 (MQ=255)
                                                                                                                 gATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTctgct                                  <  1:1532389/150‑1 (MQ=255)
                                                                                                                         tGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTgc                          <  1:1019183/150‑1 (MQ=255)
                                                                                                                         tGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTATCCTGATTGTgc                          >  2:787472/1‑150 (MQ=255)
                                                                                                                          gAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTgcg                         >  1:1277473/1‑150 (MQ=255)
                                                                                                                               gAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAAc                                                      >  2:157340/1‑116 (MQ=255)
                                                                                                                               gAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAAc                                                      <  1:157340/116‑1 (MQ=255)
                                                                                                                               gAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCgccgcc                    >  2:877269/1‑150 (MQ=255)
                                                                                                                                aGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCt                   >  2:922148/1‑150 (MQ=255)
                                                                                                                                  ccGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGcgccg                      <  1:299112/145‑1 (MQ=255)
                                                                                                                                  ccGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGcgccg                      >  2:299112/1‑145 (MQ=255)
                                                                                                                                  ccGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTg                 <  1:1980137/150‑1 (MQ=255)
                                                                                                                                    gAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGcc               <  2:58353/150‑1 (MQ=255)
                                                                                                                                     agcagcCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCg              <  1:758617/150‑1 (MQ=255)
                                                                                                                                          ccACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACt         <  1:922148/150‑1 (MQ=255)
                                                                                                                                           cACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTg        >  1:1986437/1‑150 (MQ=255)
                                                                                                                                             cTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGc      <  1:735132/150‑1 (MQ=255)
                                                                                                                                             cTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGGGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGc      <  2:1258533/150‑1 (MQ=255)
                                                                                                                                              tGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCg     <  2:282185/150‑1 (MQ=255)
                                                                                                                                               gcgcAAGAATTTCTCGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCga    >  1:1054742/1‑150 (MQ=255)
                                                                                                                                                 gcgAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCgaga  <  2:1282337/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCGAGCTGGCGCAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGA  >  NC_000913/2015546‑2015840

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