Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,223,266 G→T 12.9% intergenic (‑530/+30) PP_2823 ← / ← PP_2824 methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,223,2660GT12.9% 255.6 / 21.7 93intergenic (‑530/+30)PP_2823/PP_2824methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (35/46);  new base T (7/5);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

AGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATCA  >  NC_002947/3223117‑3223413
                                                                                                                                                     |                                                                                                                                                   
aGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAg                                                                                                                                                     <  2:2124418/150‑1 (MQ=255)
 gACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAgg                                                                                                                                                    <  1:1180720/150‑1 (MQ=255)
    gCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGaa                                                                                                                                                 <  1:96163/150‑1 (MQ=255)
      aaTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGGGGGGGTGACGCCTTCCAGGGAAgg                                                                                                                                               >  1:1178532/1‑150 (MQ=255)
       aTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGgc                                                                                                                                              >  2:1691189/1‑150 (MQ=255)
            gcCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCgg                                                                                                                                         >  1:445533/1‑150 (MQ=255)
            gcCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCgg                                                                                                                                         >  2:1579528/1‑150 (MQ=255)
                     tGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCGGGACGGCGCCAGGTCTACCgg                                                                                                                                >  2:438317/1‑150 (MQ=255)
                     tGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                >  1:2383695/1‑150 (MQ=255)
                       cAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGATGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGa                                                                                                                              >  1:240521/1‑150 (MQ=255)
                        aaaaaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGGCGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAg                                                                                                                             <  2:415989/150‑1 (MQ=255)
                         aaaaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAgg                                                                                                                            >  2:391036/1‑150 (MQ=255)
                          aaaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGc                                                                                                                           <  2:2383695/150‑1 (MQ=255)
                              cagggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACt                                                                                                                       <  2:445533/150‑1 (MQ=255)
                                gggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTa                                                                                                                     <  1:445717/150‑1 (MQ=255)
                                  gcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAgg                                                                                                                   <  1:233872/150‑1 (MQ=255)
                                  gcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAgg                                                                                                                   >  1:76231/1‑150 (MQ=255)
                                   caggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGa                                                                                                                  <  2:910772/150‑1 (MQ=255)
                                   caggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGa                                                                                                                  <  2:2054925/150‑1 (MQ=255)
                                   caggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGa                                                                                                                  <  1:181855/150‑1 (MQ=255)
                                       cgcTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTg                                                                                                              <  1:1579528/150‑1 (MQ=255)
                                         ccaTCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGc                                                                                                            <  2:740558/148‑1 (MQ=255)
                                         cTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGc                                                                                                            >  1:354004/1‑150 (MQ=255)
                                            tCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGa                                                                                                                              <  2:275562/129‑1 (MQ=255)
                                            tCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGa                                                                                                                              >  1:275562/1‑129 (MQ=255)
                                            tCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGt                                                                                                         >  1:141248/1‑150 (MQ=255)
                                               gggCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGtttt                                                                                                      <  2:1910375/150‑1 (MQ=255)
                                               gggCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGtttt                                                                                                      <  2:1499013/150‑1 (MQ=255)
                                               gggCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGtttt                                                                                                      <  1:1691189/150‑1 (MQ=255)
                                                 gCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                                    <  1:248678/150‑1 (MQ=255)
                                                 gCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                                    <  2:334758/150‑1 (MQ=255)
                                                 gCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                                    <  2:141248/150‑1 (MQ=255)
                                                 gCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCAGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                                    <  2:1253120/150‑1 (MQ=255)
                                                  ccTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCggg                                                                                                                               >  1:34060/1‑122 (MQ=255)
                                                  ccTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCggg                                                                                                                               <  2:34060/122‑1 (MQ=255)
                                                   cTGGTGCTTGGGTTGGATCGGCAAGGCGCGGCGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGCTTGGCGGGGGTGACGCCTTCCAGGGAAGACCCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGCTTCTGGGt                                                                                                  >  2:1162793/1‑150 (MQ=255)
                                                     ggTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTgg                                                                                                <  2:1434199/150‑1 (MQ=255)
                                                       tGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTAcccaagca                                                                                                                            >  1:987556/1‑114 (MQ=255)
                                                       tGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTAcccaagca                                                                                                                            <  2:987556/120‑7 (MQ=255)
                                                          ttGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACtt                                                                                                               <  2:2008596/130‑1 (MQ=255)
                                                          ttGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACtt                                                                                                               >  1:2008596/1‑130 (MQ=255)
                                                           tGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAAt                                                                                          <  2:2065152/150‑1 (MQ=255)
                                                              gTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                       >  1:2197368/1‑150 (MQ=255)
                                                                  gATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGAATTGTAGCGCAGGAGAGCTAATTAGCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACa                                                                                   >  1:1488946/1‑150 (MQ=255)
                                                                       gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCTGGAAGgcgtc                                                                                                                                           <  1:1682181/89‑3 (MQ=255)
                                                                       gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCTGGAAGgcgtc                                                                                                                                           >  2:1682181/1‑87 (MQ=255)
                                                                       gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                >  2:1285947/1‑100 (MQ=255)
                                                                       gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                <  1:1285947/100‑1 (MQ=255)
                                                                       gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGaa                                                                              <  1:973567/150‑1 (MQ=255)
                                                                        cAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAg                                                                             <  1:820287/150‑1 (MQ=255)
                                                                                   gcggcAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGtcaatcaa                                                                  >  2:821594/1‑150 (MQ=255)
                                                                                    cggcAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAg                                                                 <  2:2197368/150‑1 (MQ=255)
                                                                                                         aGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGtcaatc                                                                    >  1:2096873/1‑126 (MQ=255)
                                                                                                         aGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGtcaatc                                                                    <  2:2096873/126‑1 (MQ=255)
                                                                                                           gcggcgAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAAc                                          <  2:354004/150‑1 (MQ=255)
                                                                                                             ggcgAAGGCACGAGCGGTTGGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACgg                                        <  2:1178532/150‑1 (MQ=255)
                                                                                                                gAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                <  1:890120/59‑1 (MQ=255)
                                                                                                                gAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                >  2:890120/1‑59 (MQ=255)
                                                                                                                 aaGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTc                                    <  1:438317/150‑1 (MQ=255)
                                                                                                                  aGACACGAGCGGTGTGGTGGGGGTGATGCCCTCCAGGGAAGGCGGCGGGTGTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCg                                   <  2:240521/150‑1 (MQ=255)
                                                                                                                   gACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCgg                                  <  2:76231/150‑1 (MQ=255)
                                                                                                                     cacGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGAc                                <  1:946643/150‑1 (MQ=255)
                                                                                                                        gAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGCCTTGGCGGCTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCg                             >  2:272685/1‑150 (MQ=255)
                                                                                                                        gAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                       <  1:1932004/92‑1 (MQ=255)
                                                                                                                        gAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                       >  2:1932004/1‑92 (MQ=255)
                                                                                                                            ggTTTGGTGGGGGTGACGGCTTCCAGGGAAGGCGGCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGTAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCTCTCTGCGCTAACGGGTTCGGACTCGAAAg                         <  1:1162793/150‑1 (MQ=255)
                                                                                                                              tttGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCGAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGgc                       >  2:1739543/1‑150 (MQ=255)
                                                                                                                                  gtggGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGgcagca                   >  1:1924710/1‑150 (MQ=255)
                                                                                                                                     ggggTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGtttt                                                                                                      <  2:772274/64‑1 (MQ=255)
                                                                                                                                     ggggTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGtttt                                                                                                      >  1:772274/1‑64 (MQ=255)
                                                                                                                                      gggTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat               <  2:1488946/150‑1 (MQ=255)
                                                                                                                                        ccggcGCCTTCCCTGGAAGGCGTCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTc                                    <  1:1544341/123‑1 (MQ=255)
                                                                                                                                        ccggcGCCTTCCCTGGAAGGCGTCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTc                                    >  2:1544341/5‑127 (MQ=255)
                                                                                                                                        ccggcGCCTTCCCTGGAAGGCGTCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGa             >  2:1403281/5‑150 (MQ=255)
                                                                                                                                         tGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGa             <  1:1380600/149‑1 (MQ=255)
                                                                                                                                         tGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGa             >  2:1380600/1‑149 (MQ=255)
                                                                                                                                         cggcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGaa                                >  2:1343133/4‑129 (MQ=255)
                                                                                                                                         cggcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGAc                                <  1:1343133/127‑1 (MQ=255)
                                                                                                                                          ggcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagc           >  2:562586/3‑150 (MQ=255)
                                                                                                                                           gcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGTCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGt                                      <  2:1772310/121‑1 (MQ=255)
                                                                                                                                           gcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGTCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGt                                      >  1:1772310/2‑122 (MQ=255)
                                                                                                                                            cGCCTTCCAGGGAAGGCGCTGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagcca         <  2:2243476/150‑1 (MQ=255)
                                                                                                                                            cGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagcca         <  1:821594/150‑1 (MQ=255)
                                                                                                                                             gCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGt                                                                         >  1:1764079/1‑85 (MQ=255)
                                                                                                                                             gCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGt                                                                         <  2:1764079/85‑1 (MQ=255)
                                                                                                                                             gCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat               >  2:897755/1‑143 (MQ=255)
                                                                                                                                             gCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat               <  1:897755/143‑1 (MQ=255)
                                                                                                                                             gCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagccag        >  2:814579/1‑150 (MQ=255)
                                                                                                                                             gCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTCAGCAGACGGGTTCGGACTCGAAAGGCAGCATCATGagccag        >  1:1231907/1‑150 (MQ=255)
                                                                                                                                                ttCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAgccagcca     >  1:1359397/1‑150 (MQ=255)
                                                                                                                                                  ccAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATc   <  1:814579/150‑1 (MQ=255)
                                                                                                                                                   cAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATCa  >  1:1015356/1‑150 (MQ=255)
                                                                                                                                                   cAGGGAAGGCGCCGGGTCTACCGGAAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTaa                                           >  2:217464/1‑109 (MQ=255)
                                                                                                                                                   cAGGGAAGGCGCCGGGTCTACCGGAAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTaa                                           <  1:217464/109‑1 (MQ=255)
                                                                                                                                                   cAGGGAAGGCGCCGGGGCTACCGGGAGGCGACTTAGGACTTGGCGGCTTTGGGTGGAGAAGTCATGACAATGAAGGCGTCCATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATcc  >  2:97056/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATCA  >  NC_002947/3223117‑3223413

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