Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,015,703 T→C 8.0% L275S (TTG→TCG)  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,7030TC 8.0% 218.2 / 7.8 75L275S (TTG→TCG) fliGflagellar motor switch protein FliG
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base C (2/4);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GCGCAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGATGGTAATTG  >  NC_000913/2015555‑2015849
                                                                                                                                                    |                                                                                                                                                  
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:1176939/150‑1 (MQ=255)
             tcGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGc                                                                                                                                      >  2:876860/1‑150 (MQ=255)
             tcGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGc                                                                                                                                      <  1:1647332/150‑1 (MQ=255)
                 cgaGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCg                                                                                                                                  <  2:911055/150‑1 (MQ=255)
                 cgaGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCg                                                                                                                                  >  1:641509/1‑150 (MQ=255)
                   agaTGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTg                                                                                                                                <  2:1237529/150‑1 (MQ=255)
                   agaTGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTg                                                                                                                                >  1:1963983/1‑150 (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:1076202/150‑1 (MQ=255)
                                            ggATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGAtct                                                                                                       <  1:398447/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:641509/150‑1 (MQ=255)
                                                             gcaTTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCgg                                                                                      <  2:399366/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                                                                            <  2:283039/150‑1 (MQ=255)
                                                                       tCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTc                                                                            <  1:990880/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:1532389/150‑1 (MQ=255)
                                                                                                        gATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTctgct                                           <  1:1056124/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)
                                                                                                                                           agaAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGATgg        >  1:118467/1‑150 (MQ=255)
                                                                                                                                           agaAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGATgg        >  1:776132/1‑150 (MQ=255)
                                                                                                                                                 tCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGATGGTAATTg  <  1:1684608/150‑1 (MQ=255)
                                                                                                                                                 tCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGATGGTAATTg  >  1:27921/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGCAGAAAATCATCGACGAGATGTTCCTGTTCGAGAATCTGGTGGATGTCGACGATCGCAGCATTCAGCGTCTGTTGCAGGAAGTGGATTCCGAATCGCTGTTGATCGCGCTGAAAGGAGCCGAGCAGCCACTGCGCGAGAAATTCTTGCGCAATATGTCGCAGCGTGCCGCCGATATTCTGCGCGACGATCTCGCCAACCGTGGTCCGGTGCGTCTGTCGCAGGTGGAAAACGAACAGAAAGCGATTCTGCTGATTGTGCGCCGCCTTGCCGAAACTGGCGAGATGGTAATTG  >  NC_000913/2015555‑2015849

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: