Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,726,262 A→C 8.6% intergenic (‑186/+108) 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,2620AC 8.6% 402.7 / 16.9 144intergenic (‑186/+108)PP_5024/opgHamino acid ABC transporter substrate‑binding protein/glucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base A (69/61);  new base C (6/8);  total (75/69)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

TGTGTTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAGGGT  >  NC_002947/5726114‑5726410
                                                                                                                                                    |                                                                                                                                                     
  tgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACCAgg                                                                                                                                                    <  2:4029317/150‑1 (MQ=255)
  tgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAgg                                                                                                                                                    <  2:1695406/150‑1 (MQ=255)
  tgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAgg                                                                                                                                                    <  2:1875584/150‑1 (MQ=255)
  tgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAgg                                                                                                                                                    <  2:4069916/150‑1 (MQ=255)
  tgtTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGGACCCGCTCCAAAAAgg                                                                                                                                                    >  1:2067977/1‑150 (MQ=255)
  tgtTGGCACGGCCGTACGGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACACgg                                                                                                                                                    <  1:1126969/150‑1 (MQ=255)
    ttGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGGAACCCCACCGACAAGGGc                                                                                                                                                  >  2:383633/1‑150 (MQ=255)
    ttGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGAACCCCCACCAAAAAGGGc                                                                                                                                                  >  2:2640366/1‑150 (MQ=255)
      ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTTACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCACAATGAGGGGGTTAGGCCAAAATGCTTCGCGGGTGAACCCGGTCCTACAAGGGCAt                                                                                                                                                <  2:2012431/150‑1 (MQ=255)
      ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                                <  1:1068761/150‑1 (MQ=255)
      ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                                <  1:723643/150‑1 (MQ=255)
      ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                                <  1:436357/150‑1 (MQ=255)
      ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                                <  2:3517077/150‑1 (MQ=255)
      ggCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAt                                                                                                                                                <  2:3612069/150‑1 (MQ=255)
       gCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCGACAAGGGCagg                                                                                                                                               >  2:4028711/1‑148 (MQ=255)
       gCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCACCTACCAGGGCAtg                                                                                                                                               >  2:2188382/1‑150 (MQ=255)
       gCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGCTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCAtg                                                                                                                                               <  1:1283069/150‑1 (MQ=255)
           ggCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAAACCCCTCCAAAAAGGGCAAGTGTa                                                                                                                                           >  1:1822736/1‑150 (MQ=255)
             ccGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATc                                                                                                                                        <  2:4304202/150‑1 (MQ=255)
               gTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATccc                                                                                                                                      <  2:1317801/150‑1 (MQ=255)
                tACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCt                                                                                                                                     >  2:674716/1‑150 (MQ=255)
                tACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCt                                                                                                                                     <  2:1862403/150‑1 (MQ=255)
                tACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGGACCCCCGCCGACAACGGCCTGTGTATaccg                                                                                                                                     >  1:433412/1‑149 (MQ=255)
                 tcTGGCCGAAGCTTCGTATCGGTGAGGGATGCTAACTCATGGCGGCACGTGGACATAGAGCCGCAGTGGTACTTTGTTGGCCATAATGAGGGGTTTAGGCCAAAAAGCTTCGCGTGTGAACCCTCTCCTACAAGGGCATGTGTATCCCtt                                                                                                                                    <  2:421043/149‑2 (MQ=255)
                 aCTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTg                                                                                                                                    <  2:477564/150‑1 (MQ=255)
                  cTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTgg                                                                                                                                   <  1:3250988/150‑1 (MQ=255)
                    ggCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTa                                                                                                                                 <  2:2421352/150‑1 (MQ=255)
                    ggCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTa                                                                                                                                 <  2:294838/150‑1 (MQ=255)
                    ggCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTa                                                                                                                                 <  2:1256206/150‑1 (MQ=255)
                    ggCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTa                                                                                                                                 <  2:2104769/150‑1 (MQ=255)
                    ggCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAg                                                                                                                                <  2:2992692/150‑1 (MQ=255)
                     gCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAg                                                                                                                                <  2:1850945/150‑1 (MQ=255)
                      ccGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAgg                                                                                                                               <  2:735566/150‑1 (MQ=255)
                       cGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGa                                                                                                                              >  2:3704801/1‑150 (MQ=255)
                       cGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCGGGTAggt                                                                                                                              >  2:2565705/1‑149 (MQ=255)
                         caGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGc                                                                                                                            <  2:3061279/149‑1 (MQ=255)
                         aaGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTACCCGCTCCAACAAGGGCATGTGTATCCCGGGGAGGAGc                                                                                                                            >  1:4084484/1‑150 (MQ=255)
                         aaGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGc                                                                                                                            <  1:4428531/150‑1 (MQ=255)
                           gCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACCAGGGCATGTGTATCCCTGGTAGGAGCgg                                                                                                                          <  2:1208224/150‑1 (MQ=255)
                           gCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCgg                                                                                                                          <  1:595143/150‑1 (MQ=255)
                            cgtcgtATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCCCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCggg                                                                                                                         <  1:1841406/150‑1 (MQ=255)
                            cgtcgtATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGCTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCggg                                                                                                                         <  1:2238513/150‑1 (MQ=255)
                             gtcgtATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAACAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGt                                                                                                                        <  2:2646030/150‑1 (MQ=255)
                              tcgtATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAAACCGCTCCTACAAGGGCAAGTGTATACCAGGTAGGAGCGGGtt                                                                                                                       >  1:4419152/1‑150 (MQ=255)
                               cgtATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTc                                                                                                                      <  2:3458646/150‑1 (MQ=255)
                                    cGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCcg                                                                                                                 <  1:3565656/150‑1 (MQ=255)
                                     ggAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCCCGCCAACAACGGCATGTGTATCCCGGGGTGGAGCGGGTTCAccccc                                                                                                                >  1:1462660/1‑148 (MQ=255)
                                      gagaGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCGACAAGGGCATGTGTGTCCCTGGGAGGAGCGGGGTCACCcgcg                                                                                                               >  1:1096817/1‑150 (MQ=255)
                                       agagGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCCGGGGGACTTTGTTGGCCAAACTGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACCAGGGATACACATGCCCTTGTAGGAGCGGGTGAACCcgctc                                                                                                              <  2:1617014/150‑3 (MQ=21)
                                       agagGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACCAGGGATACACATGCCCTTGTAGGAGCGGGTGAACCcgctc                                                                                                              <  2:2222527/150‑3 (MQ=21)
                                          gggATGCTAACGCATGGCGGCCCGTGGACCTCGCGCCGCCGGGGGCCTTTGTTGGCCAACCGGAGGGGGTTCGGCCAAAAAGCTGCGCGGGTGAACCCGCTCCTACACGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAaga                                                                                                           <  2:4452588/150‑1 (MQ=38)
                                                 tAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCGCCAAAAAGGGCATGTGGGTACCGGGTAGGGGCGGGGTCACCCGCGAAGAGGCCCgg                                                                                                    >  2:1510729/1‑150 (MQ=255)
                                                 tAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGTACCCCCGCCAAAAAAGGGCGGGGTATGACCGGTAGGGGCGGGGTGGCCCGCGCAGAGGCCCgc                                                                                                    >  1:473632/1‑149 (MQ=255)
                                                 tAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGTAACCCCCCCAACAACGGCCAGGGTATCCCAGGGAGGAGCGGGGTAACCCCCGCACAGGccccc                                                                                                    >  1:1574420/1‑148 (MQ=255)
                                                            ggcCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGGTCCCCCGCCAACAACGGCCTGTGTTTACCCGGTAGGAGCGGGGTCACCCCCCCAAATGCCCCGGCAGGCGcccc                                                                                         >  1:3469862/1‑150 (MQ=255)
                                                            ggcCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGGGGACCCCCGCCGACACGGGGCGGGGGATCCCGGGGTGGGGCGGGGGGACCCGCGGAGAGGCCCCGGCCGGCGcccc                                                                                         >  2:1501360/1‑150 (MQ=255)
                                                                 gTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTAAAAGGGCCTGTGTATCCCTGGTAGGAGCGGGGTCACCCGCGCAGAGGCCCCGTCAAGCGCCCCGTAgg                                                                                    >  2:2803523/1‑150 (MQ=39)
                                                                     cccTAGAGCCGCAGGGGGACTTGGTGGGCCAAAATGTGGGGGTTCGGCCAAAAAGCTTCGCGGTTGACCCCGCTCCTACACGGGCATGTGTAGCCCTGGTCGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCgg                                                                                <  2:3905217/149‑1 (MQ=38)
                                                                       cTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCCCGCCGACAAGGGGATGTGGTTCCCGGGGCGGAGCGGGGTCACCACCGGAGAGGCCCGGGCATTCGCCCCGGAGGGCgggg                                                                              >  2:208664/1‑149 (MQ=255)
                                                                         agagCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAATGCTTCGCGGGTGAAGCCGCTGCTTCAAGGGGATGTGTATGGCTGGTTGGTGCGGGTTGACCCGCGAAGAGGGCCGGTCATGCGGCCCGTAGGGCGGGttg                                                                            >  2:2932043/1‑149 (MQ=255)
                                                                            cccgccgGGGGACTTTGTGGGCCAAACTGCGGGGGTTCGGCCAAAAAGCTGCGGGGGTGACCCCGCTCCTACACGGGCTTGTGTATCCCTGGTAGGAGCGGGTACACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATc                                                                         <  2:1291335/149‑1 (MQ=255)
                                                                             ccGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCGCCAACAAGGGCCTGGGTAACCCTGGTAGGAGCGGGTTCACCCGCCAAGAGGCCCGGTCATGCGCCCCGTGGGGCGGGGTCATCg                                                                        >  2:4307907/1‑150 (MQ=255)
                                                                              cGCTGGGGGAGTTTGTTGGCCAAAATGAGGGGGTGCGGCCAAAAAGCTTCGCGGGTGACCCGGCGCCTACACGGGCTTGTGTACCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCgg                                                                       <  1:4128967/150‑1 (MQ=38)
                                                                                cAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTAACCGCGCCTACAAGGGGCTGTGTATCCCGGGGGGGGGCGGGTTCACCCCCGAAGAGGCCCCGGCATTCGCCCCCGAGGGCGGGGTCATCgggg                                                                     >  2:423071/1‑150 (MQ=255)
                                                                                       aCTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCAACAAGGGCCTGTGTATCCCTGGTAGGGGCGGGTTCACCCGCGCAGAGGCCCGGTCCTGCGCCCCGGAGGGCGGGTTCATCGGGGCGCCGaa                                                              >  2:3858980/1‑150 (MQ=255)
                                                                                          ttgttgGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCGGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAAccg                                                           >  2:1452308/1‑150 (MQ=255)
                                                                                           tgttgGCAAAAATGAGGGGGTTAGCCCAAAAAGCTTGGGGGGTTACCCCGCCCCTACCCGGGCTTGTGTACCCCTGGTCGGCGGGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAAccgc                                                          <  2:409641/150‑1 (MQ=255)
                                                                                              tgGCCAAATTGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACcgccgc                                                       <  2:3600462/150‑1 (MQ=255)
                                                                                                 ccAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGTAACCCCACCGACAACGGCCTGTGTATCCCGGGTAGGAGCCGGGTCACCCCCCCAGAGGCCCCGGCAAGCGCCCCCGAGGGCGGGGTCAACGGGGGGCCCAACCCCcccgcc                                                    >  2:1171241/1‑148 (MQ=255)
                                                                                                     aaTGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACACGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGa                                                <  2:3094318/150‑1 (MQ=255)
                                                                                                          gggggTTAGGCCAAAAAGCTTCGCGGGTGAACCCCCTCCAACAAGGGCAAGTGTATACCTGGGAGGAGCGGGTTCAACCGCGAAGAGGCCCGGTCAAGCGCCCCCTAGGGCGGGTTCATAGGGGCGCCGAACCGCCGCGAACAGAACCGa                                           >  2:3227903/1‑150 (MQ=255)
                                                                                                           ggggTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCACCAACAAAGGCCGGTGTATCCCGGGGTGGAGCGGGTTCACCCCCGACGAGGCCCCGGGATGCGCCCCCGAGGGCGGGTTCAAAGGGGCGCCGACCCCCCGCGCACACAccccaa                                          >  2:2270578/1‑147 (MQ=255)
                                                                                                              gTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCaa                                       <  2:1676490/150‑1 (MQ=255)
                                                                                                              gTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCaa                                       <  2:2591111/150‑1 (MQ=255)
                                                                                                                                 cgGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGtt                    <  2:3300881/150‑1 (MQ=255)
                                                                                                                                  gggTGAACCCGCTCCTACCAGGGATACACATGCCCTTATAGGAGCGGGTTCACCCGCGAAGCGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTg                   >  2:3480775/1‑150 (MQ=21)
                                                                                                                                  gggTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTg                   >  1:3649686/1‑150 (MQ=255)
                                                                                                                                    gTGACCCCGCCCCTCCACGGGCTTGTGTTTCCCTGGTCGGCGCGGGTCCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGc                 <  1:2188382/150‑1 (MQ=255)
                                                                                                                                    gTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGc                 >  1:3335228/1‑150 (MQ=255)
                                                                                                                                     tGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGcc                >  1:1592708/1‑150 (MQ=255)
                                                                                                                                      gAACCCGCTCCTACCAGGGATACACATGCCCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCt               >  2:2162750/1‑150 (MQ=25)
                                                                                                                                      gAACCCGCTCCTACCAGGGATACACATGCCCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCt               >  1:2886952/1‑150 (MQ=25)
                                                                                                                                      gAACCCGCTCCTACCAGGGATACACATGCCCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCt               >  1:2887138/1‑150 (MQ=25)
                                                                                                                                         cccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCg            >  2:2852125/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCg            >  2:3949804/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCg            >  2:3949576/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCg            >  2:2318861/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCg            >  2:156518/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCg            >  2:1143546/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGt           >  2:2095984/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGt           >  2:1049284/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGGCATGTGTATCCCTGGGAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGt           >  1:2096075/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACATGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt         <  2:1264058/150‑1 (MQ=255)
                                                                                                                                            gCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt         >  1:3161428/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt         >  1:3778329/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGttt         >  1:1500916/1‑150 (MQ=255)
                                                                                                                                               ccTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCa      <  2:1822736/150‑1 (MQ=255)
                                                                                                                                                 tACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAgg    >  1:2066738/1‑150 (MQ=255)
                                                                                                                                                 tACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAgg    >  2:1379214/1‑150 (MQ=255)
                                                                                                                                                   aaaGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAGGGt  <  1:1446562/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTGTTGGCACGGCCGTACTGGCCGAAGCGTCGTATCGGAGAGGGATGCTAACGCATGGCGGCCCGTGGACCTAGAGCCGCAGGGGGACTTTGTTGGCCAAAATGAGGGGGTTAGGCCAAAAAGCTTCGCGGGTGAACCCGCTCCTACAAGGGCATGTGTATCCCTGGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGTCATGCGCCCCGTAGGGCGGGTTCATCGGGGCGCCGAACCGCCGCGAACAGACCCGATCAAGCGCCGACGAGGCTCGGTTGGCCTGCGTTTTCAGGGT  >  NC_002947/5726114‑5726410

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