Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,576,668 G→C 6.6% P235P (CCG→CCC) ‡ yneJ → LysR family DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,576,6680GC 6.6% 1155.7 / 50.1 344P235P (CCG→CCC) ‡yneJLysR family DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (156/165);  new base C (11/12);  total (167/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 353 total aligned reads displayed.

TACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAAG  >  NC_002947/3576519‑3576816
                                                                                                                                                     |                                                                                                                                                    
 aCGCCTTCCGCGCCCACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgc                                                                                                                                                     >  1:8403198/1‑150 (MQ=255)
 aCGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgc                                                                                                                                                     >  1:800284/1‑150 (MQ=255)
 aCGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgc                                                                                                                                                     >  1:3238334/1‑150 (MQ=255)
      ttCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGATGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCgcagc                                                                                                                                                <  1:2397616/150‑1 (MQ=255)
        ccGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCAt                                                                                                                                              <  1:2053161/150‑1 (MQ=255)
         cgcgCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCATg                                                                                                                                             <  2:2319082/150‑1 (MQ=255)
         cgcgCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATg                                                                                                                                             >  1:7895067/1‑150 (MQ=255)
           cgcCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCt                                                                                                                                           <  2:2709285/150‑1 (MQ=255)
            gcCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTTCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTc                                                                                                                                          <  1:6123635/150‑1 (MQ=255)
            gcCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTc                                                                                                                                          >  2:2940007/1‑150 (MQ=255)
            gcCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTc                                                                                                                                          <  2:341761/150‑1 (MQ=255)
              cAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCATggccag                                                                                                                                        <  1:1911378/150‑6 (MQ=255)
              cAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGa                                                                                                                                        <  1:965596/150‑1 (MQ=255)
              cAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGa                                                                                                                                        >  2:4934290/1‑150 (MQ=255)
              cAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGa                                                                                                                                        >  2:748776/1‑150 (MQ=255)
                  tGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAat                                                                                                                                    <  2:2880040/150‑1 (MQ=255)
                  tGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAat                                                                                                                                    <  2:4150286/150‑1 (MQ=255)
                  tGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAat                                                                                                                                    >  2:911039/1‑150 (MQ=255)
                   gTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAata                                                                                                                                   <  1:7554354/150‑1 (MQ=255)
                     tCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATg                                                                                                                                 <  1:7661494/150‑1 (MQ=255)
                      cATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGc                                                                                                                                <  1:812798/150‑1 (MQ=255)
                       aTACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGcc                                                                                                                               >  2:2190863/1‑150 (MQ=255)
                         aCCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCg                                                                                                                             >  2:7015309/1‑150 (MQ=255)
                         aCCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCg                                                                                                                             >  2:3290314/1‑150 (MQ=255)
                          ccGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCgg                                                                                                                            <  2:5079776/150‑1 (MQ=255)
                          ccGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCgg                                                                                                                            <  1:1182562/150‑1 (MQ=255)
                            gCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGt                                                                                                                          <  1:1452961/150‑1 (MQ=255)
                            gCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGt                                                                                                                          <  1:6249278/150‑1 (MQ=255)
                             ccATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTg                                                                                                                         <  1:1696488/150‑1 (MQ=255)
                                 cACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTgcagc                                                                                                                     >  1:3757784/1‑150 (MQ=255)
                                 cACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTgcagc                                                                                                                     >  1:52293/1‑150 (MQ=255)
                                    ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAcc                                                                                                                  >  2:6234847/1‑150 (MQ=255)
                                    ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAcc                                                                                                                  >  1:8786504/1‑150 (MQ=255)
                                    ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAcc                                                                                                                  >  2:1444322/1‑150 (MQ=255)
                                    ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAcc                                                                                                                  >  1:7706663/1‑150 (MQ=255)
                                    ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAcc                                                                                                                  >  1:7183034/1‑150 (MQ=255)
                                    ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAcc                                                                                                                  >  1:1664665/1‑150 (MQ=255)
                                      cGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGt                                                                                                                >  2:5982865/1‑150 (MQ=255)
                                        aaaaCTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCa                                                                                                              <  2:5751687/150‑1 (MQ=255)
                                           aCTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAgcg                                                                                                           >  1:4049634/1‑150 (MQ=255)
                                               gTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTg                                                                                                       >  2:2615292/1‑150 (MQ=255)
                                               gTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTg                                                                                                       >  1:2425042/1‑150 (MQ=255)
                                                 tCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGc                                                                                                     <  1:3036814/150‑1 (MQ=255)
                                                       aGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGt                                                                                               >  1:5880217/1‑150 (MQ=255)
                                                         gACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTTCCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                                                             <  1:6298782/150‑1 (MQ=255)
                                                         gACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                                                             >  2:7706964/1‑150 (MQ=255)
                                                         gACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                                                             >  1:7310510/1‑150 (MQ=255)
                                                          aCCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgg                                                                                            >  1:3607790/1‑150 (MQ=255)
                                                           ccAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                           <  1:106284/150‑1 (MQ=255)
                                                           ccAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                           <  1:6602844/150‑1 (MQ=255)
                                                           ccAGGCGGTACCGGGCAAGATGCACGAAATGGAGTCGTACCACGGCATGTGGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCAGGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                                                           <  2:8642509/150‑1 (MQ=255)
                                                             aGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCggagg                                                                                         <  1:3418654/150‑1 (MQ=255)
                                                             aGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCggagg                                                                                         >  2:3620202/1‑150 (MQ=255)
                                                                 ggTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACtt                                                                                     >  2:3057543/1‑150 (MQ=255)
                                                                  ggaCCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTc                                                                                    <  1:5183541/148‑1 (MQ=255)
                                                                   tACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTcc                                                                                   >  2:6726320/1‑150 (MQ=255)
                                                                    aCCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCg                                                                                  <  2:4470571/150‑1 (MQ=255)
                                                                        ggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTAt                                                                              <  2:4506114/150‑1 (MQ=255)
                                                                        ggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTAt                                                                              <  1:4888763/150‑1 (MQ=255)
                                                                          cAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCt                                                                            >  2:6839715/1‑150 (MQ=255)
                                                                                    gAAATGGGGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctgg                                                                  >  2:8232465/1‑150 (MQ=255)
                                                                                    gAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctgg                                                                  >  1:4809776/1‑150 (MQ=255)
                                                                                     aaaTGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctggc                                                                 <  2:6084159/150‑1 (MQ=255)
                                                                                      aaTGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctggct                                                                >  2:6456165/1‑150 (MQ=255)
                                                                                      aaTGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctggct                                                                <  2:8340468/150‑1 (MQ=255)
                                                                                            gTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATg                                                          >  1:3168133/1‑150 (MQ=255)
                                                                                               gTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGcg                                                       >  2:187664/1‑150 (MQ=255)
                                                                                                   cACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGcgacga                                                   <  2:3246887/150‑1 (MQ=255)
                                                                                                   cACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGcgacga                                                   >  2:5605673/1‑150 (MQ=255)
                                                                                                         aTGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTCCg                                             >  2:7424590/1‑150 (MQ=255)
                                                                                                         aTGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACg                                             <  1:6858545/150‑1 (MQ=255)
                                                                                                         aTGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACg                                             <  1:672524/150‑1 (MQ=255)
                                                                                                            ttGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGtc                                          <  2:4175053/150‑1 (MQ=255)
                                                                                                               gCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTcg                                       >  2:381045/1‑150 (MQ=255)
                                                                                                               gCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTcg                                       <  1:7436465/150‑1 (MQ=255)
                                                                                                                 cTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTcgcg                                     <  2:4254686/150‑1 (MQ=255)
                                                                                                                  tGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCgcgc                                    >  2:4915594/1‑150 (MQ=255)
                                                                                                                   gCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCa                                   >  2:7138388/1‑150 (MQ=255)
                                                                                                                       cAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGttt                               <  1:8047046/150‑1 (MQ=255)
                                                                                                                         gTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGa                             >  1:5493924/1‑150 (MQ=255)
                                                                                                                         gTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGa                             <  2:6618218/150‑1 (MQ=255)
                                                                                                                          tGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAg                            >  2:6133466/1‑150 (MQ=255)
                                                                                                                                 gCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTc                     <  2:87183/150‑1 (MQ=255)
                                                                                                                                    ggCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTg                  <  2:4841908/150‑1 (MQ=255)
                                                                                                                                      ccTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGaa                <  2:634226/150‑1 (MQ=255)
                                                                                                                                         gagcATGCTGCGCGGCAGCATGGTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGtt             >  2:8067458/4‑150 (MQ=255)
                                                                                                                                         gagcATGCTGCGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGtt             >  2:1697524/4‑150 (MQ=255)
                                                                                                                                          agcATGCTGCGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTg            >  1:4551657/3‑150 (MQ=255)
                                                                                                                                          agcATGCTGCGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTg            >  1:1401211/3‑150 (MQ=255)
                                                                                                                                           ccATGCTGCCGCGCAGCATGCTCGACAGTATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGc           >  1:3991051/1‑150 (MQ=255)
                                                                                                                                            cATGCTGCGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt          >  2:5415043/1‑150 (MQ=255)
                                                                                                                                            cATGCTGCGCGGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt          >  1:4350032/1‑150 (MQ=255)
                                                                                                                                                ctgcCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaa      <  1:2670007/150‑1 (MQ=255)
                                                                                                                                                ctgcCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaa      <  1:8819880/150‑1 (MQ=255)
                                                                                                                                                ctgcCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaa      <  1:5832216/150‑1 (MQ=255)
                                                                                                                                                  gcCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaaga    >  1:356543/1‑150 (MQ=255)
                                                                                                                                                   ccGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaagaa   <  2:4697701/150‑1 (MQ=255)
                                                                                                                                                   ccGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTgaagaa   <  1:5599332/150‑1 (MQ=255)
                                                                                                                                                   ccGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCATGAAGTTGCTTgaagaa   <  2:131418/150‑1 (MQ=255)
                                                                                                                                                    cgcgCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGaagaag  >  1:2253282/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGAAGAAG  >  NC_002947/3576519‑3576816

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: