Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,039 C→T 5.1% R403C (CGC→TGC)  imuB → DNA repair DNA linked protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,528,0390CT 5.1% 296.1 / 4.1 97R403C (CGC→TGC) imuBDNA repair DNA linked protein
Reads supporting (aligned to +/- strand):  ref base C (45/47);  new base T (2/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGGCCTGC  >  NC_002947/3527891‑3528186
                                                                                                                                                    |                                                                                                                                                   
cTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGAGCAGCCGCGCGCACAAGGCAACATGCCGGTTGCACCGGGCAgccg                                                                                                                                                    >  2:2070355/1‑150 (MQ=255)
    gcgcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCg                                                                                                                                                >  2:166081/1‑150 (MQ=255)
      gcAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCTGCCCggg                                                                                                                                              >  2:1374569/1‑150 (MQ=255)
      gcAATTGCGCGAGCGGTTGCGCGCCCGCCGGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCGGCTGCCCggt                                                                                                                                              <  1:1704546/150‑2 (MQ=255)
      cccactGCCCGAGCGGTTGCGCGCCTCCCCGCGGAATGAGGCCCTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGAGCAGCCGCCCggg                                                                                                                                              <  2:797444/145‑1 (MQ=255)
       cAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCgggg                                                                                                                                             <  1:599109/150‑1 (MQ=255)
        aaTTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGACTGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCTGCCCGGGGt                                                                                                                                            <  1:1374569/150‑1 (MQ=255)
             cgcgAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtg                                                                                                                                       <  2:969992/150‑1 (MQ=255)
                  gCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTc                                                                                                                                  <  1:228382/150‑1 (MQ=255)
                           cgcccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGAt                                                                                                                         <  1:2019156/150‑1 (MQ=255)
                              ccgccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCATGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATcgc                                                                                                                      >  2:848500/1‑150 (MQ=255)
                                 ccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCt                                                                                                                   <  1:2345078/150‑1 (MQ=255)
                                 ccTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCt                                                                                                                   <  1:166081/150‑1 (MQ=255)
                                     ggggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGCGCTTTTTGCCCAAGGCAACATGCAGGATCCACCGGGCTGTCGTTCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGAt                                                                                                               >  1:1739592/1‑150 (MQ=255)
                                       ggATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGa                                                                                                             >  2:2371283/1‑150 (MQ=255)
                                             ggCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCAGAAGGCGCCTGGCATTTGGCGGAGCCGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGCcgccgc                                                                                                       >  1:562068/1‑150 (MQ=255)
                                              gCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACgccgcc                                                                                                      >  1:2256031/1‑150 (MQ=255)
                                              gCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACgccgcc                                                                                                      >  1:690338/1‑150 (MQ=255)
                                              gCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACgccgcc                                                                                                      >  1:1274261/1‑150 (MQ=255)
                                                   cAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCATCGCAGGCCAGAGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCa                                                                                                 <  2:1871607/150‑1 (MQ=255)
                                                   cAAGGGGCTTGGCGCACCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCGGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCa                                                                                                 <  1:1343584/150‑1 (MQ=255)
                                                      ggggCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgtt                                                                                                                                     <  2:2197810/111‑1 (MQ=255)
                                                      ggggCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCtgtt                                                                                                                                     >  1:2197810/1‑111 (MQ=255)
                                                      ggggCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGt                                                                                              <  2:53654/150‑1 (MQ=255)
                                                      cacacTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGACGCCGCCCAAGGCGACTGATGTCAGCACAGCTCGCACGGGGTGTCGCGGCCGACCAgga                                                                                              >  1:1892293/5‑149 (MQ=255)
                                                           ttGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAgg                                                                                         <  2:2345582/150‑1 (MQ=255)
                                                           ttGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAgg                                                                                         <  2:690338/150‑1 (MQ=255)
                                                               cgcAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTa                                                                                     >  1:483366/1‑150 (MQ=255)
                                                                gcAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTAt                                                                                    >  1:1114334/1‑150 (MQ=255)
                                                                   gccgccGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGcc                                                                                 <  1:1655052/150‑1 (MQ=255)
                                                                      gccgATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCCGGTGCAGGGCTATGCCGaa                                                                              >  1:1061312/1‑150 (MQ=255)
                                                                       ccgATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAAc                                                                             <  2:1274261/150‑1 (MQ=255)
                                                                         gATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGc                                                                           >  2:2161763/1‑150 (MQ=255)
                                                                         gATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGATGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGAt                                                                                                               <  1:1959579/114‑1 (MQ=255)
                                                                         gATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGATGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGAt                                                                                                               >  2:1959579/1‑114 (MQ=255)
                                                                            cACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGCCGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATc                                                                        <  1:1191335/150‑1 (MQ=255)
                                                                             aCCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCg                                                                       >  2:994100/1‑150 (MQ=255)
                                                                             aCCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGAACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCg                                                                       >  2:1645563/1‑150 (MQ=255)
                                                                              ccGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCa                                                                                                                              <  2:715168/94‑1 (MQ=255)
                                                                              ccGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCa                                                                                                                              >  1:715168/1‑94 (MQ=255)
                                                                               cGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGaa                                                                     <  2:483366/150‑1 (MQ=255)
                                                                                 ccccGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAg                                                                                          >  2:1227854/1‑127 (MQ=255)
                                                                                 ccccGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAg                                                                                          <  1:1227854/127‑1 (MQ=255)
                                                                                 ccccGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATc                                                                   >  1:1486592/1‑150 (MQ=255)
                                                                                   ccGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGAAGGGCTATGCCGAACGCATCGAATCAg                                                                 <  1:2371283/150‑1 (MQ=255)
                                                                                    cGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAgg                                                                >  2:2340949/1‑150 (MQ=255)
                                                                                      aaTGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGtt                                                              <  2:2256031/150‑1 (MQ=255)
                                                                                       aTGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtg                                                             <  1:2161763/150‑1 (MQ=255)
                                                                                        tGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAg                                                                                          <  2:813960/120‑1 (MQ=255)
                                                                                        tGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAg                                                                                          >  1:813960/1‑120 (MQ=255)
                                                                                          cgcCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtg                                                          <  2:96734/150‑1 (MQ=255)
                                                                                          cgcCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtg                                                          <  2:1486592/150‑1 (MQ=255)
                                                                                           gcCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACgccgcc                                                                                                      <  2:1532553/105‑1 (MQ=255)
                                                                                           gcCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACgccgcc                                                                                                      >  1:1532553/1‑105 (MQ=255)
                                                                                             cTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGa                                                       >  2:1516980/1‑150 (MQ=255)
                                                                                              tGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCg                                                                                                                                 >  2:1539114/1‑75 (MQ=255)
                                                                                              tGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCg                                                                                                                                 <  1:1539114/75‑1 (MQ=255)
                                                                                              tGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCAc                                                                                                                             >  1:343159/1‑79 (MQ=255)
                                                                                              tGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCAc                                                                                                                             <  2:343159/79‑1 (MQ=255)
                                                                                               ggCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATg                                                     >  1:1094754/1‑150 (MQ=255)
                                                                                                  aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAAc                                                                             >  2:1393867/1‑123 (MQ=255)
                                                                                                  aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAAc                                                                             <  1:1393867/123‑1 (MQ=255)
                                                                                                  aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATggcg                                                  >  2:132315/1‑150 (MQ=255)
                                                                                                  aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATggcg                                                  >  2:425384/1‑150 (MQ=255)
                                                                                                  aTTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATggcg                                                  >  1:240713/1‑150 (MQ=255)
                                                                                                   tttGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATggcgg                                                 >  2:2042918/1‑150 (MQ=255)
                                                                                                      ggCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGa                                              <  2:1060476/150‑1 (MQ=255)
                                                                                                      ggCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGa                                              <  1:135376/150‑1 (MQ=255)
                                                                                                         ggAGCAGGGTGCCCAAGGCAACATGCGCTTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAtgt                                           <  2:562068/150‑1 (MQ=255)
                                                                                                         ggAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAtgt                                           <  1:311815/150‑1 (MQ=255)
                                                                                                           aGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGc                                         <  1:2042918/150‑1 (MQ=255)
                                                                                                              aGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTc                                      >  1:316870/1‑150 (MQ=255)
                                                                                                               gggTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTcg                                     >  1:445337/1‑150 (MQ=255)
                                                                                                                ggTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTcgc                                    <  2:2258392/150‑1 (MQ=255)
                                                                                                                ggTGCCCAAGGCAACATGCCGGTTGCACCGAGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAc                                                                                                   >  2:2039313/1‑87 (MQ=255)
                                                                                                                ggTGCCCAAGGCAACATGCCGGTTGCACCGAGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTAc                                                                                                   <  1:2039313/87‑1 (MQ=255)
                                                                                                                  tGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGa                                  <  2:1114334/150‑1 (MQ=255)
                                                                                                                     ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGacta                               >  1:944098/1‑150 (MQ=255)
                                                                                                                     ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGacta                               <  1:994100/150‑1 (MQ=255)
                                                                                                                     ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGacta                               >  1:651813/1‑150 (MQ=255)
                                                                                                                     ccAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGacta                               >  1:829282/1‑150 (MQ=255)
                                                                                                                      cAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGactac                              >  2:2241999/1‑150 (MQ=255)
                                                                                                                       aaGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGactact                             <  1:2070355/150‑1 (MQ=255)
                                                                                                                                tGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGa                                  >  1:177233/1‑136 (MQ=255)
                                                                                                                                tGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGa                                  <  2:177233/136‑1 (MQ=255)
                                                                                                                                    ggTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAAc                <  1:848500/150‑1 (MQ=255)
                                                                                                                                        gCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGt            <  1:1516980/150‑1 (MQ=255)
                                                                                                                                        accccGGGCGGCTGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGaa                                                                              <  1:2129531/81‑1 (MQ=255)
                                                                                                                                        accccGGGCGGCTGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGaa                                                                              >  2:2129531/4‑84 (MQ=255)
                                                                                                                                          aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtgg                                                         <  1:1273142/103‑1 (MQ=255)
                                                                                                                                          aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTtggtgg                                                         >  2:1273142/1‑103 (MQ=255)
                                                                                                                                          aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGa          <  2:651813/150‑1 (MQ=255)
                                                                                                                                          aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGa          <  2:445337/150‑1 (MQ=255)
                                                                                                                                          aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGa          >  2:1351824/1‑150 (MQ=255)
                                                                                                                                          aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGa          >  1:557909/1‑150 (MQ=255)
                                                                                                                                          aCCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGa          >  1:1002723/1‑150 (MQ=255)
                                                                                                                                             caacAGCCACCCCGGGCGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAt                                             >  2:1574527/4‑112 (MQ=255)
                                                                                                                                             caacAGCCACCCCGGGCGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGAt                                             <  1:1574527/109‑1 (MQ=255)
                                                                                                                                                  gccgccCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGGCCTgc  <  1:359538/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGGGCGCAATTGCGCGAGCGGTTGCGCGCCCGCCTGGGGGATGAGGCCGTCAAGGGGCTTGGCGCAGCCGCCGATCACCGCCCCGAATGCGCCTGGCATTTGGCGGAGCAGGGTGCCCAAGGCAACATGCCGGTTGCACCGGGCAGCCGCCCGGGGTGGCTGTTGTCGGCACCGATCGCGCTGGATGACGCCGCCTACCAGGTGCAGGGCTATGCCGAACGCATCGAATCAGGTTGGTGGGATGGCGGCGATGTGCGTCGCGACTACTACCGCATCGAAACCCGTGACGGCCTGC  >  NC_002947/3527891‑3528186

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