Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,690 T→A 6.9% intergenic (‑231/+35) rluB ← / ← PP_4497 23S rRNA pseudouridylate synthase B/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,112,6900TA 6.9% 245.3 / 7.1 82intergenic (‑231/+35)rluB/PP_449723S rRNA pseudouridylate synthase B/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (44/32);  new base A (3/3);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

CAGGGGAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAGGCT  >  NC_002947/5112541‑5112838
                                                                                                                                                      |                                                                                                                                                    
caGGGGAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCtt                                                                                                                                                      <  2:2260014/150‑1 (MQ=255)
 aGGGGAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTg                                                                                                                                                     >  2:395993/1‑150 (MQ=255)
  ggggAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGc                                                                                                                                                    <  1:773234/150‑1 (MQ=255)
    ggAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGccc                                                                                                                                                  <  2:2211986/150‑1 (MQ=255)
     gAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGATCGCACAACATTCGTAGCCTTGTAATCCCCGTCAGGAGCGGGCTTGcccc                                                                                                                                                 >  1:1738034/1‑149 (MQ=255)
      aaCGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCcgc                                                                                                                                                <  1:807486/150‑1 (MQ=255)
        cGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                              <  1:57228/150‑3 (MQ=255)
          cACTGGAGGGTAGAGGGTGTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCGCTTCGGGGGCACGCCCGCGCCTACATGTTCGGGCAACCGTCGGGGCCTGGGGTACCACTGTGGGAGCGGCCTTGCCCGCGacg                                                                                                                                            >  2:670185/1‑148 (MQ=255)
           aCTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                           <  1:938707/150‑1 (MQ=255)
             tGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGTACGCCCGCTCCTACATGTCCTCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAgg                                                                                                                                         <  1:1485483/150‑1 (MQ=255)
              ggAGGGTAGAGGGTTTTGTGGGGCTGGGCATTATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        <  1:395993/150‑1 (MQ=255)
              ggAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        <  2:466703/150‑1 (MQ=255)
                       ggggTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAgg                                                                                                                               <  2:878178/149‑1 (MQ=255)
                       agGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                              >  2:885761/1‑133 (MQ=255)
                       agGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                              <  1:885761/135‑3 (MQ=255)
                                   ggCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGatcaat                                                                                                                   <  2:109557/150‑1 (MQ=255)
                                           aTGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCtt                                                                                                                                                      >  2:1669190/1‑107 (MQ=255)
                                           aTGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCtt                                                                                                                                                      <  1:1669190/107‑1 (MQ=255)
                                           aTGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                           >  2:715867/1‑118 (MQ=255)
                                           aTGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                           <  1:715867/118‑1 (MQ=255)
                                                      cccAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                              <  1:1010333/104‑3 (MQ=255)
                                                      cccAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                              >  2:1010333/1‑102 (MQ=255)
                                                      cccAGTAGGCTTTGTGCTGGCCTCTTCGCGGCCACGCCCGCTCCTCCATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGAtct                                                                                                >  2:1908997/1‑150 (MQ=255)
                                                       ccAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCCCTCCCACATGATCGCGCAACCTTCGGGGCCTGGGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGAtctc                                                                                               >  1:1630719/1‑150 (MQ=255)
                                                             ggCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTg                                                                                                                                                     <  2:458472/90‑1 (MQ=255)
                                                             ggCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTg                                                                                                                                                     >  1:458472/1‑90 (MQ=255)
                                                                       tGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggc                                                                               >  1:635291/1‑150 (MQ=255)
                                                                         gCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggctc                                                                             >  2:1248121/1‑150 (MQ=255)
                                                                           ctctTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggctcag                                                                           >  1:2214598/1‑150 (MQ=255)
                                                                            tctTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggctcagg                                                                          >  1:908493/1‑150 (MQ=255)
                                                                            tctTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggctcagg                                                                          >  1:908234/1‑150 (MQ=255)
                                                                            tctTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCAACCGGaggctaagg                                                                          <  2:635291/150‑1 (MQ=255)
                                                                              ttCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggctcaggct                                                                        >  1:1776335/1‑150 (MQ=255)
                                                                              ttCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGAACTCCTCCTCCCCCGGaggcgcaggcc                                                                        >  2:51080/1‑149 (MQ=255)
                                                                                       aCGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGATGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTcccc                                                               >  1:736499/1‑150 (MQ=255)
                                                                                         gCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        <  1:422625/75‑1 (MQ=255)
                                                                                         gCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        >  2:422625/1‑75 (MQ=255)
                                                                                           ccGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAg                                                           >  2:36864/1‑150 (MQ=255)
                                                                                             gCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGat                                                         >  2:2007251/1‑150 (MQ=255)
                                                                                              cTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGatc                                                        >  1:1185632/1‑150 (MQ=255)
                                                                                               tCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGatca                                                       >  1:1944699/1‑150 (MQ=255)
                                                                                                 cTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGatcatc                                                     >  1:1532552/1‑150 (MQ=255)
                                                                                                     aTGNTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCaaaa                                                 >  2:1500636/1‑150 (MQ=255)
                                                                                                          cgcgCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGt                                            <  2:1532552/150‑1 (MQ=255)
                                                                                                           gcgcAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGtt                                           >  2:1780235/1‑150 (MQ=255)
                                                                                                           gcgcAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGtt                                           >  1:1961529/1‑150 (MQ=255)
                                                                                                            cgcAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGttt                                          >  2:1876/1‑150 (MQ=255)
                                                                                                             gcAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGTGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCg                                              <  1:1500129/145‑1 (MQ=255)
                                                                                                             gcAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGTGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCg                                              >  2:1500129/1‑145 (MQ=255)
                                                                                                               aaCCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATctcctc                                                                                            <  2:1558563/97‑1 (MQ=255)
                                                                                                               aaCCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATctcctc                                                                                            >  1:1558563/1‑97 (MQ=255)
                                                                                                               aaCCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCa                                       >  2:2006883/1‑150 (MQ=255)
                                                                                                                  cTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGa                                                                                                                        >  1:1125181/1‑66 (MQ=255)
                                                                                                                  cTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGa                                                                                                                        <  2:1125181/66‑1 (MQ=255)
                                                                                                                   ttCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAAcc                                   <  2:1961529/150‑1 (MQ=255)
                                                                                                                   ttCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAAc                                    <  2:376663/149‑1 (MQ=255)
                                                                                                                   ttCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAAc                                    >  1:376663/1‑149 (MQ=255)
                                                                                                                      ggggCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGctcct                                                                     >  2:1476748/1‑113 (MQ=255)
                                                                                                                      ggggCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGctcct                                                                     <  1:1476748/113‑1 (MQ=255)
                                                                                                                        ggCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCtt                              >  2:536697/1‑150 (MQ=255)
                                                                                                                           cTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCa                           >  1:2303963/1‑150 (MQ=255)
                                                                                                                             tgtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                                      <  2:1627958/42‑1 (MQ=255)
                                                                                                                             tgtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                                      >  1:1627958/1‑42 (MQ=255)
                                                                                                                             tgtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAATTCGGTTTTCAAACCGTCTTCCAt                          <  2:1185632/150‑1 (MQ=255)
                                                                                                                             tgtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCgg                                             >  2:775864/1‑131 (MQ=255)
                                                                                                                             tgtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCgg                                             <  1:775864/131‑1 (MQ=255)
                                                                                                                             tgtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCAt                          <  2:1630719/150‑1 (MQ=255)
                                                                                                                              gtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATg                         <  2:2214598/150‑1 (MQ=255)
                                                                                                                              gtgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATg                         <  1:36864/150‑1 (MQ=255)
                                                                                                                                gATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCACCTCCTCCACCGGAGGCTCCGGCACCTCTACCCTCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGc                       >  1:1982864/1‑150 (MQ=39)
                                                                                                                                   aCCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCg                                              <  1:2265771/124‑1 (MQ=255)
                                                                                                                                   aCCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCg                                              >  2:2265771/1‑124 (MQ=255)
                                                                                                                                     cccTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCa                  <  1:1876/150‑1 (MQ=255)
                                                                                                                                        tgtgGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTcc                                                                 <  2:775560/100‑1 (MQ=255)
                                                                                                                                        tgtgGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTcc                                                                 >  1:775560/1‑100 (MQ=255)
                                                                                                                                         gtgGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTc              <  1:1248121/150‑1 (MQ=255)
                                                                                                                                         ctcttcgCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGATCGTCTCCCCGCACATCATCAAAAGCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTc              >  1:25343/7‑150 (MQ=38)
                                                                                                                                           agtaGCGCGCCTGCCCGCGAAAACTCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGTTTTTCAAACCGTCTTCCATGGCATCCAGCTcca            <  2:736499/147‑1 (MQ=255)
                                                                                                                                              aGCTGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGACCATCACAATCGGTTTTCAAACCGTCTTCCATCGCATCCAGCTccacca         >  2:285361/1‑150 (MQ=255)
                                                                                                                                              aGCGGGCTTGCCCGCGAAGAGGCGAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGAACTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTccaccc         >  2:1184508/1‑149 (MQ=255)
                                                                                                                                              aGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGtttt                                         <  1:1422494/118‑1 (MQ=255)
                                                                                                                                              aGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGtttt                                         >  2:1422494/1‑118 (MQ=255)
                                                                                                                                                  ggCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAg     >  1:2047840/1‑150 (MQ=255)
                                                                                                                                                  ggCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAg     >  2:752074/1‑150 (MQ=255)
                                                                                                                                                   gCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAgg    >  1:497470/1‑150 (MQ=255)
                                                                                                                                                     ttGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAGGCt  >  2:632604/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAGGGGAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAGGCT  >  NC_002947/5112541‑5112838

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