Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,813,536 A→C 5.4% noncoding (416/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,813,5360AC 5.4% 1079.1 / 31.0 335noncoding (416/434 nt)PP_mr35cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (164/151);  new base C (10/10);  total (174/161)
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 349 total aligned reads displayed.

TTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2813387‑2813685
                                                                                                                                                     |                                                                                                                                                     
ttGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCca                                                                                                                                                       >  1:1036142/1‑150 (MQ=255)
ttGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCca                                                                                                                                                       >  2:3428416/1‑150 (MQ=255)
      aaaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAAACCGCCACAAggg                                                                                                                                                 >  2:7047491/1‑150 (MQ=255)
       aaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGc                                                                                                                                                <  1:5474021/150‑1 (MQ=255)
       aaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGc                                                                                                                                                <  1:2258371/150‑1 (MQ=255)
        aaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACCAGGGCg                                                                                                                                               >  1:6052917/1‑150 (MQ=255)
         aTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCgg                                                                                                                                              <  1:5675599/150‑1 (MQ=255)
          tgtgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCggg                                                                                                                                             >  2:291004/1‑150 (MQ=255)
             gAATAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttt                                                                                                                                          <  2:6198935/150‑1 (MQ=255)
             gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCTATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttt                                                                                                                                          <  2:2302278/150‑1 (MQ=255)
             gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCTATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttt                                                                                                                                          <  2:8249686/150‑1 (MQ=255)
             gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttt                                                                                                                                          >  1:7733091/1‑150 (MQ=255)
              aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttt                                                                                                                                         <  1:2839506/150‑1 (MQ=255)
              aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttt                                                                                                                                         <  2:243946/150‑1 (MQ=255)
               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:3969636/150‑1 (MQ=255)
               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:1387371/150‑1 (MQ=255)
               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCAGAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  1:1828216/150‑1 (MQ=255)
               aaaGCCTGCGCGGTCGCCTTGCAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGttttt                                                                                                                                        <  2:2804964/150‑1 (MQ=39)
                acgCCTGCGCGGTCGCCTGGGAGGAGCGGGTTCACCCGCGCACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAACAAAACCCGCCCTTGTGGCGGGtttttt                                                                                                                                       <  1:1713000/148‑1 (MQ=255)
                aaGCCTGCGCGGTCGCCTTGTCGGAGCGGGTTCACCCGCGACCACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGTtttttt                                                                                                                                       <  1:176853/150‑1 (MQ=255)
                aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACCAGGGCGGGtttttt                                                                                                                                       >  2:7729135/1‑150 (MQ=255)
                aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                       <  1:7647384/150‑1 (MQ=255)
                aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                       >  2:2582137/1‑150 (MQ=255)
                aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                       <  1:2088324/150‑1 (MQ=255)
                  gCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTa                                                                                                                                     <  2:7485482/150‑1 (MQ=255)
                    cTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACa                                                                                                                                   <  2:5013314/150‑1 (MQ=255)
                    cTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACa                                                                                                                                   <  2:5062916/150‑1 (MQ=255)
                     tGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAg                                                                                                                                  >  1:7238654/1‑150 (MQ=255)
                     tGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAg                                                                                                                                  >  2:2002753/1‑150 (MQ=255)
                     tGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAg                                                                                                                                  >  2:4835875/1‑150 (MQ=255)
                            tCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCg                                                                                                                           >  1:8068914/1‑150 (MQ=255)
                               ccTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAt                                                                                                                        >  2:3072028/1‑150 (MQ=255)
                                 ttGTAGCAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTa                                                                                                                      >  1:1289879/1‑150 (MQ=255)
                                    tAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTc                                                                                                                   >  2:6457421/1‑150 (MQ=255)
                                         gCGGGTTCACCCGCGAACACGGGCACAGCCCGTGCCATTCACCCTCCAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCCTAACAAAACCCGCCCCAAGGGCGGGTTTTTTTACAGCTATCCGACTTACTCCggtc                                                                                                              <  2:8047092/150‑3 (MQ=255)
                                               tCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATTATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCt                                                                                                        <  1:7298881/150‑1 (MQ=255)
                                               tCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCt                                                                                                        <  1:6680410/150‑1 (MQ=255)
                                                   ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGATTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGc                                                                                                    <  2:1114360/150‑1 (MQ=255)
                                                    cgcgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCa                                                                                                   >  1:2921976/1‑150 (MQ=255)
                                                           acGGGCACCGCCCGTGCCATGCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCAGAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCa                                                                                            <  2:6995494/150‑1 (MQ=255)
                                                                 aaaGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCt                                                                                      >  1:1573662/1‑150 (MQ=255)
                                                                   aGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTc                                                                                    >  1:5340448/1‑150 (MQ=255)
                                                                    gCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCt                                                                                   <  2:6397082/150‑1 (MQ=255)
                                                                    gCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCt                                                                                   <  2:4681058/150‑1 (MQ=255)
                                                                     cccGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTg                                                                                  <  1:849320/150‑1 (MQ=255)
                                                                      ccGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTgg                                                                                 <  2:722196/150‑1 (MQ=255)
                                                                        gTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTAAAGCTATACGAATTACTCCGGACGAACCTGTGCAGCCTGCATGACCTTCTGGcc                                                                               >  2:7889477/1‑150 (MQ=255)
                                                                            cATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCgcgc                                                                           >  2:4272301/1‑150 (MQ=255)
                                                                             aTTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAAc                                                                                                          >  1:4822787/1‑118 (MQ=255)
                                                                              ttCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTc                                                                         >  1:3212492/1‑150 (MQ=255)
                                                                                   ccTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCCTAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcga                                                                    <  1:6910840/150‑1 (MQ=255)
                                                                                   ccTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcga                                                                    >  1:3341482/1‑150 (MQ=255)
                                                                                       cAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGaa                                                                <  2:3337104/150‑1 (MQ=255)
                                                                                        aGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAg                                                               >  1:7052795/1‑150 (MQ=255)
                                                                                          tCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACACGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGa                                                             >  2:7644861/1‑150 (MQ=255)
                                                                                              tttCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACa                                                         >  2:6109143/1‑150 (MQ=255)
                                                                                              tttCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACa                                                         <  1:6285366/150‑1 (MQ=255)
                                                                                                 cGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTc                                                      >  1:2333679/1‑150 (MQ=255)
                                                                                                   cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAGCAGTCTg                                                    >  1:4516967/1‑150 (MQ=255)
                                                                                                    aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTgg                                                   >  1:359840/1‑150 (MQ=255)
                                                                                                    aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTgg                                                   >  2:2120176/1‑150 (MQ=255)
                                                                                                    aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTgg                                                   >  2:3247521/1‑150 (MQ=255)
                                                                                                          tCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCcttct                                             >  1:8184367/1‑150 (MQ=255)
                                                                                                            ccGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCttt                                           >  1:3385668/1‑150 (MQ=255)
                                                                                                            ccGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCttt                                           >  1:8755422/1‑150 (MQ=255)
                                                                                                            ccGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCttt                                           >  1:6465802/1‑150 (MQ=255)
                                                                                                             cGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTc                                          >  1:6563553/1‑150 (MQ=255)
                                                                                                              gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGACACAGTCTGGCCTTCTTTCa                                         >  1:6727269/1‑150 (MQ=255)
                                                                                                              gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         <  1:2027094/150‑1 (MQ=255)
                                                                                                              gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTGTCa                                         >  2:3806753/1‑150 (MQ=255)
                                                                                                               ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAg                                        >  1:829227/1‑150 (MQ=255)
                                                                                                                cccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAgg                                       >  1:8097987/1‑150 (MQ=255)
                                                                                                                   ttCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGctc                                    >  2:1518258/1‑150 (MQ=255)
                                                                                                                    tCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGctct                                   <  2:6012076/150‑1 (MQ=255)
                                                                                                                        atgatgCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGaa                               >  1:4036926/1‑150 (MQ=255)
                                                                                                                         tgatgCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAg                              >  1:888714/1‑150 (MQ=255)
                                                                                                                            tgCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCg                           >  1:1677932/1‑150 (MQ=255)
                                                                                                                               aTCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTCGAAGCCGTCa                        >  1:2369111/1‑150 (MQ=255)
                                                                                                                                tCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAg                       >  1:7464866/1‑150 (MQ=255)
                                                                                                                                  aCCCATAAAAAATCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGAt                     >  2:4883979/1‑150 (MQ=255)
                                                                                                                                  aCCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGAt                     >  2:6718099/1‑150 (MQ=255)
                                                                                                                                  aCCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGAt                     >  2:168710/1‑150 (MQ=255)
                                                                                                                                   cccATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGAtt                    <  2:4453957/150‑1 (MQ=255)
                                                                                                                                    ccATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTc                   <  2:1289879/150‑1 (MQ=255)
                                                                                                                                    ccATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTc                   >  2:7420683/1‑150 (MQ=255)
                                                                                                                                      gtAAAAAATCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCAGAATTACTCCGGAAGAAACTGTGCAGCATGCATTACATTCTGTCCGCGATCGGCGACTAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGa                 >  2:8467938/2‑150 (MQ=255)
                                                                                                                                      aTAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGa                 <  2:7296427/150‑1 (MQ=255)
                                                                                                                                       tAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCGGGCTCTTGAAGCCGTCAGATTCGAt                <  1:2582137/150‑1 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATgg              <  1:3428416/150‑1 (MQ=255)
                                                                                                                                            aaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCTGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAGGCCGTCAGATTCGATGGCtt           >  2:5202955/1‑150 (MQ=255)
                                                                                                                                            aaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           >  1:3799421/1‑150 (MQ=255)
                                                                                                                                               cccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGa        >  2:4213678/1‑150 (MQ=255)
                                                                                                                                                ccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa       <  2:5904944/150‑1 (MQ=255)
                                                                                                                                                ccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa       >  2:555692/1‑150 (MQ=255)
                                                                                                                                                ccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa       <  2:8200729/150‑1 (MQ=255)
                                                                                                                                                  gCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAgt     <  2:8643959/150‑1 (MQ=255)
                                                                                                                                                   ccACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCg                  <  2:356548/136‑1 (MQ=255)
                                                                                                                                                   ccACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCg                  >  1:356548/1‑136 (MQ=255)
                                                                                                                                                     acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  <  2:6764797/150‑1 (MQ=255)
                                                                                                                                                     acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  <  1:5838782/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2813387‑2813685

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: