Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,272 G→C 10.1% intergenic (+48/+94) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,2720GC10.1% 434.3 / 31.9 139intergenic (+48/+94)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (57/68);  new base C (7/7);  total (64/75)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 145 total aligned reads displayed.

TGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGA  >  NC_002947/1707124‑1707421
                                                                                                                                                    |                                                                                                                                                      
tgTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcggg                                                                                                                                                       <  2:710968/150‑1 (MQ=255)
 gTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggc                                                                                                                                                      >  2:1699717/1‑150 (MQ=255)
  ttCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggcg                                                                                                                                                     >  2:4595274/1‑150 (MQ=255)
      acacCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                                 <  2:1981138/150‑1 (MQ=255)
       cacCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                                <  2:2394354/150‑1 (MQ=255)
       cacCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                                <  2:2394425/150‑1 (MQ=255)
        acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               <  2:149596/150‑1 (MQ=255)
        acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                               <  1:2288824/150‑1 (MQ=39)
          ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                                  >  1:5103313/1‑145 (MQ=255)
          ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                             >  1:4388648/1‑150 (MQ=255)
          ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                             <  1:3929860/150‑1 (MQ=255)
           cagcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                            >  1:4010136/1‑150 (MQ=255)
            agcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                           <  1:929667/150‑1 (MQ=255)
             gcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                          >  1:4719491/1‑150 (MQ=255)
             gcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                          >  1:5333110/1‑150 (MQ=255)
               agAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                        <  2:3729704/150‑1 (MQ=255)
                 aaaGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                      <  2:614322/150‑1 (MQ=255)
                     gCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcg                                                                                                                                  >  2:1787220/1‑150 (MQ=255)
                     gCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcg                                                                                                                                  >  2:2022876/1‑150 (MQ=255)
                     gCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcg                                                                                                                                  >  1:5143910/1‑150 (MQ=255)
                                  gCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGGTCGATCTCACAGGCGGCACAACATTATc                                                                                                                     >  1:5112333/1‑150 (MQ=255)
                                    tggtggACAACCGCTAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCTCAGGCGGCACAACATTATCgg                                                                                                                   >  1:3428856/1‑150 (MQ=255)
                                    tggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcgg                                                                                                                                 <  1:3584705/136‑1 (MQ=255)
                                    tggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcgg                                                                                                                                 >  2:3584705/1‑136 (MQ=255)
                                    tggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCgg                                                                                                                   >  1:3542535/1‑150 (MQ=255)
                                        ggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGaa                                                                                                               <  1:249230/150‑1 (MQ=255)
                                         gACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACCCCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAAc                                                                                                              >  2:2392834/1‑150 (MQ=255)
                                         gACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAAc                                                                                                              >  2:3675260/1‑150 (MQ=255)
                                           cAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaa                                                                                                            >  2:22635/1‑150 (MQ=255)
                                            aaCCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaaa                                                                                                           >  2:3389282/1‑150 (MQ=255)
                                            aaCCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaaa                                                                                                           <  2:3874210/150‑1 (MQ=255)
                                            aaCCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaaa                                                                                                           <  1:4595274/150‑1 (MQ=255)
                                              ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGa                                                                                                         >  1:4167136/1‑150 (MQ=255)
                                              ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGa                                                                                                         >  2:1402613/1‑150 (MQ=255)
                                              ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGa                                                                                                         >  2:4960427/1‑150 (MQ=255)
                                               cGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               >  1:5007211/1‑111 (MQ=255)
                                                gCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGaaa                                                                                                       <  1:1299756/150‑1 (MQ=255)
                                                 cAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGaaaa                                                                                                      >  2:3076209/1‑150 (MQ=255)
                                                 cAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGaaaa                                                                                                      <  2:4435611/150‑1 (MQ=255)
                                                 cAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGaaaa                                                                                                      <  1:1990578/150‑1 (MQ=255)
                                                   aGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGc                                                                                                    >  2:4983021/1‑150 (MQ=255)
                                                   aGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGc                                                                                                    >  2:4983089/1‑150 (MQ=255)
                                                    gcgcgcGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGcc                                                                                                   >  2:1359050/1‑150 (MQ=255)
                                                    gcgcgcGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGcc                                                                                                   >  1:4742605/1‑150 (MQ=255)
                                                           gTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTc                                                                                           <  2:4167136/150‑1 (MQ=255)
                                                               gcTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATc                                                                                       >  1:1586768/1‑150 (MQ=255)
                                                                     tGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgggg                                                                                 <  2:2406628/150‑1 (MQ=255)
                                                                        gCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCtt                                                                              >  1:1760625/1‑150 (MQ=255)
                                                                           aGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTc                                                                           <  1:3774857/150‑1 (MQ=255)
                                                                            ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACAGTATCGGCGAACAAAGAAAAGCCCCGCGATCGACGCGGGGCTTTTct                                                                          <  2:4682802/150‑2 (MQ=38)
                                                                             gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGACGCGGGGc                                                                                >  2:2806462/1‑143 (MQ=255)
                                                                             gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                         >  1:359012/1‑150 (MQ=255)
                                                                             gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                         >  2:1829559/1‑150 (MQ=255)
                                                                              aaaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCtt                                                                        <  2:5418848/150‑1 (MQ=255)
                                                                              aaaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCtt                                                                        >  2:2925604/1‑150 (MQ=255)
                                                                               aaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCttt                                                                       <  2:5264865/150‑1 (MQ=255)
                                                                               aaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCttt                                                                       <  1:1554530/150‑1 (MQ=255)
                                                                                 cacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCa                                                                     <  1:2922409/150‑1 (MQ=255)
                                                                                 cacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCa                                                                     <  2:2155525/150‑1 (MQ=255)
                                                                                  acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                          >  1:3221906/1‑81 (MQ=255)
                                                                                  acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                          <  2:3221906/81‑1 (MQ=255)
                                                                                  acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                                    >  1:4351532/1‑150 (MQ=255)
                                                                                  acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                                    <  1:137761/150‑1 (MQ=255)
                                                                                          aGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccag                                                            <  2:2616641/150‑1 (MQ=255)
                                                                                           ggTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                           <  1:719863/150‑1 (MQ=255)
                                                                                           ggTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                           <  1:719610/150‑1 (MQ=255)
                                                                                            gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagcc                                                          <  2:928116/150‑1 (MQ=255)
                                                                                               cAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTa                                                       <  1:1892304/150‑1 (MQ=255)
                                                                                               cAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTa                                                       <  1:3675260/150‑1 (MQ=255)
                                                                                                 gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCGTTAgg                                                     >  2:4663351/1‑149 (MQ=255)
                                                                                                 gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgc                                                     <  2:4508898/150‑1 (MQ=255)
                                                                                                  tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                   <  2:323736/119‑1 (MQ=255)
                                                                                                  tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                   >  1:323736/1‑119 (MQ=255)
                                                                                                    aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCt                                                  >  2:3670528/1‑150 (MQ=255)
                                                                                                       tttCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGaa                                               <  2:2727563/150‑1 (MQ=255)
                                                                                                            gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                          <  1:5376901/150‑1 (MQ=255)
                                                                                                            gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                          >  2:3416289/1‑150 (MQ=255)
                                                                                                            gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGAGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTCGCCAGCCTTAGCGCTGAACCTgg                                          >  1:4178046/1‑150 (MQ=255)
                                                                                                                 gcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAc                                     >  2:2794695/1‑150 (MQ=255)
                                                                                                                   aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCg                                   >  1:4286404/1‑150 (MQ=255)
                                                                                                                   aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCg                                   <  1:2794733/150‑1 (MQ=255)
                                                                                                                   aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCg                                   <  1:2794695/150‑1 (MQ=255)
                                                                                                                    attttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGa                                  <  2:1307140/149‑1 (MQ=255)
                                                                                                                        tGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgt                              <  1:3819965/150‑1 (MQ=255)
                                                                                                                         ggCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtg                             <  1:3325673/150‑1 (MQ=255)
                                                                                                                           cgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                           <  1:3670531/150‑1 (MQ=255)
                                                                                                                                 cgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                          >  1:804004/1‑129 (MQ=255)
                                                                                                                                       cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTgg               <  2:827765/150‑1 (MQ=255)
                                                                                                                                       cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTgg               >  2:3056392/1‑150 (MQ=255)
                                                                                                                                       cGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                           <  1:4858116/122‑1 (MQ=255)
                                                                                                                                       cGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                           >  2:4858116/1‑122 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCGGGt              >  2:4830936/1‑150 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt              >  1:4459744/1‑150 (MQ=255)
                                                                                                                                         aGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg                >  2:3015032/1‑147 (MQ=255)
                                                                                                                                         aGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg                <  1:3015032/147‑1 (MQ=255)
                                                                                                                                            gccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTc          <  1:2171976/150‑1 (MQ=255)
                                                                                                                                            gccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTc          <  1:1699717/150‑1 (MQ=255)
                                                                                                                                             ccgcccgcGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCg         >  2:5009312/1‑150 (MQ=255)
                                                                                                                                               gcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGACCACAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt       <  2:5017455/150‑1 (MQ=255)
                                                                                                                                                    ggcggcGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGgaga  <  1:1137747/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGA  >  NC_002947/1707124‑1707421

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