Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,982,561 T→A 5.8% intergenic (‑33/+38) PP_4387 ← / ← flgE hypothetical protein/flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,982,5610TA 5.8% 626.6 / 17.7 188intergenic (‑33/+38)PP_4387/flgEhypothetical protein/flagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base T (72/105);  new base A (5/6);  total (77/111)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

CACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAG  >  NC_002947/4982413‑4982709
                                                                                                                                                    |                                                                                                                                                    
cacaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCACTTTg                                                                                                                                                     <  2:1639266/150‑1 (MQ=255)
 acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                    <  1:1686222/150‑1 (MQ=255)
 acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                    <  2:4541167/150‑1 (MQ=255)
  caGGTCGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCg                                                                                                                                                   <  1:545932/150‑1 (MQ=255)
  caGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCg                                                                                                                                                   <  1:1084336/150‑1 (MQ=255)
   aGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCgg                                                                                                                                                  <  2:3439397/150‑1 (MQ=255)
    ggTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGc                                                                                                                                                 <  2:4097218/150‑1 (MQ=255)
      tAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGccc                                                                                                                                               <  1:2033645/150‑1 (MQ=255)
      gaGCGGCCGGTGACCAGGCTTTGCCAGTGATCGCACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGccc                                                                                                                                               <  2:5401044/149‑1 (MQ=255)
       aGCGGCCGGTGACCAGGCTTTGCCATTGATAGACCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGcccc                                                                                                                                              <  2:4633416/150‑1 (MQ=255)
       aGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGcccc                                                                                                                                              <  1:736247/150‑1 (MQ=255)
         cGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCaa                                                                                                                                            <  1:212855/150‑1 (MQ=255)
          ggCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAt                                                                                                                                           <  1:2688987/150‑1 (MQ=255)
           gCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaaga                                                                                                                                          <  1:3085541/150‑3 (MQ=255)
             cGGTGACCAGGCTTTGCCATTGATAGACCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTATTGGGGCCGCAAAGCGGCCCCaagaat                                                                                                                                        <  1:3187877/150‑5 (MQ=255)
             cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAc                                                                                                                                        <  1:1079192/150‑1 (MQ=255)
             cGGTCACCAGGATTTGCAATTGATACAACGGGGCGGGGGCAGCGGGGGCAGCCACAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTACACATTATTGGGGCCGCAAAGCGGCCCCaagaat                                                                                                                                        <  1:538425/150‑5 (MQ=255)
              ggTGACCAGGCTTTGCCATTGATAGACCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACa                                                                                                                                       <  2:2643465/150‑1 (MQ=255)
                tGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACaaa                                                                                                                                     <  2:3576679/150‑1 (MQ=255)
                tGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACaaa                                                                                                                                     <  2:5344759/150‑1 (MQ=255)
                  aCCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAaca                                                                                                                                   >  2:3749941/1‑150 (MQ=255)
                   ccAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAacac                                                                                                                                  <  1:4477114/150‑1 (MQ=255)
                   ccAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAacac                                                                                                                                  <  2:1047415/150‑1 (MQ=255)
                         tttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                            >  2:2190011/1‑150 (MQ=255)
                         tttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                            <  2:1614741/150‑1 (MQ=255)
                          ttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGcc                                                                                                                           <  2:765473/150‑1 (MQ=255)
                            gCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                                         <  1:4563864/150‑1 (MQ=255)
                            gCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                                         <  2:2292838/150‑1 (MQ=255)
                            gCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                                         >  2:2990645/1‑150 (MQ=255)
                               aTTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATAc                                                                                                                      >  1:3371151/1‑150 (MQ=255)
                               aTTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATAc                                                                                                                      >  1:5130623/1‑150 (MQ=255)
                                 tgttAGACCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGGAAGGGGATCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACaa                                                                                                                    <  2:2433959/147‑1 (MQ=255)
                                 tGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACaa                                                                                                                    <  2:1069293/150‑1 (MQ=255)
                                 tGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACaa                                                                                                                    <  1:1963966/150‑1 (MQ=255)
                                   aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatc                                                                                                                  >  1:5071410/1‑150 (MQ=255)
                                     aGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatcat                                                                                                                >  2:4627652/1‑150 (MQ=255)
                                      gAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAaacatc                                                                                                               >  2:1215513/1‑150 (MQ=255)
                                       aaCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAtcatca                                                                                                              >  1:1324682/1‑150 (MQ=255)
                                        aCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAg                                                                                                             <  2:4161635/150‑1 (MQ=255)
                                        aCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAg                                                                                                             <  1:244204/150‑1 (MQ=255)
                                        aCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAg                                                                                                             <  2:3445265/150‑1 (MQ=255)
                                        aCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAg                                                                                                             <  2:1501137/150‑1 (MQ=255)
                                          ggggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGt                                                                                                           >  2:941135/1‑150 (MQ=255)
                                          ggggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaaga                                                                                                                                          >  2:3184786/1‑117 (MQ=255)
                                             gcgggggcagcgggggcaGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCAt                                                                                                        <  1:429437/150‑1 (MQ=255)
                                             gcgggggcagcgggggcaGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCAt                                                                                                        <  1:3158162/150‑1 (MQ=255)
                                              cgggggctgcgggggcagGCATAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATc                                                                                                       <  2:1407095/150‑1 (MQ=255)
                                               gggggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATa                                                                                                                       >  2:1215959/1‑133 (MQ=255)
                                                  ggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa                                                                                                   <  1:4241381/150‑1 (MQ=255)
                                                   gcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAat                                                                                                  >  1:3913590/1‑150 (MQ=255)
                                                     agcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatga                                                                                                >  2:4008829/1‑150 (MQ=255)
                                                     agcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatga                                                                                                >  2:3500189/1‑150 (MQ=255)
                                                      gcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat                                                                                               >  1:4682845/1‑150 (MQ=255)
                                                           ggcaggcaGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCt                                                                                          <  2:4377204/150‑1 (MQ=255)
                                                            gcaggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTg                                                                                         >  2:1811147/1‑150 (MQ=255)
                                                             caggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGc                                                                                        <  1:4900636/150‑1 (MQ=255)
                                                                gcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAATACAAACACCAAAGCCAATACAATCAACAGGTCATCTGAATGATGGTCTGCatg                                                                                     >  1:4709709/1‑150 (MQ=255)
                                                                  agAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTGTGCGGCCCCCAGTACAAACCCCAAAGCCAATACCCTCATCAGGTCATCTGAATGAGGGGCTGCatgag                                                                                   >  2:4672850/1‑149 (MQ=255)
                                                                   gAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                                  <  1:3627651/150‑1 (MQ=255)
                                                                    aaCCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACa                                                                                                                                       >  2:2855734/1‑96 (MQ=255)
                                                                    aaCCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACa                                                                                                                                       <  1:2855734/96‑1 (MQ=255)
                                                                      ccagcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTTCATGATGGTg                                                                               >  1:3818412/1‑150 (MQ=255)
                                                                       cagcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGc                                                                              >  1:3706423/1‑150 (MQ=255)
                                                                            ccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCa                                                                         <  2:3335491/150‑1 (MQ=255)
                                                                            ccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGATTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCa                                                                         >  2:137423/1‑150 (MQ=255)
                                                                                  aGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGaa                                                                   <  1:1170392/150‑1 (MQ=255)
                                                                                      gAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                                  <  2:3843582/131‑1 (MQ=255)
                                                                                      gAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                                  >  1:3843582/1‑131 (MQ=255)
                                                                                        aCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTc                                                             <  1:2893016/150‑1 (MQ=255)
                                                                                          aaGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCtt                                                           <  1:3992876/150‑1 (MQ=255)
                                                                                            ggCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTgg                                                         <  2:4770077/150‑1 (MQ=255)
                                                                                            ggCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTAgg                                                         >  2:2220102/1‑150 (MQ=255)
                                                                                             gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                                        >  2:2747094/1‑150 (MQ=255)
                                                                                             gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                                        >  2:19331/1‑150 (MQ=255)
                                                                                               cAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCGGTGCGGCCCCAATTACAAACACCAAAGCCAATACAAACATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAAAGGTCTTGGCAt                                                      >  1:2527265/1‑150 (MQ=255)
                                                                                                aGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAtt                                                     >  2:5214146/1‑150 (MQ=255)
                                                                                                    aCGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGcc                                                 >  2:2306704/1‑150 (MQ=255)
                                                                                                     cgtcAGCCGCTTCATTCGCTTGAACCTTCAAATTCTTGGTGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGTATGATGGTCTGAATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCt                                                <  1:2939832/147‑1 (MQ=255)
                                                                                                      gACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTg                                               <  2:4517645/150‑1 (MQ=255)
                                                                                                      gACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTg                                               <  1:3564612/150‑1 (MQ=255)
                                                                                                       aCAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTgg                                              <  2:286757/150‑1 (MQ=255)
                                                                                                         aGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTa                                            <  1:247366/150‑1 (MQ=255)
                                                                                                              cTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGc                                       >  2:4294198/1‑150 (MQ=255)
                                                                                                              cTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAATACAAACACCAAAGCCAAAACAAACAACAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAATGGCCTGGTAGtttc                                       >  1:14998/1‑148 (MQ=255)
                                                                                                                       cgTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCt                              <  1:2243265/150‑1 (MQ=255)
                                                                                                                          gAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGGTGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGa                           >  2:2531958/1‑150 (MQ=255)
                                                                                                                          gAACCTTCACATTCTTGGGGCCGCGTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCGGGTAGTTGCTCTGCGCCTTGa                           >  2:513890/1‑150 (MQ=255)
                                                                                                                           aaCCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGAt                          <  1:4953290/150‑1 (MQ=255)
                                                                                                                              cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGa                                                                                                    <  1:1564969/73‑1 (MQ=255)
                                                                                                                              cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                       >  1:887814/1‑150 (MQ=255)
                                                                                                                              cTTCACATTCTTGGGGCCGCGTTGCGGCCCCAATTACAAAAAACAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATTGTCATGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                       >  1:1797985/1‑150 (MQ=255)
                                                                                                                                      tCTTGGGGCCGCTGTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCa               <  1:2587669/150‑1 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTc           <  2:1501343/150‑1 (MQ=255)
                                                                                                                                          ggggCCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTc           >  2:3902907/1‑150 (MQ=255)
                                                                                                                                             gCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGTTCAGGTTTACCAGCTCGcc        >  2:2162933/1‑150 (MQ=255)
                                                                                                                                              ccGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACACTCAGCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCg       <  1:2193211/150‑1 (MQ=255)
                                                                                                                                                gCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGt     >  1:1776307/1‑150 (MQ=255)
                                                                                                                                                 cTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTc    <  1:4627652/150‑1 (MQ=255)
                                                                                                                                                  tttGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCCGCTCGCCGGTCa   >  2:3240702/1‑150 (MQ=255)
                                                                                                                                                   ttGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAg  >  2:2333985/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAG  >  NC_002947/4982413‑4982709

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