Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 812,056 G→A 5.9% L67L (CTG→TTG)  PP_0698 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947812,0560GA 5.9% 344.3 / 6.6 114L67L (CTG→TTG) PP_0698LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (55/51);  new base A (2/5);  total (57/57)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

CATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCGGTC  >  NC_002947/811909‑812204
                                                                                                                                                   |                                                                                                                                                    
cATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAagg                                                                                                                                                    <  2:158258/150‑1 (MQ=255)
 aTGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAgggg                                                                                                                                                   >  2:1184348/1‑149 (MQ=255)
      aGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCTTGGGTCGTGCCCTGCAAGGcaccac                                                                                                                                              <  1:68953/150‑1 (MQ=255)
      aGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCTGCAAGGcaccac                                                                                                                                              <  2:891362/150‑1 (MQ=255)
       gcgcGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGAcc                                                                                                                                             >  2:1661177/1‑150 (MQ=255)
           ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcg                                                                                                                                         <  1:815381/150‑1 (MQ=255)
           ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcg                                                                                                                                         >  2:367471/1‑150 (MQ=255)
           ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcg                                                                                                                                         >  2:815381/1‑150 (MQ=255)
           ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccacccagg                                                                                                                                         <  2:1438280/150‑5 (MQ=255)
            gCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcgg                                                                                                                                        >  2:1331990/1‑150 (MQ=255)
              aGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                      >  2:418682/1‑150 (MQ=255)
              aGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                      <  2:1385253/150‑1 (MQ=255)
              aGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                      <  2:767325/150‑1 (MQ=255)
                gcgcGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGt                                                                                                                                    <  2:1586599/150‑1 (MQ=255)
                  gcgGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGa                                                                                                                                  <  1:403930/150‑1 (MQ=255)
                   cgGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAg                                                                                                                                 <  2:142019/150‑1 (MQ=255)
                     ggATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcgg                                                                                                                                        <  1:905891/141‑1 (MQ=255)
                        tCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGATCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCTGCAAGGCAcca                                                                                                                                               <  1:1141955/131‑3 (MQ=255)
                        tCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCTGCAAGGCAcca                                                                                                                                               >  2:1141955/1‑129 (MQ=255)
                            ccagcGCATGGCGATGGTCAGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGc                                                                                                                        >  2:1860915/1‑150 (MQ=255)
                                 gcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGt                                                                                                                   >  2:186021/1‑150 (MQ=255)
                                           ggTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCTGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                         <  2:502014/150‑1 (MQ=255)
                                           ggTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                         <  2:2081917/150‑1 (MQ=255)
                                              cGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                      <  1:767475/150‑1 (MQ=255)
                                              cGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                      <  1:612623/150‑1 (MQ=255)
                                                gTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACa                                                                                                    <  1:1661177/150‑1 (MQ=255)
                                                  tCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGt                                                                                                  <  1:1796616/150‑1 (MQ=255)
                                                    gCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTc                                                                                                <  2:335654/150‑1 (MQ=255)
                                                      aaGACAGCACCAGCACATTGTCGCGAGTGTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                              >  1:2130825/1‑150 (MQ=255)
                                                      aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                              >  1:1065820/1‑150 (MQ=255)
                                                      aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                              >  1:163260/1‑150 (MQ=255)
                                                      aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTAGAGCTGTTCGAAGATGGCCTGGGGGGGGCCTTGGAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                              >  1:899543/1‑150 (MQ=255)
                                                        gACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACg                                                                                            >  1:2036016/1‑150 (MQ=255)
                                                        gACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccaccca                                                                                                                                           >  1:299017/1‑101 (MQ=38)
                                                        gACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccaccca                                                                                                                                           <  2:299017/103‑3 (MQ=38)
                                                             caccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAg                                                                                       >  2:12293/1‑150 (MQ=255)
                                                              accagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATg                                                                                                                          >  2:583528/1‑114 (MQ=255)
                                                              accagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATg                                                                                                                          <  1:583528/114‑1 (MQ=255)
                                                                cagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCt                                                                                    >  1:1420380/1‑150 (MQ=255)
                                                                cagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCt                                                                                    >  1:2243605/1‑150 (MQ=255)
                                                                  gcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCttt                                                                                  <  2:1183062/150‑1 (MQ=255)
                                                                      aTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGATTAGCATGGCACGACCGGCGGGGTTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                              <  1:1496386/150‑1 (MQ=255)
                                                                       ttGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTGCTTcc                                                                             >  1:104335/1‑150 (MQ=255)
                                                                       ttGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGACCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTcc                                                                             >  1:1335561/1‑150 (MQ=255)
                                                                            gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGTAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCaagcag                                                                        >  1:274714/1‑150 (MQ=255)
                                                                            gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCaa                                                                                        >  2:1968821/1‑134 (MQ=255)
                                                                            gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCaa                                                                                        <  1:1968821/134‑1 (MQ=255)
                                                                            gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcagcag                                                                        <  1:367471/150‑1 (MQ=255)
                                                                            gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTCATGATCAGGTCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcagcag                                                                        <  1:150356/150‑1 (MQ=255)
                                                                             cgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAgg                                                                       <  1:12293/150‑1 (MQ=255)
                                                                             cgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAgg                                                                       <  2:2402911/150‑1 (MQ=255)
                                                                              gAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCggg                                                                                                                                       <  2:35131/85‑1 (MQ=255)
                                                                              gAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCggg                                                                                                                                       >  1:35131/1‑85 (MQ=255)
                                                                                 cTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTa                                                                   >  1:1654508/1‑150 (MQ=255)
                                                                                  ttgttgATGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTAt                                                                  >  2:31435/1‑150 (MQ=255)
                                                                                    gttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCCAGGCTTTCTTCCGGCAGGCGTAttt                                                                <  2:216262/150‑1 (MQ=255)
                                                                                    gttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTAttt                                                                >  1:192762/1‑150 (MQ=255)
                                                                                     ttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTg                                                               >  2:581654/1‑150 (MQ=255)
                                                                                     ttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTg                                                               >  2:2007563/1‑150 (MQ=255)
                                                                                      tgctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTgg                                                              <  1:581654/150‑1 (MQ=255)
                                                                                         tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                           >  1:706447/1‑150 (MQ=255)
                                                                                         tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                           <  2:405081/150‑1 (MQ=255)
                                                                                           cgggCGCTGGACGGCGCCGTCGAGCCGTTCGAGGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCt                                                         <  2:802739/148‑1 (MQ=255)
                                                                                             ggCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGa                                                       >  2:362123/1‑150 (MQ=255)
                                                                                              gcgcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCAc                                                                                             >  1:711723/1‑111 (MQ=255)
                                                                                                gcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACcg                                                    >  1:2112379/1‑150 (MQ=255)
                                                                                                     aCGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGcc                                                                                                                       <  1:1906658/78‑1 (MQ=255)
                                                                                                        gcgcCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                            <  1:1331990/150‑1 (MQ=255)
                                                                                                            cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                           <  1:1751156/131‑1 (MQ=255)
                                                                                                            cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                           >  2:1751156/1‑131 (MQ=255)
                                                                                                            cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtgagtga                                        >  1:519461/1‑150 (MQ=255)
                                                                                                               cGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGt                                     <  2:378701/150‑1 (MQ=255)
                                                                                                                gAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                      >  2:2362251/1‑52 (MQ=255)
                                                                                                                gAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                      <  1:2362251/52‑1 (MQ=255)
                                                                                                                  gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGa                                                                                                                                               <  2:2293475/41‑1 (MQ=255)
                                                                                                                  gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGa                                                                                                                                               >  1:2293475/1‑41 (MQ=255)
                                                                                                                  gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGTCGt                                                                    <  2:1915223/116‑1 (MQ=255)
                                                                                                                  gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGTCGt                                                                    >  1:1915223/1‑116 (MQ=255)
                                                                                                                         gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                           <  1:186021/150‑1 (MQ=255)
                                                                                                                         gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                           >  1:1127349/1‑150 (MQ=255)
                                                                                                                         gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGACCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                           <  2:1335561/150‑1 (MQ=255)
                                                                                                                          aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcggc                                                                                                             >  2:759700/1‑67 (MQ=255)
                                                                                                                          aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcggc                                                                                                             <  1:759700/67‑1 (MQ=255)
                                                                                                                                 ccTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAAc                                                                                                     >  2:1434154/1‑68 (MQ=255)
                                                                                                                                 ccTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAAc                                                                                                     <  1:1434154/68‑1 (MQ=255)
                                                                                                                                  cTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAGATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                            >  2:504505/1‑124 (MQ=255)
                                                                                                                                  cTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAGATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                            <  1:504505/124‑1 (MQ=255)
                                                                                                                                        ggtgCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGat            <  2:1654508/150‑2 (MQ=255)
                                                                                                                                        ggtgCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGaa            <  2:706447/150‑1 (MQ=255)
                                                                                                                                         gtgCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGgcgc                   >  2:1998934/1‑142 (MQ=255)
                                                                                                                                         gtgCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGgcgc                   <  1:1998934/142‑1 (MQ=255)
                                                                                                                                          tgCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAACTGGCCCGGTTGCGGCGTTCGAAGAGTTCCCCGCCAGCGCGTGCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGc          >  1:295141/1‑150 (MQ=255)
                                                                                                                                           gCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                      >  1:211481/1‑57 (MQ=255)
                                                                                                                                              ttGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCtctc      <  2:1065820/150‑1 (MQ=255)
                                                                                                                                              ttGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCtctc      >  1:310328/1‑150 (MQ=255)
                                                                                                                                               tGCAGGGCCCGCCAGGCGGGGGTGCGGAAAATCGCTCTTGTGCTGTGGTGGAACTGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCg     <  2:2070607/150‑1 (MQ=255)
                                                                                                                                               tGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcagcag                                                                        <  2:1079105/83‑1 (MQ=255)
                                                                                                                                               tGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcagcag                                                                        >  1:1079105/1‑83 (MQ=255)
                                                                                                                                                  aGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCGGTc  >  2:1808655/1‑150 (MQ=255)
                                                                                                                                                  aGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCGGTc  >  2:1239848/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCGGTC  >  NC_002947/811909‑812204

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: