Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,837 G→T 6.8% noncoding (433/451 nt) PP_mr14 ← cspA regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,261,8370GT 6.8% 989.3 / 44.4 296noncoding (433/451 nt)PP_mr14cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base G (143/133);  new base T (12/8);  total (155/141)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 298 total aligned reads displayed.

CATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTT  >  NC_002947/1261693‑1261986
                                                                                                                                                |                                                                                                                                                     
cATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGc                                                                                                                                                  >  2:2623908/1‑150 (MQ=255)
     gCGAAGCCTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtt                                                                                                                                             <  1:1111195/150‑1 (MQ=38)
      cGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttt                                                                                                                                            <  2:2377161/150‑1 (MQ=255)
      cGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttt                                                                                                                                            >  2:4168084/1‑150 (MQ=255)
      cGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttt                                                                                                                                            >  1:2294956/1‑150 (MQ=255)
       gAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGtttt                                                                                                                                           >  2:2904278/1‑150 (MQ=255)
       gAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGtttt                                                                                                                                           >  1:2467780/1‑150 (MQ=255)
        aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttt                                                                                                                                          <  1:5165407/150‑1 (MQ=255)
        aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttt                                                                                                                                          <  1:2610824/150‑1 (MQ=255)
        aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                                                                                                                          <  1:3459306/150‑1 (MQ=255)
             tGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaa                                                                                                                                     >  2:3196396/1‑150 (MQ=255)
             tGACCTGCCAGACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaa                                                                                                                                     <  2:4634897/150‑1 (MQ=255)
              gACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaaa                                                                                                                                    <  2:2465498/150‑1 (MQ=255)
                  tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAATCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                >  2:3062184/1‑150 (MQ=255)
                  tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                >  2:568619/1‑150 (MQ=255)
                  tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                >  1:1870020/1‑150 (MQ=255)
                  tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                >  1:1788434/1‑150 (MQ=255)
                  tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                >  1:829158/1‑150 (MQ=255)
                          cgcgTTCGCCCGCACACAGCGACTGCGCCGTACTTGACGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTCTATCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTa                                                                                                                        <  2:3042894/147‑1 (MQ=39)
                          cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                                                                                                                        >  1:4047867/1‑134 (MQ=255)
                             gTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGCGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGAt                                                                                                                     <  1:2157829/150‑1 (MQ=255)
                                 gCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGAc                                                                                                                 >  2:4051886/1‑150 (MQ=255)
                                 ccgctccCACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGAc                                                                                                                 <  1:1706530‑M1/143‑1 (MQ=255)
                                   ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                               >  1:3239409/1‑150 (MQ=255)
                                   ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                               >  1:1970880/1‑150 (MQ=255)
                                   ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                               >  2:1647598/1‑150 (MQ=255)
                                    cGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAATGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTa                                                                                                              <  1:5342758/150‑1 (MQ=255)
                                     gCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAg                                                                                                             <  1:4520525/150‑1 (MQ=255)
                                     gCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAg                                                                                                             <  2:4205944/150‑1 (MQ=255)
                                     gCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAg                                                                                                             <  1:3700445/150‑1 (MQ=255)
                                       acacaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGc                                                                                                           >  2:3621957/1‑150 (MQ=255)
                                          CAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGAt                                                                                                        >  1:2004415‑M1/1‑150 (MQ=255)
                                           AGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATc                                                                                                       <  2:2141589‑M1/150‑1 (MQ=255)
                                           AGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATc                                                                                                       >  2:2740497‑M1/1‑150 (MQ=255)
                                            GCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCt                                                                                                      <  2:715564‑M1/150‑1 (MQ=255)
                                             cgagggcgcCGTCCTTGATGCCTTGGCTTTCCCTTGGGGTGGGGGGTAGCCCGGGATGTGGCCTGCGCCTATTGCGTGCAGGCATAAAAAAACCCGCCAGATGGGGGTTTTTTTTGAAAGCTAGGAAGCTAGATCGCCTAAGTCGATCta                                                                                                     <  1:5090744/145‑2 (MQ=255)
                                             cGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTg                                                                                                     >  2:5507495/1‑150 (MQ=255)
                                              gACTGCGCCGTACTTGAAGCCTGTGCTGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGa                                                                                                    <  2:168239/150‑1 (MQ=255)
                                              gACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCCAGATCGCCTAAGTCGATCtag                                                                                                    <  1:850593/150‑3 (MQ=255)
                                               aCTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGaa                                                                                                   <  1:3337853/150‑1 (MQ=255)
                                                 tGCGCCGTACTTGAAGCCTATGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTAAGGCGATCTGAAct                                                                                                 >  2:4001386/1‑150 (MQ=255)
                                                    gcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGAt                                                                                                        >  2:3281739/1‑140 (MQ=255)
                                                     ccGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtctt                                                                                             <  2:1706530/150‑1 (MQ=255)
                                                      cGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACttcttc                                                                                            >  1:1308926/1‑150 (MQ=255)
                                                       gTACTGGCAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCa                                                                                           <  2:3747259/150‑1 (MQ=255)
                                                         aCTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGc                                                                                         <  1:1585089/150‑1 (MQ=255)
                                                           ttGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCt                                                                                       >  1:1017951/1‑150 (MQ=255)
                                                            tGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTg                                                                                      >  1:2624604/1‑150 (MQ=255)
                                                             gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGc                                                                                     >  2:1144714/1‑150 (MQ=255)
                                                              aaGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCa                                                                                    <  2:2201396/150‑1 (MQ=255)
                                                                gCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATg                                                                                  >  2:4966476/1‑150 (MQ=255)
                                                                   tgtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCt                                                                               <  1:2904278/150‑1 (MQ=255)
                                                                      gCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTc                                                                            >  2:5601011/1‑150 (MQ=255)
                                                                      gCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTc                                                                            >  2:123545/1‑150 (MQ=255)
                                                                      gCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTc                                                                            <  1:2200997/150‑1 (MQ=255)
                                                                        aGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                          >  1:3810776/1‑150 (MQ=255)
                                                                        aGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                          >  2:1028160/1‑150 (MQ=255)
                                                                        aGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                          >  1:2977685/1‑150 (MQ=255)
                                                                          tCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGc                                                                        >  1:784848/1‑150 (MQ=255)
                                                                            ccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcg                                                                      >  2:474157/1‑150 (MQ=255)
                                                                            ccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcg                                                                      >  1:37575/1‑150 (MQ=255)
                                                                              tgtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgcg                                                                    <  2:4482881‑M2/146‑1 (MQ=255)
                                                                                tcttcgCGGGCTCGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCggg                                                                  >  1:5212844/6‑150 (MQ=255)
                                                                                   gAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGACGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAg                                                               <  2:3453559/150‑1 (MQ=255)
                                                                                       gggctcgCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgac                                                           >  1:4254546/7‑150 (MQ=255)
                                                                                         gCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgacga                                                         >  1:3262392/1‑150 (MQ=255)
                                                                                         gCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgacga                                                         >  1:415760/1‑150 (MQ=255)
                                                                                           gAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGaaa                                                       <  2:5581030/150‑1 (MQ=255)
                                                                                               ccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                   <  1:5444706/150‑1 (MQ=255)
                                                                                                cgcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAAc                                                  >  2:443366/1‑150 (MQ=255)
                                                                                                 gcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACa                                                 <  1:2749257/150‑1 (MQ=255)
                                                                                                 gcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTGTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCGTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACa                                                 >  2:5005445/1‑150 (MQ=255)
                                                                                                  cgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCATTTCTGGCCGCGGGTAGCGACGAAAGTAACAg                                                <  2:2450829/150‑1 (MQ=255)
                                                                                                      gagGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTg                                            <  2:5114033/150‑1 (MQ=255)
                                                                                                         gCCCGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGa                                                                                                                  >  1:3112464/1‑77 (MQ=255)
                                                                                                         gCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGAGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGa                                                                                                                  <  2:3112464/77‑1 (MQ=255)
                                                                                                           cAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCctt                                       >  1:1010147/1‑150 (MQ=255)
                                                                                                             gcgcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCcttct                                     <  2:247837/150‑1 (MQ=255)
                                                                                                               gcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCttt                                   >  1:4329015/1‑150 (MQ=255)
                                                                                                                 cTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCa                                 <  2:637769/150‑1 (MQ=255)
                                                                                                                 cTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCa                                 <  2:1436976/150‑1 (MQ=255)
                                                                                                                 cTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCa                                 >  1:1943901/1‑150 (MQ=255)
                                                                                                                 cTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCa                                 <  2:3926269/150‑1 (MQ=255)
                                                                                                                   aTTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                               <  2:1870020/150‑1 (MQ=255)
                                                                                                                    tttCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAggg                              >  1:4557809/1‑150 (MQ=255)
                                                                                                                      tCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTCGCGACGAAGGTAACAGCCTGGCCTTCTTTCGGGGtt                            >  2:5314607/1‑150 (MQ=255)
                                                                                                                      tCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGtt                            <  1:3883613/150‑1 (MQ=255)
                                                                                                                              tttcaaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                    >  1:3608535/5‑150 (MQ=255)
                                                                                                                              tttcaaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                    >  1:1031934/5‑150 (MQ=255)
                                                                                                                              gCATACAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                    <  1:4630120/150‑1 (MQ=255)
                                                                                                                              gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                    >  1:4906225/1‑150 (MQ=255)
                                                                                                                                 tAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGAAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTc                 >  1:4187333/1‑150 (MQ=255)
                                                                                                                                  aaaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCa                >  1:1074711/1‑150 (MQ=255)
                                                                                                                                   aaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAg               >  1:4824529/1‑150 (MQ=255)
                                                                                                                                      aaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCtt            <  1:2517942/150‑1 (MQ=255)
                                                                                                                                        aCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTgg          <  2:197533/150‑1 (MQ=255)
                                                                                                                                        aCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTgg          >  2:1951940/1‑150 (MQ=255)
                                                                                                                                         cccGCCAGAAGGCGGGTTTTTTTGAACGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGa         <  2:825959/150‑1 (MQ=255)
                                                                                                                                           cGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATg       >  1:3212316/1‑150 (MQ=255)
                                                                                                                                              cAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCt    >  2:4826164/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTT  >  NC_002947/1261693‑1261986

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: ‑