Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,974,082 T→A 7.0% K523* (AAG→TAG) ‡ flgK ← flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,974,0820TA 7.0% 343.0 / 12.6 112K523* (AAG→TAG) ‡flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base T (44/60);  new base A (3/5);  total (47/65)
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 124 total aligned reads displayed.

GCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCGGTG  >  NC_002947/4973933‑4974230
                                                                                                                                                     |                                                                                                                                                     
gCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACtt                                                                                                                                                       <  1:2454773/150‑1 (MQ=255)
   gTCGGTATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                    <  2:2032740/150‑1 (MQ=255)
   gTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                    <  2:2273416/150‑1 (MQ=255)
    tCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgttg                                                                                                                                                   <  1:2880057/150‑1 (MQ=255)
    tCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGATTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgttg                                                                                                                                                   <  2:1293859/150‑1 (MQ=255)
         aTTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTc                                                                                                                                                 <  1:2238130/147‑1 (MQ=255)
         aTTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTg                                                                                                                                              >  1:928129/1‑150 (MQ=255)
         aTTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTg                                                                                                                                              >  2:1729765/1‑150 (MQ=255)
          ttGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTccatc                                                                                                                                             <  1:1935575/150‑4 (MQ=255)
          ttGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGAACGATTTTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGc                                                                                                                                             <  1:2832444/150‑1 (MQ=255)
           tGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCTCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTccatcg                                                                                                                                            <  2:1206407/150‑5 (MQ=255)
           tGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGcc                                                                                                                                            <  2:2032274/150‑1 (MQ=255)
             aaCGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCa                                                                                                                                          <  2:1680884/150‑1 (MQ=255)
              aCGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAg                                                                                                                                         <  1:102967/150‑1 (MQ=255)
                     ggTAAAGCTGTCATTGGCCGCAGTGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGcaccac                                                                                                                                  <  2:1365602/150‑1 (MQ=255)
                     ggTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGcaccac                                                                                                                                  <  1:2018060/150‑1 (MQ=255)
                      gTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACg                                                                                                                                 >  1:1258287/1‑150 (MQ=255)
                       tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGc                                                                                                                                             <  2:1934824/137‑1 (MQ=255)
                           gCTGTCATTGGCCGCAGGGCTGCCACCTATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCa                                                                                                                                          <  1:2366512/136‑1 (MQ=255)
                           gCTGTCATTGGCCGCAGGGCTGCCACCTATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCa                                                                                                                                          >  2:2366512/1‑136 (MQ=255)
                              gTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGt                                                                                                                         >  2:865556/1‑150 (MQ=255)
                               tCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                        <  2:475845/150‑1 (MQ=255)
                               tCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                        >  2:1562628/1‑150 (MQ=255)
                                 aTTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGcc                                                                                                                                            >  2:2403269/1‑128 (MQ=255)
                                  ttGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCa                                                                                                                     >  2:2056748/1‑150 (MQ=255)
                                       cGCAGGGCTGCCACCGATTGTGGTCTCGACTGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGAGAACTTGTTGTCCTGCCCAGGCACCATGCTGCGGGTACCAATCTg                                                                                                                <  2:765281/150‑1 (MQ=255)
                                       cGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGcc                                                                                                                                            >  2:1164686/1‑122 (MQ=255)
                                       cGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTccatcg                                                                                                                                            >  1:1533728/1‑118 (MQ=255)
                                       cGCAGGGCTGCCACCGATTCTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTg                                                                                                                <  2:2018001/150‑1 (MQ=255)
                                        gCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCt                                                                                                                                               >  1:1107507/1‑118 (MQ=255)
                                        gCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCt                                                                                                                                               <  2:1107507/118‑1 (MQ=255)
                                        gCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                              <  2:1560239/150‑1 (MQ=255)
                                        gCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                              >  1:239346/1‑150 (MQ=255)
                                                 ccACCTATTGTGGTATCGACGGCAAAGGTACGCGGTTTACCGCTACCATCCAGACTCGGGTCGCCGGCTGCATCCAGCATCGTCACGGCGATGGTCAACTTGTTGTCCTGTCCCGGCACCACGCTGACGGTACCAATCTGGGTACCCttt                                                                                                     <  2:969525/150‑2 (MQ=255)
                                                 ccACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                     <  1:664416/150‑1 (MQ=255)
                                                 ccACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTcc                                                                                                                                                <  2:773619/108‑1 (MQ=255)
                                                 ccACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTcc                                                                                                                                                >  1:773619/1‑108 (MQ=255)
                                                  cACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTgg                                                                                                    <  1:1690479/150‑1 (MQ=255)
                                                   aCCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGc                                                                                                   <  2:2961880/150‑1 (MQ=255)
                                                    ccGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCg                                                                                                  >  1:904068/1‑150 (MQ=255)
                                                      gATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGCCCGCATCCAGCATCGGCACGGCGACGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTc                                                                                                <  2:3023015/150‑1 (MQ=255)
                                                       aTTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCg                                                                                               <  2:53875/150‑1 (MQ=255)
                                                       aTTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCg                                                                                               >  1:1055279/1‑150 (MQ=255)
                                                              tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCt                                                                                        <  2:734894/150‑1 (MQ=255)
                                                              tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCt                                                                                        >  2:3000189/1‑150 (MQ=255)
                                                                tcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACtt                                                                                                                                                       <  2:2873104/86‑1 (MQ=255)
                                                                tcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACtt                                                                                                                                                       >  1:2873104/1‑86 (MQ=255)
                                                                  gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCg                                                                                                                           >  2:744685/1‑112 (MQ=255)
                                                                  gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCg                                                                                                                           <  1:744685/112‑1 (MQ=255)
                                                                  gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                         >  2:474272/1‑145 (MQ=255)
                                                                  gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTa                                                                                    <  1:865556/150‑1 (MQ=255)
                                                                   aCGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAg                                                                                   <  2:679850/150‑1 (MQ=255)
                                                                       cAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                     <  2:2768721/128‑1 (MQ=255)
                                                                       cAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                     >  1:2768721/1‑128 (MQ=255)
                                                                       cAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCt                                                                               <  1:3000189/150‑1 (MQ=255)
                                                                          aGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTggg                                                                            >  1:2590225/1‑150 (MQ=255)
                                                                           ggTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                         <  1:359076/136‑1 (MQ=255)
                                                                           ggTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                         >  2:359076/1‑136 (MQ=255)
                                                                           ggTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCCAACAGCTTGGAGCCCTGGGa                                                                           <  2:1442651/150‑1 (MQ=255)
                                                                                     ttACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                                 >  1:327903/1‑150 (MQ=255)
                                                                                     ttACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGCGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                                 <  2:433113/150‑1 (MQ=255)
                                                                                     ttACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATATGGGTAACCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                                 >  1:2636030/1‑150 (MQ=255)
                                                                                      tACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTgg                                                                <  2:327903/150‑1 (MQ=255)
                                                                                       aCCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGc                                                               <  2:928129/150‑1 (MQ=255)
                                                                                         cGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCt                                                                               >  2:2735799/1‑132 (MQ=255)
                                                                                           cTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGc                                                                                                                                       >  1:1017879/1‑75 (MQ=255)
                                                                                           cTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGc                                                                                                                                       <  2:1017879/75‑1 (MQ=255)
                                                                                            tACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                         >  1:1102470/1‑119 (MQ=255)
                                                                                            tACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                         <  2:1102470/119‑1 (MQ=255)
                                                                                            tACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGt                                                          <  1:2056748/150‑1 (MQ=255)
                                                                                             aCCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCAGCCCGCAGTTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTc                                                         <  2:1951899/150‑1 (MQ=255)
                                                                                              ccATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCg                                                        >  2:2599277/1‑150 (MQ=255)
                                                                                                aTCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                              >  2:2222854/1‑94 (MQ=255)
                                                                                                 tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAac                                                     <  2:2326120/150‑1 (MQ=255)
                                                                                                 tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAac                                                     >  1:85861/1‑150 (MQ=255)
                                                                                                           ggTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTc                                                                                                <  1:1241695/97‑1 (MQ=255)
                                                                                                           ggTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTc                                                                                                >  2:1241695/1‑97 (MQ=255)
                                                                                                            gTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAaca                                                    <  1:350037/140‑1 (MQ=255)
                                                                                                              cGCCGGCCGAACCCAGTATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCCGGGCGGGCGTGCCGATGTGTAGGCGCGGTACCCTCATGGTCGAACAGCTTGTAGAATTGGGAGCCACCGGTGGCAGCGTCGAACACCACGCGTACCg                                        <  1:700859/150‑1 (MQ=255)
                                                                                                                     cgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTa                                           >  2:607592/1‑140 (MQ=255)
                                                                                                                     cgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTa                                           <  1:607592/140‑1 (MQ=255)
                                                                                                                      gCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaa                                <  1:1562628/150‑1 (MQ=255)
                                                                                                                         tCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcg                             >  2:2414590/1‑150 (MQ=255)
                                                                                                                          ccAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGaa                                                      <  2:2866194/124‑1 (MQ=255)
                                                                                                                          ccAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGaa                                                      >  1:2866194/1‑124 (MQ=255)
                                                                                                                                      cGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActgct                <  2:56197/150‑1 (MQ=255)
                                                                                                                                       ggCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActgctg               >  1:989620/1‑150 (MQ=255)
                                                                                                                                           aTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATc           >  1:1245316/1‑150 (MQ=255)
                                                                                                                                             ggACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                              <  1:503680/49‑1 (MQ=255)
                                                                                                                                             ggACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActg                  >  2:2368101/1‑141 (MQ=255)
                                                                                                                                             ggACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActg                  <  1:2368101/141‑1 (MQ=255)
                                                                                                                                              gACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTg                                                                              <  2:283650/80‑1 (MQ=255)
                                                                                                                                                cAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGt        >  2:881180/1‑148 (MQ=255)
                                                                                                                                                 aaCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCg     >  1:1674411/1‑150 (MQ=255)
                                                                                                                                                  aCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTg       >  1:2592226/1‑147 (MQ=255)
                                                                                                                                                  aCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTg       <  2:2592226/147‑1 (MQ=255)
                                                                                                                                                   cTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGc                                                            >  2:1999316/1‑93 (MQ=255)
                                                                                                                                                   cTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGc                                                            <  1:1999316/93‑1 (MQ=255)
                                                                                                                                                    ttgttgTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatc                              >  2:1483255/1‑122 (MQ=255)
                                                                                                                                                    ttgttgTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCGgtg  >  1:1216616/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCGGTG  >  NC_002947/4973933‑4974230

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