Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,800 G→A 5.5% G213G (GGC→GGT PP_3503 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,8000GA 5.5% 512.6 / 12.7 163G213G (GGC→GGTPP_3503sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (79/75);  new base A (7/2);  total (86/77)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 169 total aligned reads displayed.

CCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCGCTCTCCCCGCGG  >  NC_002947/3975651‑3975947
                                                                                                                                                     |                                                                                                                                                   
cccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcg                                                                                                                                                     >  2:1850655/1‑150 (MQ=255)
cccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcg                                                                                                                                                     >  2:2538668/1‑150 (MQ=255)
 ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                                                                                                    >  1:2239164/1‑150 (MQ=255)
 ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                                                                                                    >  1:1229768/1‑150 (MQ=255)
    cccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCt                                                                                                                                                 >  1:1811512/1‑150 (MQ=255)
    cccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCt                                                                                                                                                 >  2:3096815/1‑150 (MQ=255)
     ccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCtt                                                                                                                                                <  2:4973011/150‑1 (MQ=255)
       gCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtct                                                                                                                                              >  2:3485753/1‑150 (MQ=255)
         tCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCa                                                                                                                                            >  2:2124629/1‑150 (MQ=255)
          cGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCAt                                                                                                                                           <  1:347065/150‑1 (MQ=255)
           gATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATg                                                                                                                                          >  1:4566574/1‑150 (MQ=255)
           gATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATg                                                                                                                                          <  2:1481773/150‑1 (MQ=255)
             tgtgCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGc                                                                                                                                        >  2:4706289/1‑150 (MQ=255)
                gCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGa                                                                                                                                     >  1:140543/1‑150 (MQ=255)
                 cTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGaa                                                                                                                                    <  1:3562036/150‑1 (MQ=255)
                   tctcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCGAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACa                                                                                                                                  <  2:4393121/150‑1 (MQ=255)
                     tcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGt                                                                                                                                >  1:856875/1‑150 (MQ=255)
                      ccTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGtt                                                                                                                               >  2:3756960/1‑150 (MQ=255)
                        tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                                                                                                             >  1:782007/1‑150 (MQ=255)
                        tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                                                                                                             <  2:4111686/150‑1 (MQ=255)
                        tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                                                                                                             >  1:2638531/1‑150 (MQ=255)
                        tGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCg                                                                                                                             >  1:2124354/1‑150 (MQ=255)
                        tGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCg                                                                                                                             >  2:437604/1‑150 (MQ=255)
                           aacacaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGAc                                                                                                                          <  1:2823338/145‑1 (MQ=255)
                           aaGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGAc                                                                                                                          <  2:1548444/150‑1 (MQ=255)
                           aaGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGAc                                                                                                                          <  1:1078259/150‑1 (MQ=255)
                             gAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTg                                                                                                                        <  1:1982189/150‑1 (MQ=255)
                             gAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTg                                                                                                                        >  1:249167/1‑150 (MQ=255)
                               aaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGa                                                                                                                      >  1:3494680/1‑150 (MQ=255)
                                  gcagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCa                                                                                                                   >  2:1218966/1‑150 (MQ=255)
                                  gcagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCa                                                                                                                   >  1:2419573/1‑150 (MQ=255)
                                   cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                                                                  >  1:3113905/1‑150 (MQ=255)
                                   cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                                                                  <  2:1968957/150‑1 (MQ=255)
                                    agcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAg                                                                                                                 <  1:130801/150‑1 (MQ=255)
                                       aaGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGc                                                                                                              >  1:4822096/1‑150 (MQ=255)
                                         gTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTc                                                                                                            <  2:930780/150‑1 (MQ=255)
                                            ggCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCggg                                                                                                         >  1:2990506/1‑150 (MQ=255)
                                                tGCAATTCCAACGGCAGGTCGCCGATTTCTTCAGGTTAGAATGTGCCGCCAGCGGCCGCTTCGTTTCGCCGGGTCGTGCGCTGGTGAGCACCGGTCGTGGCGCCCGGCGCATGGCCGTACGGTGCGGACTGGATGGGGTGCTCggggggc                                                                                                     >  2:720779/1‑147 (MQ=255)
                                                   aaTTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcc                                                                                                  >  1:876291/1‑150 (MQ=255)
                                                     ttCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcccc                                                                                                >  1:1966370/1‑150 (MQ=255)
                                                          aCGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGt                                                                                           <  2:2670777/150‑1 (MQ=255)
                                                           cGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGtt                                                                                          <  1:2551567/150‑1 (MQ=255)
                                                            ggCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTg                                                                                         >  1:3269981/1‑150 (MQ=255)
                                                               aGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATc                                                                                      <  1:1330020/150‑1 (MQ=255)
                                                                ggTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATcg                                                                                     <  2:261731/150‑1 (MQ=255)
                                                                ggTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATcg                                                                                     <  1:978778/150‑1 (MQ=255)
                                                                 gTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGATCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATcgc                                                                                    <  2:133791/150‑1 (MQ=255)
                                                                   cGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGa                                                                                  <  2:1431586/150‑1 (MQ=255)
                                                                   cGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGa                                                                                  <  1:2626022/150‑1 (MQ=255)
                                                                     ccGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATg                                                                                >  2:2217996/1‑150 (MQ=255)
                                                                      cGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGa                                                                               >  2:1775868/1‑150 (MQ=255)
                                                                      cGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGa                                                                               >  2:3567119/1‑150 (MQ=255)
                                                                        aTTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAAc                                                                             >  1:2579043/1‑150 (MQ=255)
                                                                          ttCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACgg                                                                           <  2:3263917/150‑1 (MQ=255)
                                                                           tcatcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgc                                                                          <  1:1311083/150‑1 (MQ=255)
                                                                            catcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgct                                                                         <  1:4149267/150‑1 (MQ=255)
                                                                            catcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgct                                                                         <  1:4149238/150‑1 (MQ=255)
                                                                             atcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgctg                                                                        <  1:4167076/150‑1 (MQ=255)
                                                                             atcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTAAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgctg                                                                        >  1:61696/1‑150 (MQ=255)
                                                                               caAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgctgct                                                                      >  2:4460693/1‑150 (MQ=255)
                                                                                aaGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTc                                                                     >  2:830070/1‑150 (MQ=255)
                                                                                    aagaaTGTGGGGCCATCGGGCGCTTCGATACGCCCGGTCTTGCGGTGGTGAGCACCGGTAATGGCGCCGTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCGCCGCTGTTGATCGCGATGAACGGGTGCTCGCTg                                                                 >  2:3388507/4‑150 (MQ=255)
                                                                                     aCAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGc                                                                <  2:3609471/150‑1 (MQ=255)
                                                                                      cAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCg                                                               <  1:3756960/150‑1 (MQ=255)
                                                                                       aaTGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGt                                                              >  2:1023593/1‑150 (MQ=255)
                                                                                        ctgtgCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCGGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGtt                                                             <  2:3868219/149‑1 (MQ=255)
                                                                                        aTGTGCCGCCATCGGCCGGTTCGATACGCCCGGTCTTGCGCTGGTGTGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACTGTTCGGCCTGGATCAAGGGCTGGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGtt                                                             <  1:653814/150‑1 (MQ=255)
                                                                                                ccATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTgg                                                     <  1:4612366/150‑1 (MQ=255)
                                                                                                  aTCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCg                                                   <  1:2124629/150‑1 (MQ=255)
                                                                                                  aTCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCg                                                   <  1:915473/150‑1 (MQ=255)
                                                                                                    cGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGt                                                 <  1:1573875/150‑1 (MQ=255)
                                                                                                      gCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGc                                               >  1:2861966/1‑150 (MQ=255)
                                                                                                         gCTTCGATACGCCCGGTCTTGCGCTGGTGTGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCTGTTGTTCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg                                            >  1:3590835/1‑150 (MQ=255)
                                                                                                         gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg                                            <  1:1712448/150‑1 (MQ=255)
                                                                                                         gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg                                            <  1:3908569/150‑1 (MQ=255)
                                                                                                         gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg                                            >  1:3485719/1‑150 (MQ=255)
                                                                                                          cTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAgg                                           <  2:1080307/150‑1 (MQ=255)
                                                                                                           ttCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGt                                          >  2:2621267/1‑150 (MQ=255)
                                                                                                           ttCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGt                                          >  2:3327884/1‑150 (MQ=255)
                                                                                                              gATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTgcg                                       >  1:587447/1‑150 (MQ=255)
                                                                                                               aTACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTgcgc                                      >  2:1621658/1‑150 (MQ=255)
                                                                                                               aTACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTgcgc                                      <  1:277202/150‑1 (MQ=255)
                                                                                                                  cGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGcc                                   <  2:3047464/150‑1 (MQ=255)
                                                                                                                   gCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCATTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGcca                                  <  2:1981135/150‑1 (MQ=255)
                                                                                                                    cccGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGccac                                 >  2:3789802/1‑150 (MQ=255)
                                                                                                                    cccGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATAAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGccac                                 <  1:150897/150‑1 (MQ=255)
                                                                                                                      cGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGccacca                               >  1:3886205/1‑150 (MQ=255)
                                                                                                                             gcgcTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCttt                        <  2:4084295/150‑1 (MQ=255)
                                                                                                                               gcTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGc                      <  1:1041606/150‑1 (MQ=255)
                                                                                                                                   gTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGt                  >  2:4554276/1‑150 (MQ=255)
                                                                                                                                   gTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGt                  <  2:938817/150‑1 (MQ=255)
                                                                                                                                     gAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGc                <  1:1163717/150‑1 (MQ=255)
                                                                                                                                     gAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGc                <  1:3096815/150‑1 (MQ=255)
                                                                                                                                       gCACCGGTAAAGGCGCCCTTCTCCTGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCg              <  2:2045493/150‑1 (MQ=255)
                                                                                                                                       gCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCg              <  1:3797736/150‑1 (MQ=255)
                                                                                                                                       gCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACCGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCg              <  1:1218966/150‑1 (MQ=255)
                                                                                                                                           cGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCGctct          <  2:140543/150‑1 (MQ=255)
                                                                                                                                             gTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCGCTCTcc        >  2:2413087/1‑150 (MQ=255)
                                                                                                                                                aaGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCGCTCTCCCcg     <  2:1797962/150‑1 (MQ=255)
                                                                                                                                                   gcgcCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCGCTCTCCCCGCgg  >  1:1879441/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCGCGCCACCAGCTCTTTGCCGGTGCCGCTCTCCCCGCGG  >  NC_002947/3975651‑3975947

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