Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,726,259 T→A 14.6% intergenic (‑183/+111) PP_5024 ← / ← opgH amino acid ABC transporter substrate‑binding protein/glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,726,2590TA14.6% 218.3 / 4.4 96intergenic (‑183/+111)PP_5024/opgHamino acid ABC transporter substrate‑binding protein/glucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base T (26/50);  new base A (12/2);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 6.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

TTATTGTGTTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAGG  >  NC_002947/5726110‑5726408
                                                                                                                                                      |                                                                                                                                                      
ttATTGTGTTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCt                                                                                                                                                        <  2:4992165/150‑1 (MQ=255)
  aTTGTGTTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTAc                                                                                                                                                      <  1:535875/150‑1 (MQ=255)
     gtgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAg                                                                                                                                                   <  1:4653518/150‑1 (MQ=255)
     gtgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTGAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAg                                                                                                                                                   <  1:2963729/150‑1 (MQ=255)
      tgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAgg                                                                                                                                                  >  1:3202956/1‑150 (MQ=255)
        ttGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGc                                                                                                                                                <  2:2672374/150‑1 (MQ=255)
        ttGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGc                                                                                                                                                <  2:2672242/150‑1 (MQ=255)
          ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGCGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTGAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                              <  1:452334/150‑1 (MQ=255)
          ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                              <  2:91134/150‑1 (MQ=255)
           gCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAtg                                                                                                                                             <  2:2161221/150‑1 (MQ=255)
           gCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAtg                                                                                                                                             <  1:2086362/150‑1 (MQ=255)
           gCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCCCGCCAAAAAGGGCaag                                                                                                                                             >  2:243890/1‑148 (MQ=255)
           gCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTCGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAtg                                                                                                                                             <  1:671595/150‑1 (MQ=255)
            cACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAtgt                                                                                                                                            <  2:2293500/150‑1 (MQ=255)
            cACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAtgt                                                                                                                                            <  1:5094718/150‑1 (MQ=255)
             aCGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTACCCGCACCAACAAGGGCCtgtt                                                                                                                                           >  1:1239352/1‑149 (MQ=255)
             aCGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAgggatacac                                                                                                                                           <  2:2708865/150‑7 (MQ=255)
             aCGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                              >  1:443853/1‑147 (MQ=255)
             aCGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                              <  2:443853/147‑1 (MQ=255)
             aCGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAAACCCCACCAACAAGGGCAggtg                                                                                                                                           >  2:1321815/1‑150 (MQ=255)
               ggCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTa                                                                                                                                         <  2:3202956/150‑1 (MQ=255)
                ggcGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTAt                                                                                                                                        <  2:3174269/148‑1 (MQ=255)
                gCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCGCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTAt                                                                                                                                        <  2:1733281/150‑1 (MQ=255)
                gCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCGAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTAt                                                                                                                                        <  2:1733167/150‑1 (MQ=255)
                   gTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATccc                                                                                                                                    <  2:1062106/150‑1 (MQ=255)
                     aCTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTg                                                                                                                                  <  2:2857202/150‑1 (MQ=255)
                     aCTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTg                                                                                                                                  <  1:2272415/150‑1 (MQ=255)
                      cTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTgg                                                                                                                                 <  2:4631436/150‑1 (MQ=255)
                        ggCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTa                                                                                                                               <  2:687989/150‑1 (MQ=255)
                         gCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGATAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAg                                                                                                                              <  1:3510442/150‑1 (MQ=255)
                          ccGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAgg                                                                                                                             <  1:4871180/150‑1 (MQ=255)
                           cGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGa                                                                                                                            <  2:4313177/150‑1 (MQ=255)
                           cGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCAACAAGGGCCTGTGTATCCCggggggga                                                                                                                            >  1:3305121/1‑145 (MQ=255)
                           cGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAAACCCCTACAAAAAGGGCAAGTGTATCCCAGGTAgtg                                                                                                                            >  1:3305166/1‑148 (MQ=255)
                           cGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAAACCCCGCCAACAAGGGCATGTGTATACCGGGTAGGa                                                                                                                            >  1:3305088/1‑150 (MQ=255)
                             aaGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGc                                                                                                                          <  2:1422353/150‑1 (MQ=255)
                               gCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCgg                                                                                                                        <  2:1141541/150‑1 (MQ=255)
                               gCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGCTTCGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACACGGGCATGTGTATCCCTGGTAGGAGCgg                                                                                                                        <  1:1459437/150‑1 (MQ=255)
                                cgtcgtATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCggg                                                                                                                       <  1:4563772/150‑1 (MQ=255)
                                         ggAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAATGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGGAGGAGCGGGTTCACCcgc                                                                                                              >  2:2981538/1‑150 (MQ=255)
                                               ggATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGGACCCCCGCCGACCAGGGCCGGTGTATACCGGGGGGGGGCGGGTTCACCcccggcgag                                                                                                        >  2:57119/1‑145 (MQ=255)
                                               ggATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGAACCCCCGCCAACAACGGCAGGTGTATCCCGGGGAGGAGCGGGTGCACCcccgccgcg                                                                                                        >  2:4791119/1‑145 (MQ=255)
                                                  tGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTAACCCCGCCAACAAAGGCATGTGTATAACAGGTAGGAGCGGGTTGAACCCCGACGAGGcc                                                                                                     >  1:3213178/1‑150 (MQ=255)
                                                   gCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCTGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGccc                                                                                                    <  2:300291/150‑1 (MQ=255)
                                                    cTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCGCCTACAAGGGCATGTGTATCCCGGGTAGGAGCGGGGTCACCCGCGAAGAGGCCCg                                                                                                   >  2:1185341/1‑150 (MQ=255)
                                                     taccGCTTGGGGGCCGGTGCACTTAGCGCCGCTGGGGGATTTTGTGGCCCAATTGGTGGGTTTTGGCCAAAAATGTTCGGGGTGTCACCCGGCACCCACCCGGGCTTGTTTATCCGTGTTTGGCGGGGGTACACCCGCGAAGAGGCCCgg                                                                                                  <  2:1730554/147‑1 (MQ=255)
                                                           aTGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTACCCGCACCTAAAACGGGAAGTGTATACCGGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATgc                                                                                            >  1:3779514/1‑150 (MQ=255)
                                                           aTGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGGACCCCCACCAACAACGGCATGTGTATCCCGGGGAGGAGCGGGTTGACCCGCGCAGAGGCCCGGGCATgc                                                                                            >  1:3779513/1‑150 (MQ=255)
                                                                   ccGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTAACCGCTCCTACAAGGGCATGTGTATCCCTGGGAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTa                                                                                    >  2:419851/1‑150 (MQ=255)
                                                                      tGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGAAACCGCGCCAACAAAGGCATGTGTATCCCGGGGGGGAGCCGGTTCACCCGCGAAGTGGCCCGGGCATGCGcccccgtggg                                                                                 >  1:366434/1‑144 (MQ=255)
                                                                        gACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCg                                                                               >  2:4618184/1‑150 (MQ=255)
                                                                         ccctcgcgccgccGGGGGACTTTGTGGGCAAAAATGCGGGGTTTAGGCAAAAAAGCTTCGGGGGTGAACCCGCTCCAACACGGGATTGTGTAGCCCGGGTCGGAGCGGGTACACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCgg                                                                              <  2:618191/143‑1 (MQ=255)
                                                                         cccTAGCGCCGCAGGGGGACTTTTTTGGCCAAATTGCGGGGTTTTGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACATGGGCATGTTTAGCCCTGGTAGGCGGGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCgg                                                                              <  2:2653479/149‑1 (MQ=38)
                                                                           ctcgcgCCGCAGGGGGACTTTGTTGGCCAAAAGGCGGGGGTTCGGCCAAAAAGCTTCGCGGGTGACCCCGCTCCTACCCGGGCATGTGTAGCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGt                                                                            <  2:4826954/145‑1 (MQ=38)
                                                                            tAGCGCCGCCGGGTGACTTTTTTGGCAAAAATGAGGGGGTTAGGCAAAAACGCTTCGGGGGTGACCCCGCTCCTACAAGGGCATGTGTATCCGTGTTTGGAGCGGTTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGtt                                                                           <  2:1819349/150‑1 (MQ=255)
                                                                             agagCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGGACCCCCGCCTAAAACGGCCTGGGTATCCCTGGTAGGAGCGGGTTCCCCCGCGAAGGGGCCCGGTCCTGCGCCCCGGAGGGCGGGTTc                                                                          >  1:4626213/1‑150 (MQ=255)
                                                                                   gccgggggACGTGTTTGGCCAACCGGCGGGGGTTCGGCAAAAACGCTGCGGGGGGGACCCCGCGCCCCCACGGGCTGGTGGAGCCGTGGGTGGCGGGGGTCCCCCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCggg                                                                    <  2:610198/147‑1 (MQ=255)
                                                                                    ctgggggACTTTGTTGGCAAAAATGCGGGGGTTCGGCCAAAAAGCCGCGCGGGTGAACCCGCTCCTACACGGGCATGTGAATCCCTGGTCGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCgggg                                                                   <  2:3496540/148‑1 (MQ=38)
                                                                                       ggggACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCGCCGAAAAGGGCATGTGTTTCCCGGGTAGGGGCGGGTTGACCCGCGCAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGGTCCTCGGGgcgc                                                                >  1:3369356/1‑150 (MQ=255)
                                                                                       ggggACTTGTTTGCCCAAATTGTGGGGGTTAGTCAAAAACGCTTGGGGGGTTACCCCGCGCCTACCCGGGCTTGTGTATCCCTGGTTGGAGCGGTTACACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGgcgg                                                                <  2:185474/150‑2 (MQ=255)
                                                                                         ggACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAATGCTTCGCGGGTGTACCCGCACCAAAAAGGGCAGGTGTATACCAGGGAGGAGCGGGGTCACCCGCGAAGAGGCCCCGGCAATCGCCCCGGTGGGCGGGGTAATCGGGGCGCCg                                                              >  1:4650168/1‑150 (MQ=255)
                                                                                             tttGTTGGCCAAAATGCGGGGGTTCGCCCAAAAAGCTTCGGGGGGCAACCCGCCCCAACACGGGCATGTGTACCCCGGGTAGGCGCGGGTACACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAAcc                                                          <  2:674609/150‑1 (MQ=255)
                                                                                             tttGTGGCCCAACATGCGGGGGTTCGGCCAAAACGCGTCGGGGGTGAACCCGCCCCCACACGGGCTTGTGTATCCCGGGTCGGCGGGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAAcc                                                          <  2:541465/150‑1 (MQ=255)
                                                                                                  tgGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGAACCCGCTCCTACAAGGGAATGTGTAACCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACcgccgc                                                     <  2:2698960/150‑1 (MQ=255)
                                                                                                  tgGCCAAAATGAGGGGGACCGCCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACcgccgc                                                     <  1:4188509/150‑1 (MQ=255)
                                                                                                   ggCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCCCGCCAACAAGGGCAGGTGTATCCCCGGGAGGAGCCGGGTCACCCCCCCAGAGGCCCCGGCCTGCGCCCCGGAGGGCGGGGTCATCGGGGCGCCGCACccccgcc                                                    >  2:1063199/1‑149 (MQ=255)
                                                                                                     caaaaaTGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACACGGGCTTGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGaa                                                  <  2:3173407/148‑1 (MQ=255)
                                                                                                         aaggcgggggTTCGGCCAAAACGCTGCGCGGGTGAACCCGCGCCTACACGGGCCTGTGTATCCCTGGTCGGAGCGGGTTCCCCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGa                                              <  2:3104676/145‑1 (MQ=255)
                                                                                                           tgtgggggTTCGGCAAAAACGCTTGGGGGTTCACCCCGCTCCTAAACGGGATTGTGTATCCCTGTTTGGTGGGGTTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGAcc                                            <  1:1185341/147‑1 (MQ=255)
                                                                                                            gtgggggTTAGGCCAAAAAGCTTCGGGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTACACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGAccc                                           <  2:1423512/148‑1 (MQ=255)
                                                                                                              gggggTTCGGCAAAAACGCTGCGGGGGTGAACCCGCTCCTACCCGGGCTTGTGTAGCCGTGGTAGGAGGGGGTCCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGa                                         <  2:1008778/150‑1 (MQ=255)
                                                                                                              gggggTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCCCCTACAAGGGATTGTGTATCCCTGGTAGTAGGGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGa                                         <  2:552139/150‑1 (MQ=39)
                                                                                                                        ccAAAAAGCTTCGCGGGTGAACCCCCGCCGACAAGGGCAGGTGTATACCGGGTGGGCGCGGGGTAACCCCCGAAGAGGCCCCCTCAATCGCCCCCTAGGGCGGGTTAAACGGGGCGCCGCACCCCCCCGAACACACCCCAACAAGCGCcg                               >  2:2641311/1‑150 (MQ=255)
                                                                                                                              aGCTTCGCGGGTGAATCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAgg                         <  2:4462968/150‑1 (MQ=255)
                                                                                                                                    gcgGGTGAACCCGCTCCTACCAGGGATACACATGCCCGTGGAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGt                   >  1:4561234/1‑150 (MQ=25)
                                                                                                                                          tGAACCCGCTCCTACCAGGGATACACATGCCCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGcc              >  1:616938/1‑150 (MQ=25)
                                                                                                                                          tGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGcc              >  2:1573738/1‑150 (MQ=255)
                                                                                                                                            aaCCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCACGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTg            >  1:3241207/1‑150 (MQ=255)
                                                                                                                                            aaCCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCACGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTg            >  1:2692127/1‑150 (MQ=255)
                                                                                                                                             aCCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGc           >  2:4915682/1‑150 (MQ=255)
                                                                                                                                             aCCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGc           >  1:4665880/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGt         >  1:2894866/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGt         >  1:3025132/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGt         >  1:2953028/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCTACAAGGGCATGTGTATCACTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGCCTCGGTTGGCCTGCGt         >  1:4065446/1‑150 (MQ=255)
                                                                                                                                                 gccccTACATGGGCTTTTTTTTCCGTGTTTGGTGGGGGTTCCCCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt       <  1:1917526/147‑1 (MQ=255)
                                                                                                                                                 gCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt       >  2:3146906/1‑150 (MQ=255)
                                                                                                                                                 gCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGAGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt       >  2:1008142/1‑150 (MQ=255)
                                                                                                                                                 gCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGAGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt       >  2:2579858/1‑150 (MQ=255)
                                                                                                                                                  ccccTACCGGGGCTTGTTAACCCCTGTTTGGCGGGGGTACCCCCGCGACGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGtttt      <  1:1051930/148‑1 (MQ=255)
                                                                                                                                                    ccTACACGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCa    <  2:3284208/150‑1 (MQ=255)
                                                                                                                                                    ccTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCa    >  1:4282250/1‑150 (MQ=255)
                                                                                                                                                     cTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCTTAGGGCGGGTTCATCGGGGCGCCGAACCCCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAg   >  1:3646079/1‑150 (MQ=255)
                                                                                                                                                     cTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCggg                                                                             >  1:3839269/1‑76 (MQ=255)
                                                                                                                                                     cTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCggg                                                                             <  2:3839269/76‑1 (MQ=255)
                                                                                                                                                     cTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCACGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTATTCAg   >  2:2008127/1‑150 (MQ=255)
                                                                                                                                                      tACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAgg  <  2:4785146/150‑1 (MQ=255)
                                                                                                                                                      tACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAgg  >  2:4784671/1‑150 (MQ=255)
                                                                                                                                                      tACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAgg  >  1:4276627/1‑150 (MQ=255)
                                                                                                                                                      tACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAgg  >  2:3537182/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TTATTGTGTTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAGG  >  NC_002947/5726110‑5726408

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