Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,795,164 A→T 10.6% intergenic (+38/+30) ynfM → / ← PP_5073 efflux transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,795,1640AT10.6% 756.7 / 56.9 249intergenic (+38/+30)ynfM/PP_5073efflux transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (101/121);  new base T (11/16);  total (112/137)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

TTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAGGTCCTCA  >  NC_002947/5795016‑5795311
                                                                                                                                                       |                                                                                                                                                   
  cTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCg                                                                                                                                                  <  2:6218094/150‑1 (MQ=39)
  cTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCg                                                                                                                                                  <  1:1915855/150‑1 (MQ=255)
  cTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACACCCGCTCCCACAGGTTGTGGGTGCg                                                                                                                                                  <  2:3120966/150‑1 (MQ=255)
   tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                                 >  1:3480001/1‑150 (MQ=255)
   tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                                 >  2:1669861/1‑150 (MQ=255)
   tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                                 <  2:4975704/150‑1 (MQ=255)
   tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                                 <  1:2720904/150‑1 (MQ=255)
   tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                                 >  2:4121820/1‑150 (MQ=255)
   tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                                 >  1:2256480/1‑150 (MQ=255)
   tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATGGGCAGCGGGCGGGCGGTGGCCAGGACGGTGGCGTTGCACCTGAGCAAGTTGCCCCCCAAACCGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                                 <  2:3873536/150‑1 (MQ=255)
    ggCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCggg                                                                                                                                                >  1:1000752/1‑150 (MQ=255)
     gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgt                                                                                                                                               >  1:983272/1‑150 (MQ=255)
     gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgt                                                                                                                                               <  1:4037121/150‑1 (MQ=255)
     gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGAGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCgggc                                                                                                                                               <  2:6030176/150‑2 (MQ=38)
       aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                             <  1:2245569/150‑1 (MQ=255)
       aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                             <  1:6182911/150‑1 (MQ=255)
       aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                             <  2:4962762/150‑1 (MQ=255)
          aGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGcgtgcc                                                                                                                                          <  1:4622556/150‑3 (MQ=37)
          aGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGAACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGAc                                                                                                                                          <  2:1999975/150‑1 (MQ=255)
           gTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGAcc                                                                                                                                         >  2:3412292/1‑150 (MQ=255)
            tGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCa                                                                                                                                        <  1:1081996/150‑1 (MQ=255)
            tGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCa                                                                                                                                        <  2:6242923/150‑1 (MQ=255)
            tGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGAGGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCa                                                                                                                                        <  2:4024300/150‑1 (MQ=255)
              gggCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACg                                                                                                                                      <  2:5030165/150‑1 (MQ=255)
                 cTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGaat                                                                                                                                   <  2:1406438/150‑2 (MQ=37)
                  tGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAaga                                                                                                                                  <  2:5824596/150‑1 (MQ=255)
                   ggAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTCTGATGGCCTCTTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGaa                                                                                                                                    <  1:6749758/150‑1 (MQ=18)
                     aaCGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGaatgcat                                                                                                                               <  2:3497571/150‑6 (MQ=25)
                       cGGGATCGGGCGTTTCATTGGCAGCTGGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGCTG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCa                                                                                                                             <  2:186352/150‑1 (MQ=255)
                           aTCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcag                                                                                                                         >  2:1088431/1‑150 (MQ=255)
                               ggCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCACGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatca                                                                                                                     <  1:5471930/150‑1 (MQ=255)
                               ggCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatca                                                                                                                     >  1:2378238/1‑150 (MQ=255)
                                 cTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGt                                                                                                                   >  1:5851935/1‑150 (MQ=255)
                                    ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtga                                                                                                                <  1:5070378/150‑1 (MQ=255)
                                    ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtga                                                                                                                <  2:6814162/150‑1 (MQ=255)
                                        ttGGCAGCTTGCAGGCGGTGGCCTTGATGGTGGCCTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCAGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTgatgat                                                                                                            >  2:301991/1‑150 (MQ=255)
                                          ggCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                                          <  1:3304885/150‑1 (MQ=255)
                                           gCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtctc                                                                                                         >  1:5281631/1‑150 (MQ=255)
                                              gCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGAc                                                                                                      <  2:3751440/150‑1 (MQ=255)
                                                  gCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAata                                                                                                  <  2:3182981/150‑1 (MQ=255)
                                                  gCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAata                                                                                                  <  2:2289401/150‑1 (MQ=255)
                                                    tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatc                                                                                                >  1:4215629/1‑150 (MQ=255)
                                                    tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatc                                                                                                >  1:4113573/1‑150 (MQ=255)
                                                     ggCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcc                                                                                               <  1:6333565/150‑1 (MQ=255)
                                                      gCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcca                                                                                              <  2:4726824/150‑1 (MQ=255)
                                                        ggTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccaca                                                                                            >  2:5052182/1‑150 (MQ=255)
                                                         gTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacat                                                                                           <  1:2142550/150‑1 (MQ=255)
                                                             ccTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacatccac                                                                                       >  2:2188756/1‑150 (MQ=255)
                                                              cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                                                      >  1:2881158/1‑150 (MQ=255)
                                                              cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                                                      >  2:6675926/1‑150 (MQ=255)
                                                              cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                                                      >  1:6004747/1‑150 (MQ=255)
                                                               ttGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACtt                                                                                     >  2:1282748/1‑150 (MQ=255)
                                                               ttGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAACAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACtt                                                                                     >  1:677465/1‑150 (MQ=255)
                                                                    ggtggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCgg                                                                                <  2:3005583/150‑1 (MQ=255)
                                                                    ggtggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCgg                                                                                <  1:255927/150‑1 (MQ=255)
                                                                     gtggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGc                                                                               >  1:5567628/1‑150 (MQ=255)
                                                                     gtggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGc                                                                               <  1:5109011/150‑1 (MQ=255)
                                                                          gTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                                          <  1:252941/150‑1 (MQ=255)
                                                                           ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatca                                                                         <  1:6464017/150‑1 (MQ=255)
                                                                           ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatca                                                                         <  1:1457541/150‑1 (MQ=255)
                                                                                ccTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGc                                                                    <  1:4362693/150‑1 (MQ=255)
                                                                                 cTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCTCCCACAGGTTGTGGGTGCgg                                                                                                                                                 >  1:454563/1‑72 (MQ=21)
                                                                                   gAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCttt                                                                 <  1:5426978/150‑1 (MQ=255)
                                                                                   gAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCttt                                                                 <  2:573977/150‑1 (MQ=255)
                                                                                   gAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCttt                                                                 <  1:4987899/150‑1 (MQ=255)
                                                                                   gAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCttt                                                                 >  2:915683/1‑150 (MQ=255)
                                                                                       aaGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCaa                                                             >  1:1861189/1‑150 (MQ=255)
                                                                                        aGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAc                                                            <  2:6113395/150‑1 (MQ=255)
                                                                                        aGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAc                                                            <  2:387231/150‑1 (MQ=255)
                                                                                        aGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAc                                                            <  1:2297739/150‑1 (MQ=255)
                                                                                         gTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAcc                                                           >  2:2880867/1‑150 (MQ=255)
                                                                                              cACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAcc                                                      >  1:5822717/1‑150 (MQ=255)
                                                                                                  cAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTc                                                  >  1:529825/1‑150 (MQ=255)
                                                                                                        ggTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTc                                            <  2:4884960/150‑1 (MQ=255)
                                                                                                          tgtgATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCAc                                          >  1:6666118/1‑150 (MQ=255)
                                                                                                           gtgATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCAcc                                         <  1:3647178/150‑1 (MQ=255)
                                                                                                             gATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCtt                                       <  2:3285129/150‑1 (MQ=255)
                                                                                                              aTG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                                      <  1:307991/150‑1 (MQ=255)
                                                                                                              aTG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                                      <  1:3729486/150‑1 (MQ=255)
                                                                                                              aTG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTCTTCAACCTCACCGGTCAGCTTCACCTTg                                      <  2:5440028/150‑1 (MQ=255)
                                                                                                                   ctcttCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCAACAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGt                                    >  1:4394501/4‑150 (MQ=25)
                                                                                                                     aTTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCt                                  >  1:1235051/1‑150 (MQ=255)
                                                                                                                      ttCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCtt                                 <  2:4548743/150‑1 (MQ=255)
                                                                                                                      ttCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCtt                                 <  1:731696/150‑1 (MQ=255)
                                                                                                                          tGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCg                             >  2:446374/1‑150 (MQ=37)
                                                                                                                          tGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCg                             <  1:1810688/150‑1 (MQ=255)
                                                                                                                             gCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTg                          >  1:6139309/1‑150 (MQ=255)
                                                                                                                               acacccGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGaa                        >  2:917568/4‑150 (MQ=38)
                                                                                                                                 acccGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACg                      >  2:5618501/2‑150 (MQ=38)
                                                                                                                                 acccGCTCCCACAGGTTGTGGGTGCGGGCGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACg                      >  1:2153734/2‑150 (MQ=39)
                                                                                                                                   ccGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGgt                    >  2:1745601/1‑150 (MQ=255)
                                                                                                                                       aCCCACAACCTGTGGGAGCGGGAGTGACCACGAAGAGGCCATCAGATCAGTTGATGCTCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTcg                <  1:3785328/150‑1 (MQ=255)
                                                                                                                                         ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTcggc              >  1:4431511/1‑150 (MQ=255)
                                                                                                                                          cacaggttgtgGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAACTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTcggcg             >  1:6401577/8‑150 (MQ=255)
                                                                                                                                          cacaACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTcggcg             >  1:1148937/1‑150 (MQ=255)
                                                                                                                                             aaCCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCa          <  2:2640963/150‑1 (MQ=255)
                                                                                                                                             aaCCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAAGATCCACATCCACGTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCa          <  1:6463588/150‑1 (MQ=255)
                                                                                                                                                 tgtgGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAGGTc      >  1:1978140/1‑150 (MQ=255)
                                                                                                                                                    gggAGCGGGTGTGACCAAGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAGGTCCTc   <  1:762197/150‑1 (MQ=255)
                                                                                                                                                     ggAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAGGTCCTCa  <  1:4877332/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAGGTCCTCA  >  NC_002947/5795016‑5795311

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