Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,525,527 T→C 10.2% R89R (CGT→CGC)  PP_4858 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,525,5270TC10.2% 160.1 / 4.3 74R89R (CGT→CGC) PP_4858hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base C (2/6);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAGGCGG  >  NC_002947/5525379‑5525676
                                                                                                                                                    |                                                                                                                                                     
aGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGCTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGGTCGAACGTTCTGGCCgtg                                                                                                                                                      >  2:1828328/1‑150 (MQ=255)
    ccTGACCACCGTAGCTGGAGGCGGCTCCGATCGACTGATCGAACGTTCTGTCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGACGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGc                                                                                                                                                  <  2:1047391/149‑1 (MQ=14)
       gAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGa                                                                                                                                               >  1:953333/1‑150 (MQ=255)
        aaCACCGTAGCTGAAGGCGGCTCCGCTGGCCTGGTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGaa                                                                                                                                              <  2:1934826/150‑1 (MQ=255)
        aaCACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGaa                                                                                                                                              >  1:1588696/1‑150 (MQ=255)
         acacCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAg                                                                                                                                             <  2:1588696/150‑1 (MQ=255)
          cacCGTAGCTGAATGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAgg                                                                                                                                            <  2:1875160/150‑1 (MQ=255)
           acCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                                                                                           <  2:2142067/150‑1 (MQ=255)
             cGTAGCTGGAGGCGGCGCGGCGCGGCTGGTCGAACGGTCGGGCAGCGTAGCGGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTgg                                                                                                                                         <  2:364436/150‑1 (MQ=14)
             cGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTgg                                                                                                                                         <  1:2074858/150‑1 (MQ=255)
                aGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTc                                                                                                                                      <  1:1644836/150‑1 (MQ=255)
                     aaGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCCGAAGGTGGCTCCGATCGTCTGATCGAACGTGCTGGCCGTGTTGCAGAAGGCGGCTCCggat                                                                                                                                 >  2:1219599/1‑147 (MQ=11)
                      aGGCGGCTCCGATCGCCTGATCGAACGTTCTGCCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCg                                                                                                                                >  1:220071/1‑150 (MQ=255)
                        gcggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGTCGCGTTGCTGAAGGTGTCTCCGTTCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTc                                                                                                                              <  1:293379/150‑1 (MQ=255)
                        gcggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGCTCGTc                                                                                                                              >  2:1554785/1‑150 (MQ=38)
                           gcTCCGATCGCATGATCGAACGTGCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGa                                                                                                                           <  1:1075904/150‑1 (MQ=255)
                                       tGATCCAACCTCCGGGCCGCGTTGCGGAAGGCGGCTCCGATCGCCTGATCGAGCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTg                                                                                                               <  2:1432683/150‑1 (MQ=11)
                                       ccgccGAACGTCCGGGCCCCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTg                                                                                                               <  2:1351129/146‑1 (MQ=255)
                                            gagcGTTCGGGCGGCGTCGCTGACGGTGGCGCCGACGGCCTGATCGAACGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCcgc                                                                                                          <  2:2196990/147‑1 (MQ=11)
                                             aaCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                                                                                                 <  1:1689315/126‑1 (MQ=40)
                                             aaCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                                                                                                 >  2:1689315/1‑126 (MQ=39)
                                                     ggCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                                                                                           <  1:1436098/108‑1 (MQ=255)
                                                     ggCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                                                                                           >  2:1436098/1‑108 (MQ=40)
                                                      gCCGCGTTGGTGAAGGTGGCTCCGATCGCCGGGTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGGTGAAGGTGGCGCCGATCGTCTGCGCGAACGTTGGGGCCGGGTTGCTGAAg                                                                                                >  2:1970472/1‑150 (MQ=11)
                                                      gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATGGCTTGCTCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAg                                                                                                >  1:1992389/1‑150 (MQ=255)
                                                       ccGCGTTTCTGAAGGTGGCGCCGAGGGCCCGCTCGCACGTTCTGGCAGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAgg                                                                                               <  2:953333/150‑1 (MQ=14)
                                                       ccGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTGTGGCCGGGTGGcgggggg                                                                                               >  2:1742307/1‑144 (MQ=11)
                                                                tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATTGAACGTTCTGGTGGGGTTGCTGAAGGTTGCTCCGATCGTCTCATTTCTTATACTGGCCGTGTTGCTGACGGTGGCCCCgc                                                                                      >  1:1438478/1‑149 (MQ=11)
                                                                tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCGGAGGGTGGCTCCgg                                                                                      >  2:2020040/1‑149 (MQ=255)
                                                                    ggcggCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTg                                                                                                               <  2:328373/118‑1 (MQ=21)
                                                                    ggTGGCTCTGATCGACTGGTCGAACGTTCTGGCAGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGt                                                                                  <  2:954824/150‑1 (MQ=11)
                                                                       ggCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTg                                                                               >  2:1703546/1‑150 (MQ=14)
                                                                        gCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGGGTTGCCGAAGGTGGCTCCGATCGTCTGATCGAACACTCTGGCCGCGTTGCGGGAGGCGGCTCGGATCGTCTGa                                                                              >  2:377981/1‑150 (MQ=11)
                                                                         cTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAg                                                                                                                                             >  2:1979494/1‑86 (MQ=35)
                                                                         cTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAg                                                                                                                                             <  1:1979494/86‑1 (MQ=35)
                                                                          tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTg                                                                                                                                                <  2:1948748/82‑1 (MQ=35)
                                                                          tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTg                                                                                                                                                >  1:1948748/1‑82 (MQ=35)
                                                                          tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                                                                                           <  1:1026374/87‑1 (MQ=35)
                                                                          tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                                                                                           >  2:1026374/1‑87 (MQ=33)
                                                                          tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGACCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTGGTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCtagaa                                                                            >  2:13380/1‑146 (MQ=11)
                                                                                 gCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGtt                                                                     <  1:1703546/150‑1 (MQ=17)
                                                                                   cTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGGTCGAACGTTCTGGCCGCGTGGCCGGGGGGGGCTCCGATCGTCTGATAGATCGGTCTGGCCGCGTGGCGGGAGGGAGCTCAGATGGTCGGatcgatcgcggt                                                                   >  2:2104680/1‑146 (MQ=2)
                                                                                       tCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGcc                                                               >  1:1431036/1‑150 (MQ=38)
                                                                                        cGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCg                                                              >  1:1103905/1‑150 (MQ=38)
                                                                                         gAACGTTCTGGCCGCGCTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGa                                                                                                                           <  1:2241318/88‑1 (MQ=255)
                                                                                         gAACGTTCTGGCCGCGCTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGa                                                                                                                           >  2:2241318/1‑88 (MQ=255)
                                                                                              ttCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                                                                                                 >  1:2110763/1‑77 (MQ=35)
                                                                                              ttCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                                                                                                 <  2:2110763/77‑1 (MQ=35)
                                                                                                  ggCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGaa                                                                                                                      <  1:2281709/84‑1 (MQ=39)
                                                                                                     cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGt                                                 <  1:2412138/150‑1 (MQ=255)
                                                                                                            cTGAAGGTGGCTCCGGTCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCCGAAGGTGGCTCCGATCGTCAGCTCGAACGTGCCGGCCGTGTTGCCGAAGGAGGCTccg                                          >  2:223893/1‑148 (MQ=11)
                                                                                                            cTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCgg                                          <  1:442524/150‑1 (MQ=255)
                                                                                                             tGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGa                                         <  2:340651/150‑1 (MQ=21)
                                                                                                                 ggTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCCGAAGGTGGCTCGGATCGt                                     <  2:1625872/150‑1 (MQ=255)
                                                                                                                  gTGGCTCCTATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTGCGTGCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTc                                    <  1:1970472/150‑1 (MQ=255)
                                                                                                                    ggCTCTGATGGCGTGATCGAACGTTCTGTGAGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTg                                  <  2:2357907/150‑1 (MQ=14)
                                                                                                                     ccgcagATCGCCTGATCGAACGTGCTGGCAGTGTTGATGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGa                                 <  1:2104680/145‑1 (MQ=255)
                                                                                                                       tCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCGGATCGTCTGATCGAACGCTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGCCTGATc                               >  2:451713/1‑150 (MQ=255)
                                                                                                                       tCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGt                                                                                              <  1:577643/87‑1 (MQ=40)
                                                                                                                       tCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGt                                                                                              >  2:577643/1‑87 (MQ=39)
                                                                                                                         cGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTc                                                                                         <  1:2379754/90‑1 (MQ=40)
                                                                                                                            tCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGATTCCACAGATTCACGGATCGGTTTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACa                          <  1:1742307/150‑1 (MQ=11)
                                                                                                                               ccTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                                                                                         >  2:837316/1‑52 (MQ=255)
                                                                                                                               ccTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGc                       >  1:1277441/1‑150 (MQ=255)
                                                                                                                               ccTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGc                       >  1:1423979/1‑150 (MQ=255)
                                                                                                                                cTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCt                                                                   <  1:2203019/105‑1 (MQ=17)
                                                                                                                                  gATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTgg                    <  1:2020040/150‑1 (MQ=255)
                                                                                                                                      gAACGTTCTGGCCGTGTTGCTGAAGGTTGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATc                                       <  2:1254166/127‑1 (MQ=255)
                                                                                                                                      gAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATc                                       >  1:1254166/1‑127 (MQ=255)
                                                                                                                                         cctccgggCCCTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGtt             <  2:2280358/144‑1 (MQ=255)
                                                                                                                                          gTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGaa                                                    >  2:246795/1‑110 (MQ=255)
                                                                                                                                          gTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGaa                                                    <  1:246795/110‑1 (MQ=255)
                                                                                                                                               ggCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTg                                                <  2:298495/109‑1 (MQ=255)
                                                                                                                                               ggCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTg                                                >  1:298495/1‑109 (MQ=255)
                                                                                                                                                  cGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACACAATGACCGCGTTGCTGAAGGTGGCCCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAggc    >  2:1902102/1‑150 (MQ=17)
                                                                                                                                                    tgtTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACaa                         >  1:1617529/1‑127 (MQ=255)
                                                                                                                                                    tgtTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACaa                         <  2:1617529/127‑1 (MQ=255)
                                                                                                                                                    tgtTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAggcgg  <  1:1554785/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAGGCGG  >  NC_002947/5525379‑5525676

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