Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 812,051 T→C 7.2% Q68Q (CAA→CAG PP_0698 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947812,0510TC 7.2% 326.6 / 10.5 111Q68Q (CAA→CAGPP_0698LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (52/50);  new base C (2/6);  total (54/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 = 9.50e-01

Only 100 of 113 total aligned reads displayed.

GTGGCATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTC  >  NC_002947/811905‑812200
                                                                                                                                                  |                                                                                                                                                     
gtggCATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGc                                                                                                                                                    >  1:2188779/1‑150 (MQ=255)
    cATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCGTTTCGATGGTCGGTTCGCAAGACAGCACCAGTCCACAGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAggg                                                                                                                                                <  1:1621736/150‑1 (MQ=255)
    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                                                                                                                                     >  2:367471/1‑150 (MQ=255)
               ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcg                                                                                                                                     <  1:815381/150‑1 (MQ=255)
               ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcg                                                                                                                                     >  2:815381/1‑150 (MQ=255)
                gCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcgg                                                                                                                                    >  2:1331990/1‑150 (MQ=255)
                  aGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                  <  2:1385253/150‑1 (MQ=255)
                  aGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                  >  2:418682/1‑150 (MQ=255)
                  aGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                  <  2:767325/150‑1 (MQ=255)
                    gcgcGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGt                                                                                                                                <  2:1586599/150‑1 (MQ=255)
                      gcgGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGa                                                                                                                              <  1:403930/150‑1 (MQ=255)
                      gcgGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCATCAAGTGGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGa                                                                                                                              <  1:667296/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)
                                ccagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGc                                                                                                                    >  2:1143869/1‑150 (MQ=255)
                                ccagcGCATGGCGATGGTCAGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGc                                                                                                                    >  2:1860915/1‑150 (MQ=255)
                                   gcgcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCg                                                                                                                 <  2:1376145/150‑1 (MQ=255)
                                     gcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGt                                                                                                               >  2:186021/1‑150 (MQ=255)
                                               ggTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCTGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                     <  2:502014/150‑1 (MQ=255)
                                               ggTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                     <  2:2081917/150‑1 (MQ=255)
                                                  cagTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                  <  1:1860915/148‑1 (MQ=255)
                                                  cGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                  <  1:612623/150‑1 (MQ=255)
                                                    gTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACa                                                                                                <  1:1661177/150‑1 (MQ=255)
                                                        gCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTc                                                                                            <  2:335654/150‑1 (MQ=255)
                                                          aaGACAGCACCAGCACATTGTCGCGAGTGTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                          >  1:2130825/1‑150 (MQ=255)
                                                          aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                          >  1:163260/1‑150 (MQ=255)
                                                          aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                          >  1:1065820/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)
                                                             acggcaccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGc                                                                                       >  1:2402911/4‑150 (MQ=255)
                                                                 caccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAg                                                                                   >  2:12293/1‑150 (MQ=255)
                                                                  accagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATg                                                                                                                      <  1:583528/114‑1 (MQ=255)
                                                                  accagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAgg                                                                                  <  1:1239840/150‑1 (MQ=255)
                                                                    cagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCt                                                                                >  1:2243605/1‑150 (MQ=255)
                                                                          aTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                          <  2:163260/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)
                                                                                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                                                                   <  2:2402911/150‑1 (MQ=255)
                                                                                 cgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAgg                                                                   <  1:12293/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)
                                                                                        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                                                       <  2:405081/150‑1 (MQ=255)
                                                                                             tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                       >  1:706447/1‑150 (MQ=255)
                                                                                               cgggCGCTGGACGGCGCCGTCGAGCCGTTCGAGGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCt                                                     <  2:802739/148‑1 (MQ=255)
                                                                                                  gcgcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCAc                                                                                         >  1:711723/1‑111 (MQ=255)
                                                                                                  gcgcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCAc                                                                                         <  2:711723/111‑1 (MQ=255)
                                                                                                         aCGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGcc                                                                                                                   >  2:1906658/1‑78 (MQ=255)
                                                                                                         aCGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGcc                                                                                                                   <  1:1906658/78‑1 (MQ=255)
                                                                                                         aCGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTCACCGCGCCg                                           >  1:391767/1‑150 (MQ=255)
                                                                                                            gcgcCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                        <  1:1331990/150‑1 (MQ=255)
                                                                                                              gccctggAGCTGTTCGAAGACCGCCTGGGTGGTGCCTTCCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAAAAGTTGCACGCCAAGGCTTTCTTCCATCGGGCGTATTTGGCGGCTGACCTGGCCgggagg                                      >  2:41300/7‑149 (MQ=255)
                                                                                                                cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                       <  1:1751156/131‑1 (MQ=255)
                                                                                                                cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                       >  2:1751156/1‑131 (MQ=255)
                                                                                                                   cGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGt                                 <  2:378701/150‑1 (MQ=255)
                                                                                                                    gAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                  <  1:2362251/52‑1 (MQ=255)
                                                                                                                    gAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCgggg                                                                                                                                  >  2:2362251/1‑52 (MQ=255)
                                                                                                                      gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGa                                                                                                                                           >  1:2293475/1‑41 (MQ=255)
                                                                                                                      gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGa                                                                                                                                           <  2:2293475/41‑1 (MQ=255)
                                                                                                                      gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGTCGt                                                                >  1:1915223/1‑116 (MQ=255)
                                                                                                                      gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGTCGt                                                                <  2:1915223/116‑1 (MQ=255)
                                                                                                                             gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                       >  1:1127349/1‑150 (MQ=255)
                                                                                                                             gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                       >  1:823670/1‑150 (MQ=255)
                                                                                                                             gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                       <  1:186021/150‑1 (MQ=255)
                                                                                                                             gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGACCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                       <  2:1335561/150‑1 (MQ=255)
                                                                                                                              aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcggc                                                                                                         <  1:759700/67‑1 (MQ=255)
                                                                                                                              aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcggc                                                                                                         >  2:759700/1‑67 (MQ=255)
                                                                                                                                     ccTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAAc                                                                                                 <  1:1434154/68‑1 (MQ=255)
                                                                                                                                     ccTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAAc                                                                                                 >  2:1434154/1‑68 (MQ=255)
                                                                                                                                      cTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAGATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                        <  1:504505/124‑1 (MQ=255)
                                                                                                                                      cTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAGATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                        >  2:504505/1‑124 (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                                                                                                  <  2:211481/57‑1 (MQ=255)
                                                                                                                                               gCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                  >  1:211481/1‑57 (MQ=255)
                                                                                                                                                ccTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGctc    >  2:326169/1‑150 (MQ=255)
                                                                                                                                                  ttGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCtctc  >  1:310328/1‑150 (MQ=255)
                                                                                                                                                  ttGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCtctc  <  2:1065820/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGGCATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTC  >  NC_002947/811905‑812200

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: