Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,541,800 A→T 5.6% D149V (GAC→GTC)  PP_3125 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,541,8000AT 5.6% 863.2 / 26.6 251D149V (GAC→GTC) PP_3125Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (133/104);  new base T (9/5);  total (142/109)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 253 total aligned reads displayed.

GCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAACGCC  >  NC_002947/3541651‑3541948
                                                                                                                                                     |                                                                                                                                                    
 cgcAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGAc                                                                                                                                                     >  2:3256164/1‑150 (MQ=255)
 cgcAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGAc                                                                                                                                                     <  2:1558546/150‑1 (MQ=255)
    aGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGcc                                                                                                                                                  >  2:5786963/1‑150 (MQ=255)
      ccTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACTCCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCt                                                                                                                                                <  2:1289942/150‑1 (MQ=255)
      ccTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCt                                                                                                                                                >  2:6026950/1‑150 (MQ=255)
        tgctgcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTGc                                                                                                                                              <  1:1113803/150‑1 (MQ=255)
        tgctgcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGc                                                                                                                                              <  2:763113/150‑1 (MQ=255)
         gctgcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTGACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCa                                                                                                                                             <  2:6198326/150‑1 (MQ=255)
         gctgcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCa                                                                                                                                             <  1:4279438/150‑1 (MQ=255)
          ctgcCGTTTCTGAGTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTGctc                                                                                                                                            <  2:2944818/150‑3 (MQ=255)
           tgcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCAtt                                                                                                                                           <  1:1570928/150‑1 (MQ=255)
            gcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCAttt                                                                                                                                          <  1:2013037/150‑1 (MQ=255)
             ccGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTc                                                                                                                                         >  1:6036259/1‑150 (MQ=255)
               gTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCaa                                                                                                                                       >  2:874251/1‑150 (MQ=255)
                tttCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAAc                                                                                                                                      <  1:4438086/150‑1 (MQ=255)
                    tGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGccc                                                                                                                                  >  1:4827072/1‑150 (MQ=255)
                     gATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCa                                                                                                                                 <  1:5500616/150‑1 (MQ=255)
                     gATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCa                                                                                                                                 <  1:6060734/150‑1 (MQ=255)
                                 gTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGAcaccgcaccg                                                                                                                     >  1:6318745/1‑150 (MQ=255)
                                  tCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCAGGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCgg                                                                                                                    >  1:110473/1‑150 (MQ=255)
                                    cGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCt                                                                                                                  <  1:3789357/150‑1 (MQ=255)
                                    cGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCACGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCt                                                                                                                  <  2:5696416/150‑1 (MQ=255)
                                     gAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTg                                                                                                                 >  1:5339895/1‑150 (MQ=255)
                                      aGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTgc                                                                                                                <  2:6436844/150‑1 (MQ=255)
                                        ttCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTgcgc                                                                                                              >  1:4368803/1‑150 (MQ=255)
                                           aGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCg                                                                                                           >  2:5091386/1‑150 (MQ=255)
                                           aGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCg                                                                                                           <  1:2986714/150‑1 (MQ=255)
                                            gCGACCCTACGTTCAAGGAACACCAGGGCGACGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCgg                                                                                                          >  2:6395700/1‑150 (MQ=255)
                                              gACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTa                                                                                                        <  1:2312725/150‑1 (MQ=255)
                                                 ccTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgt                                                                                                     <  1:6444730/150‑1 (MQ=255)
                                                 ccTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgt                                                                                                     >  1:2681065/1‑150 (MQ=255)
                                                 ccTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgt                                                                                                     >  2:4177921/1‑150 (MQ=255)
                                                  cTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgtg                                                                                                    >  2:4930815/1‑150 (MQ=255)
                                                   tACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgtgt                                                                                                   >  1:5571470/1‑150 (MQ=255)
                                                   tACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgtgt                                                                                                   <  2:1246006/150‑1 (MQ=255)
                                                   tACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgtgt                                                                                                   <  1:2881334/150‑1 (MQ=255)
                                                   aaCGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgtgt                                                                                                   <  2:1827526/149‑1 (MQ=255)
                                                     cGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGa                                                                                                 <  2:4855928/150‑1 (MQ=255)
                                                     cGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGa                                                                                                 <  1:6270100/150‑1 (MQ=255)
                                                      gTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGAc                                                                                                >  2:2894780/1‑150 (MQ=255)
                                                       ttCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACg                                                                                               >  2:5167554/1‑150 (MQ=255)
                                                        tCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGc                                                                                              >  1:347530/1‑150 (MQ=255)
                                                           aGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCgg                                                                                           >  1:2206375/1‑150 (MQ=255)
                                                            ggAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGgc                                                                                          >  2:5285878/1‑150 (MQ=255)
                                                             gAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGgcg                                                                                         >  1:1152860/1‑150 (MQ=255)
                                                               acacCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCa                                                                                       >  1:6000386/1‑150 (MQ=255)
                                                                cacCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAg                                                                                      <  2:3007419/150‑1 (MQ=255)
                                                                 acCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGt                                                                                     <  1:1382817/150‑1 (MQ=255)
                                                                  ccAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGtt                                                                                    >  2:6021769/1‑150 (MQ=255)
                                                                   cAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTa                                                                                   >  1:1535603/1‑150 (MQ=255)
                                                                    aGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTAc                                                                                  <  1:6653908/150‑1 (MQ=255)
                                                                      ggCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTg                                                                                <  2:5533010/150‑1 (MQ=255)
                                                                        cgtcgagTTCCTGTTTGTTCACACCGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACTCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgt                                                                              <  2:2237755/146‑1 (MQ=255)
                                                                        cGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgt                                                                              >  1:3918781/1‑150 (MQ=255)
                                                                         gAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgtt                                                                             <  2:739871/150‑1 (MQ=255)
                                                                         gAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgtt                                                                             <  2:2064472/150‑1 (MQ=255)
                                                                            gagTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgtt                                                                          <  1:726740/150‑1 (MQ=255)
                                                                             agTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttg                                                                         <  2:5326991/150‑1 (MQ=255)
                                                                             agTTCCTGTCTGTTCACACAGCCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttg                                                                         <  1:6501159/150‑1 (MQ=255)
                                                                               ttCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttgtt                                                                       >  1:1059775/1‑150 (MQ=255)
                                                                                 ccTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTca                                                                     <  2:3814186/150‑1 (MQ=255)
                                                                                 ccTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTca                                                                     <  1:3794225/150‑1 (MQ=255)
                                                                                  cTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcac                                                                    <  2:6219940/150‑1 (MQ=255)
                                                                                   tgtttgttCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcacc                                                                   >  1:6695243/1‑150 (MQ=255)
                                                                                     tttgttCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcaccac                                                                 >  2:4975821/1‑150 (MQ=255)
                                                                                          tCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGa                                                            <  1:2481797/150‑1 (MQ=255)
                                                                                            acacaGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACg                                                          >  2:2169399/1‑150 (MQ=255)
                                                                                            acacaGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACg                                                          >  1:5943530/1‑150 (MQ=255)
                                                                                              acaGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGaa                                                        >  1:4801396/1‑150 (MQ=255)
                                                                                              acaGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGaa                                                        <  1:100539/150‑1 (MQ=255)
                                                                                              acaGGCCAGGTGGAAGTGGATTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAAt                                                       >  1:3293774/1‑150 (MQ=255)
                                                                                                 ggCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGa                                                     <  1:2422749/150‑1 (MQ=255)
                                                                                                  gCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGaa                                                    <  1:3585185/150‑1 (MQ=255)
                                                                                                  gCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGaa                                                    >  2:1093132/1‑150 (MQ=255)
                                                                                                    cAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGGAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAg                                                  >  1:5126014/1‑150 (MQ=255)
                                                                                                     aGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAgc                                                 >  1:769902/1‑150 (MQ=255)
                                                                                                       gTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAgcgc                                               >  1:1942249/1‑150 (MQ=255)
                                                                                                         ggAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCaa                                             >  1:1216604/1‑150 (MQ=255)
                                                                                                         ggAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCaa                                             >  2:1481047/1‑150 (MQ=255)
                                                                                                            aGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAggg                                          <  1:3944847/150‑1 (MQ=255)
                                                                                                               ggATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGc                                       <  1:5324322/150‑1 (MQ=255)
                                                                                                                 aTTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCaa                                     >  1:507398/1‑150 (MQ=255)
                                                                                                                 aTTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCaa                                     >  1:3390993/1‑150 (MQ=255)
                                                                                                                  tttCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAAc                                    >  1:1327588/1‑150 (MQ=255)
                                                                                                                  tttCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACAACGCAGACGAATGAAAGCGCAAGGGGGCAAc                                    >  2:6673288/1‑150 (MQ=255)
                                                                                                                     cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                                 <  2:6256130/150‑1 (MQ=255)
                                                                                                                     cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                                 >  1:6749685/1‑150 (MQ=255)
                                                                                                                              gCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCcac                        <  2:3587115/150‑1 (MQ=255)
                                                                                                                               cGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCcaca                       >  2:192697/1‑150 (MQ=255)
                                                                                                                                    gtTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTtc                  <  1:4126977/150‑1 (MQ=255)
                                                                                                                                     ttGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGGGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTtct                 >  2:4315513/1‑150 (MQ=255)
                                                                                                                                       gCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCa               >  2:5665537/1‑150 (MQ=255)
                                                                                                                                       gCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCa               <  2:3506326/150‑1 (MQ=255)
                                                                                                                                        cTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGaaa                                                   >  1:716514/1‑113 (MQ=255)
                                                                                                                                        cTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAg              >  2:5661498/1‑150 (MQ=255)
                                                                                                                                           atgCAGGGCGTCGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTCGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCg           >  1:4698907/3‑150 (MQ=255)
                                                                                                                                            tgCAGGGCGTCGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGt          >  1:4133475/2‑150 (MQ=255)
                                                                                                                                             gCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTg         >  1:4542043/1‑150 (MQ=255)
                                                                                                                                                 ggCGACGCCATGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTATGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCACGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAAc     >  2:4744289/1‑150 (MQ=255)
                                                                                                                                                   cGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAACGc   >  2:2674412/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAACGCC  >  NC_002947/3541651‑3541948

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: