Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,810,512 A→G 6.9% noncoding (413/434 nt) PP_mr35 ← cspA regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,810,5120AG 6.9% 692.8 / 27.3 230noncoding (413/434 nt)PP_mr35cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (100/114);  new base G (8/8);  total (108/122)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

Only 100 of 234 total aligned reads displayed.

CGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810363‑2810658
                                                                                                                                                      |                                                                                                                                                  
cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACTCGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                    <  2:4684632/150‑1 (MQ=255)
cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                    >  2:1714171/1‑150 (MQ=255)
  cAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAgg                                                                                                                                                  <  2:703140/150‑1 (MQ=255)
   aaaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAggg                                                                                                                                                 >  2:5798704/1‑150 (MQ=255)
    aaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGc                                                                                                                                                <  2:8103305/150‑1 (MQ=255)
       tgtgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCggg                                                                                                                                             <  2:6008137/150‑1 (MQ=255)
         tgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAACGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCACTCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCCAAGGGCGGGtt                                                                                                                                           >  2:6487128/1‑150 (MQ=255)
           aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGtttt                                                                                                                                         >  2:2366936/1‑150 (MQ=39)
           aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAAACCGCCACAAGGGCGGGtttt                                                                                                                                         >  1:5632943/1‑150 (MQ=255)
           aaaaGACTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttt                                                                                                                                         >  2:6290619/1‑150 (MQ=255)
            aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGttttt                                                                                                                                        <  2:3915594/150‑1 (MQ=39)
            aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  1:4612709/150‑1 (MQ=255)
            aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:1960129/150‑1 (MQ=255)
            aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:693400/150‑1 (MQ=255)
            aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:7662745/150‑1 (MQ=255)
            aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:5553383/150‑1 (MQ=255)
            aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  1:4002522/150‑1 (MQ=255)
             aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGtttttt                                                                                                                                       <  1:932701/150‑1 (MQ=39)
             aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                       <  1:2886013/150‑1 (MQ=255)
             aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                       <  2:3103989/150‑1 (MQ=255)
               gCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTa                                                                                                                                     >  2:9122298/1‑150 (MQ=255)
                 cTGCGCGGTGGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACa                                                                                                                                   <  1:1853876/150‑1 (MQ=255)
                 cTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACa                                                                                                                                   <  1:5761263/150‑1 (MQ=255)
                  tGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAg                                                                                                                                  >  1:5971853/1‑150 (MQ=255)
                  tGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAg                                                                                                                                  >  2:1470696/1‑150 (MQ=255)
                       ggTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCACTCCCGCCCCTTCAGATGTTGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATc                                                                                                                             >  1:2641552/1‑150 (MQ=255)
                           gCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGaa                                                                                                                         >  2:1073828/1‑150 (MQ=255)
                            ccTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAt                                                                                                                        <  2:4969573/150‑1 (MQ=255)
                                gTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACt                                                                                                                    <  1:5993303/150‑1 (MQ=255)
                                        gggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGa                                                                                                            <  1:5431417/150‑1 (MQ=255)
                                        gggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGa                                                                                                            <  1:6489807/150‑1 (MQ=255)
                                         ggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGaa                                                                                                           <  2:331373/150‑1 (MQ=255)
                                            tCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCt                                                                                                        <  2:3835119/150‑1 (MQ=255)
                                              aCCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCtgt                                                                                                      <  2:6559564/150‑1 (MQ=255)
                                                ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATTCATCACCCAAAAAAAAAACCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGc                                                                                                    >  2:2018895/1‑150 (MQ=255)
                                                ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACACGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGc                                                                                                    >  2:5603620/1‑150 (MQ=255)
                                                    tggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTTCAGCTATCCGAATTACTCCGGACGAAc                                                                                                          >  1:6604462/4‑140 (MQ=255)
                                                    tggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAAc                                                                                                          <  2:6604462/137‑1 (MQ=255)
                                                    tggcacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGcc                                                                                                >  2:8251125/4‑150 (MQ=255)
                                                    gAACACGGGCAAAGCCCGTGCCAGTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCCTAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGcc                                                                                                <  2:1013590/150‑1 (MQ=255)
                                                      acacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAATCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTg                                                                                              >  1:3744922/1‑150 (MQ=255)
                                                        acGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCCGATGATGCATCACCCATAAAAAAACCCGCCCCAAGGGCGGGTTTTTTTACAGCTCTCCGAATTACTCCGGACGAACCTGTGCAGCCTGCa                                                                                            >  1:7866200/1‑150 (MQ=255)
                                                         cGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCAt                                                                                           <  2:6090428/150‑1 (MQ=255)
                                                          gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCgg                                                                                                                                              >  1:635044/1‑98 (MQ=255)
                                                          gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCgg                                                                                                                                              <  2:635044/98‑1 (MQ=255)
                                                          gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATg                                                                                          <  2:2259121/150‑1 (MQ=255)
                                                            gCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGcc                                                                                        >  2:7588216/1‑150 (MQ=255)
                                                             cAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCCTCACCCATAAAAAAACCCGCCCCAAGGGCGGGTTTTTTTACAGCTATCCGACTTACTCCGGACGAACCTGTGCAGCCTGCATGccc                                                                                       >  2:5432012/1‑150 (MQ=255)
                                                              aaaGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGcccg                                                                                      >  1:3267264/1‑149 (MQ=255)
                                                                              aCCCGACAGTCGGAGTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCAAAACAAAACCCGCCCCCAGGGCGGGTTTTTTTACAGCAATCCGAATTACTCCGGACGACCCCGTGCAGCCTGCCTGCCCTTCTGGCCGCGCTCGGc                                                                      <  1:3019538/150‑1 (MQ=255)
                                                                                  tACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcga                                                                    <  1:3891714/148‑1 (MQ=255)
                                                                                  tACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCAACCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCAGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgacg                                                                  >  2:4524934/1‑150 (MQ=255)
                                                                                    cAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGaa                                                                <  1:8918452/150‑1 (MQ=255)
                                                                                     aGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCAGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAg                                                               >  2:3021804/1‑150 (MQ=255)
                                                                                       tCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGa                                                             >  1:3280074/1‑150 (MQ=255)
                                                                                         ggTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAACAAACCCGCCACAAGGGCGGGTTTTTTTACAGCGATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaaa                                                           <  1:8334970/150‑1 (MQ=255)
                                                                                          gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaaa                                                           <  1:2788183/150‑1 (MQ=255)
                                                                                               gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACCGTCt                                                     >  1:6485432/1‑150 (MQ=255)
                                                                                               gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  2:3469128/1‑150 (MQ=255)
                                                                                               gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  1:1038673/1‑150 (MQ=255)
                                                                                               gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  2:6758936/1‑150 (MQ=255)
                                                                                               gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCGGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  1:4475468/1‑150 (MQ=255)
                                                                                                cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCAGAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    <  2:1391307/150‑1 (MQ=255)
                                                                                                cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    <  2:2944586/150‑1 (MQ=255)
                                                                                                cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    >  1:1927956/1‑150 (MQ=255)
                                                                                                 aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTgg                                                   <  2:4758361/150‑1 (MQ=255)
                                                                                                  ggCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGc                                                  <  1:8381477/150‑1 (MQ=255)
                                                                                                    cAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACCGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCct                                                <  2:1264106/150‑1 (MQ=255)
                                                                                                     aaTCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCctt                                               <  1:5164025/150‑1 (MQ=255)
                                                                                                     aaTCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAGTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGACACAGTCTGGCctt                                               >  2:7293942/1‑150 (MQ=255)
                                                                                                        cccGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCcttctt                                            >  1:2101160/1‑150 (MQ=255)
                                                                                                           gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         >  1:2244062/1‑150 (MQ=255)
                                                                                                           gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         <  1:6119738/150‑1 (MQ=255)
                                                                                                           gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         <  1:2047640/150‑1 (MQ=255)
                                                                                                           gCCCCGTCAGATGATGCATCACCCACACCACACCCCGCCCCAAGGGCGGGTTTTTTTACAGCTCTCCGAATTACGCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         <  1:8494332/150‑1 (MQ=255)
                                                                                                            ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAg                                        <  1:1711208/150‑1 (MQ=255)
                                                                                                            ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTATTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAg                                        <  1:2509837/150‑1 (MQ=255)
                                                                                                            ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTCTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAg                                        <  1:8512920/150‑1 (MQ=255)
                                                                                                                ttCAGATGATGCATCACCCATAGAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGctc                                    <  2:4606345/150‑1 (MQ=255)
                                                                                                                     atgatgCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGaa                               <  1:5498176/150‑1 (MQ=255)
                                                                                                                            aTCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCa                        <  2:1498675/150‑1 (MQ=255)
                                                                                                                                 ccATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTc                   <  1:3509216/150‑1 (MQ=255)
                                                                                                                                   aTAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGa                 >  1:5109579/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATgg              >  1:4179089/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGcc                                                                               >  1:6323625/1‑84 (MQ=37)
                                                                                                                                        aaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  1:1405536/1‑150 (MQ=39)
                                                                                                                                        aaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTCTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  2:2790461/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  1:5534055/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCcttc                                              >  1:2559640/1‑116 (MQ=255)
                                                                                                                                         aaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCcttc                                              <  2:2559640/116‑1 (MQ=255)
                                                                                                                                         aaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCt            <  1:2867782/150‑1 (MQ=255)
                                                                                                                                         aaaaCCCGCCACAAGGGCGGGTTGTTTTACAGCTATACGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGTTGGCt            >  2:2210773/1‑150 (MQ=255)
                                                                                                                                          aatcccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           >  2:3779371/4‑150 (MQ=255)
                                                                                                                                          aaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           <  2:6865836/150‑1 (MQ=255)
                                                                                                                                          aaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           <  2:8123590/150‑1 (MQ=255)
                                                                                                                                          aaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           >  1:4103516/1‑150 (MQ=255)
                                                                                                                                              ccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa       <  1:8000155/150‑1 (MQ=255)
                                                                                                                                                   acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCTGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  >  2:8737281/1‑150 (MQ=255)
                                                                                                                                                   acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  <  1:3656921/150‑1 (MQ=255)
                                                                                                                                                   acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  <  1:1470696/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810363‑2810658

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: