Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 798,666 C→G 6.8% intergenic (+18/+9) fklB‑I → / ← PP_0685 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947798,6660CG 6.8% 272.4 / 10.0 91intergenic (+18/+9)fklB‑I/PP_0685FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/42);  new base G (2/4);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

AGCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCAC  >  NC_002947/798518‑798805
                                                                                                                                                    |                                                                                                                                           
aGCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCt                                                                                                                                            >  1:1438268/1‑150 (MQ=255)
   tgatgaCCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGcc                                                                                                                                         <  1:2304885/150‑1 (MQ=255)
   taatcaaCGCTGGCAGCAAATGACGCCTGTACGTGCCAATAGCGCTGGCTTACAGCGACAAAGACGTTGGCAGCATCCCGCCGCGAATCGTGCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGCAATGCCCTGGCCGCCGGACCTGcc                                                                                                                                         <  2:850889/145‑1 (MQ=255)
      tgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGc                                                                                                                                      <  1:342358/150‑1 (MQ=255)
        aaCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCt                                                                                                                                    <  2:1467989/150‑1 (MQ=255)
        aaCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCt                                                                                                                                    >  1:496567/1‑150 (MQ=255)
         aCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTc                                                                                                                                   >  2:1491028/1‑150 (MQ=255)
           gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                 <  2:1457498/150‑3 (MQ=255)
           gCTGGCAGCAAATGGCGCCCGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                 <  2:608663/150‑3 (MQ=255)
            cTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTCaag                                                                                                                                <  1:82062/150‑2 (MQ=255)
              ggCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTc                                                                                                                              <  2:686616/150‑1 (MQ=255)
              ggCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTc                                                                                                                              <  2:1685170/150‑1 (MQ=255)
               gcagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTcc                                                                                                                             >  1:1484668/1‑150 (MQ=255)
               gcagcaAATGGCGCCTGTACGTGCCAAGAGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTcc                                                                                                                             >  1:38614/1‑150 (MQ=255)
                cagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                            <  2:1484668/150‑1 (MQ=255)
                cagcaAATGGCGCCTGGACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCGCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                            <  2:247394/150‑1 (MQ=255)
                  gcaAATGGCGCCTGTAAGTGCCAAGCGAGCTGGCTTACGGCGCCCGGAATGGTGGCAGGATCCCGCCGCACAGCGTTCTGGTGTTCGACGTTTATTTGCTGGACTTTCTGTAATGCCCGTGTGGGCGGCCTTTCTGGCCGCTATTCTAGt                                                                                                                          >  2:327332/1‑150 (MQ=255)
                   caAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTc                                                                                                                         >  2:765857/1‑150 (MQ=255)
                      aTGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGt                                                                                                                      >  2:379036/1‑150 (MQ=255)
                        ggCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGGCGAGGTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTcc                                                                                                                    <  2:1218955/150‑1 (MQ=255)
                         gcgcCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGCGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTccc                                                                                                                   >  1:2076482/1‑150 (MQ=255)
                           gcCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTccccc                                                                                                                 <  1:379036/150‑1 (MQ=255)
                                    tGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCa                                                                                                        >  2:1286943/1‑150 (MQ=255)
                                    tGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCa                                                                                                        >  1:758249/1‑150 (MQ=255)
                                     gCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCaa                                                                                                       >  1:377178/1‑150 (MQ=255)
                                      ccAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAAc                                                                                                      >  1:1886895/1‑150 (MQ=255)
                                        cagagAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGc                                                                                                    >  2:1440805/5‑150 (MQ=255)
                                        aaGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGc                                                                                                    >  2:1633286/1‑150 (MQ=255)
                                               cTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCa                                                                                             >  1:528883/1‑150 (MQ=255)
                                                  gCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAg                                                                                          <  2:528883/150‑1 (MQ=255)
                                                    ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGca                                                                                        >  2:224544/1‑150 (MQ=255)
                                                    ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGca                                                                                        <  2:92968/150‑1 (MQ=255)
                                                        ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaac                                                                                    >  2:883180/1‑150 (MQ=255)
                                                          cgcCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCATCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaacag                                                                                  <  2:472176/150‑1 (MQ=255)
                                                            cccAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaacagaa                                                                                <  1:2316768/150‑1 (MQ=255)
                                                             ccAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaa                                                                               >  1:2325169/1‑150 (MQ=255)
                                                              cAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCTACGTTCTGTAATGCCCTGGGGGCCGGCCCTGCCGTCCTCAATTCCAGTCCGTTCTCCCAGGGCGCAACGCGCGCGCATAGCAGaacagaaac                                                                              >  1:1653657/1‑150 (MQ=255)
                                                              cAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaac                                                                              >  1:53284/1‑150 (MQ=255)
                                                              cAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaac                                                                              >  2:1126262/1‑150 (MQ=255)
                                                              cAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGAAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACCCCCGCGCATAGCAGaacagaaac                                                                              >  2:155689/1‑150 (MQ=255)
                                                               aaGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAAc                                                                                                      <  1:487613/125‑1 (MQ=255)
                                                               aaGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAAc                                                                                                      >  2:487613/1‑125 (MQ=255)
                                                                  gCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTccagg                                                                                                                               >  1:847855/1‑93 (MQ=38)
                                                                  gCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTccagg                                                                                                                               <  2:847855/97‑5 (MQ=38)
                                                                    gTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTc                                                                        <  2:2325169/150‑1 (MQ=255)
                                                                           gcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCACCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCa                                                                 <  2:2314506/150‑1 (MQ=255)
                                                                               cccGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTc                                                             <  1:765857/150‑1 (MQ=255)
                                                                                 cgccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCt                                                           >  2:1344911/1‑150 (MQ=255)
                                                                                    cgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGa                                                        <  1:1491028/150‑1 (MQ=255)
                                                                                    cgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGa                                                        <  2:377178/150‑1 (MQ=255)
                                                                                      cacaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGc                                                      <  1:1633286/150‑1 (MQ=255)
                                                                                       acaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCa                                                     >  2:1851580/1‑150 (MQ=255)
                                                                                        caGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                            >  1:1870213/1‑78 (MQ=255)
                                                                                        caGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                            <  2:1870213/78‑1 (MQ=255)
                                                                                         aGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCAcc                                                   <  2:53284/150‑1 (MQ=255)
                                                                                             ttCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCg                                               <  1:1665068/150‑1 (MQ=255)
                                                                                              tCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCAAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCgg                                              >  2:1687642/1‑150 (MQ=255)
                                                                                                 ggTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTc                                                                                                                     <  1:1780587/76‑1 (MQ=255)
                                                                                                 ggTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTc                                                                                                                     >  2:1780587/1‑76 (MQ=255)
                                                                                                   tgtTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCAc                                         <  1:883180/150‑1 (MQ=255)
                                                                                                    gtTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACt                                        >  1:956049/1‑150 (MQ=255)
                                                                                                    gtTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACt                                        >  1:1079418/1‑150 (MQ=255)
                                                                                                     ttCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGGCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTc                                       >  1:1120465/1‑150 (MQ=255)
                                                                                                     ttCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTc                                       <  2:1886895/150‑1 (MQ=255)
                                                                                                     cagaaCGTCGAGCAGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTcc                                                            >  2:808803/5‑129 (MQ=255)
                                                                                                     cagaaCGTCGAGCAGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTcc                                                            <  1:808803/125‑1 (MQ=255)
                                                                                                      agaaCGTCGAGCAGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCg                                      >  2:2120933/4‑150 (MQ=255)
                                                                                                       gaaCGTCGAGCAGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGa                                     >  1:417351/3‑150 (MQ=255)
                                                                                                       cGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGa                                     >  2:962147/1‑150 (MQ=255)
                                                                                                        gACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                    <  1:1126262/150‑1 (MQ=255)
                                                                                                            gcGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCa                                <  1:1361158/149‑1 (MQ=255)
                                                                                                             cGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAg                               <  2:2034717/150‑1 (MQ=255)
                                                                                                               aGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGcc                             <  2:2050786/150‑1 (MQ=255)
                                                                                                                    ctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTc                        >  2:463114/1‑150 (MQ=255)
                                                                                                                           gtCTGTAATGCCCTGGAGGACGGCGATGACCGACTCAAATACAGCCTATCCCACCAAGGCGCAACGCCCGCGCCTACCAAAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCa                 <  2:1653657/149‑1 (MQ=255)
                                                                                                                            tCTGTACTGCCCTGGCGCCCGCCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                <  2:1120465/150‑1 (MQ=255)
                                                                                                                            tCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                <  2:2076482/150‑1 (MQ=255)
                                                                                                                            tCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                <  1:1344911/150‑1 (MQ=255)
                                                                                                                            tCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                <  1:327332/150‑1 (MQ=255)
                                                                                                                              tGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAg                                                       <  1:958508/109‑1 (MQ=255)
                                                                                                                              tGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAg                                                       >  2:958508/1‑109 (MQ=255)
                                                                                                                                      ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCa                                >  1:85843/1‑124 (MQ=255)
                                                                                                                                      ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCa                                <  2:85843/124‑1 (MQ=255)
                                                                                                                                      ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGc                          <  1:6145/130‑1 (MQ=255)
                                                                                                                                      ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGc                          >  2:6145/1‑130 (MQ=255)
                                                                                                                                       cTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGt     <  1:1440805/150‑1 (MQ=255)
                                                                                                                                        tGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTc    <  2:417351/150‑1 (MQ=255)
                                                                                                                                         ggAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAAAAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCa   >  2:555341/1‑150 (MQ=255)
                                                                                                                                         ggAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTCTCCCTCTACGGCGCTACGCCCGCGCATAGCAGAACAGACCCAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCa   <  1:746252/150‑1 (MQ=255)
                                                                                                                                          gAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCACCGGTTCAGTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCAc  >  2:362663/1‑150 (MQ=255)
                                                                                                                                          gAGGCCGGCAGGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATa                                                                                           >  2:1492269/1‑61 (MQ=37)
                                                                                                                                                    |                                                                                                                                           
AGCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCAC  >  NC_002947/798518‑798805

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