Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 563,007 T→G 6.7% P30P (CCT→CCG) ‡ rpoA → DNA‑directed RNA polymerase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947563,0070TG 6.7% 228.5 / 5.5 74P30P (CCT→CCG) ‡rpoADNA‑directed RNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (32/37);  new base G (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAATT  >  NC_002947/562860‑563154
                                                                                                                                                   |                                                                                                                                                    
aaCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCctc                                                                                                                                                     >  8:119884/1‑149 (MQ=255)
 aCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCctct                                                                                                                                                    >  4:160529/1‑149 (MQ=255)
   gAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGATCGAGCCtctc                                                                                                                                                   >  4:167042/1‑148 (MQ=255)
      aaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCAAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCACCGACACATTGAAGTGCAGGTAGACAGTCCGACCCGCGCCAAAATTACTCTCGATCATCTCGAGc                                                                                                                                               >  4:201129/1‑149 (MQ=255)
        ccTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGt                                                                                                                                             <  5:2989/149‑1 (MQ=255)
         cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           <  1:108666/149‑1 (MQ=255)
         cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           <  3:160529/149‑1 (MQ=255)
         cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           >  2:5292/1‑149 (MQ=255)
            aTTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGtt                                                                                                                                        <  1:164287/149‑1 (MQ=255)
                tCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCgg                                                                                                                                    >  3:146146/1‑149 (MQ=255)
                     ctctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATa                                                                                                                               >  7:212138/1‑149 (MQ=255)
                       ctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAGATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAcc                                                                                                                             >  7:159331/1‑149 (MQ=255)
                        tACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAccc                                                                                                                            >  1:282991/1‑149 (MQ=255)
                          cTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCGCATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCt                                                                                                                           <  7:77787/148‑1 (MQ=255)
                          cTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTg                                                                                                                          <  8:212138/149‑1 (MQ=255)
                           tCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTgg                                                                                                                         >  2:250430/1‑149 (MQ=255)
                           tCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTgg                                                                                                                         >  2:278525/1‑149 (MQ=255)
                             cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTgg                                                                                                                         <  3:145058/147‑1 (MQ=255)
                             cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGc                                                                                                                       <  2:99929/149‑1 (MQ=255)
                               aGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCaa                                                                                                                     <  6:191193/149‑1 (MQ=255)
                                   aGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAAcgcg                                                                                                                 <  4:146146/149‑1 (MQ=255)
                                          gAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGTTGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGcgcc                                                                                                          <  1:250430/149‑1 (MQ=255)
                                          gAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGcgcc                                                                                                          >  3:11862/1‑149 (MQ=255)
                                          caaaaTAGTTGAATCCATGAAGAATTCGGTAAATCAGTACCTGACTACCCCCCACATTGAAGTGCAGGAAGTCAGTCCGACCCGCGCCAAAATAACGATCTAGCATCTCGACCGAGTATACGGCAATACCCTGCTCAACGAGAAGcgcc                                                                                                          >  8:223477/2‑149 (MQ=255)
                                                  gttcATCCATGCAGATTTCGTTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGCACCGCGCCAAAATTACGCTCGAGACTCTCGAGCGTGGTTTCGGCCAGACCCTGGGCAACGCGCTGCGCCGCATccag                                                                                                  >  6:270962/4‑147 (MQ=255)
                                                  gTGCATCCATGCAGATTTCGGTACATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtg                                                                                                  <  8:159331/149‑1 (MQ=255)
                                                    gCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAAAGCGCTGCGCCGCATCCtgtt                                                                                                >  1:99909/1‑149 (MQ=255)
                                                     catgcatgCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgtt                                                                                                <  5:248655/148‑1 (MQ=255)
                                                     cATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgtt                                                                                                <  7:119884/148‑1 (MQ=255)
                                                      aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttgt                                                                                              >  7:14690/1‑149 (MQ=255)
                                                      aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttgt                                                                                              >  5:45161/1‑149 (MQ=255)
                                                       tCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtc                                                                                             <  4:41773/149‑1 (MQ=255)
                                                         cATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcct                                                                                           <  1:278525/149‑1 (MQ=255)
                                                              aGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTTAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATg                                                                                      >  7:235940/1‑149 (MQ=255)
                                                                  ttCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTTAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCt                                                                                   <  8:235940/148‑1 (MQ=255)
                                                                  ttCGGTAAATGAGTTCCTGACTCACCTCCACATTGATTTGCAGGTTGTCAGTCCGACCCGCGCCAAAATTACGCTGGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCt                                                                                   <  3:201129/148‑1 (MQ=255)
                                                                     ggTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGc                                                                                >  4:156532/1‑148 (MQ=255)
                                                                     ggTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGc                                                                                <  3:156532/148‑1 (MQ=255)
                                                                         aaTGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCAAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgtg                                                                            <  6:45161/148‑1 (MQ=255)
                                                                          aTGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCa                                                                          <  4:11862/149‑1 (MQ=255)
                                                                             aGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCagta                                                                       <  3:167042/149‑1 (MQ=255)
                                                                              gTTCCTTACTCCACGCCACATTTCTGTGCAGGTAGTCAGTCCGACCCTCGCCAAAATTACGCTGGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCGTGGCTGTGCagtag                                                                      <  7:223477/149‑1 (MQ=255)
                                                                                 ccTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                    >  6:150994/1‑148 (MQ=255)
                                                                                  cTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                    >  4:145142/1‑147 (MQ=255)
                                                                                         cccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTg                                                                                  <  4:244969/126‑1 (MQ=255)
                                                                                         cccGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTg                                                                                  >  3:244969/1‑126 (MQ=255)
                                                                                               acaTTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGAt                                                     >  1:132874/1‑149 (MQ=255)
                                                                                                caTTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATg                                                    <  3:145142/149‑1 (MQ=255)
                                                                                                               agtcagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGAGTACTcatggag                                      >  8:112497/1‑142 (MQ=255)
                                                                                                                 tcagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTAc                                   >  1:172292/1‑149 (MQ=255)
                                                                                                                   agtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTc                                 >  4:141459/1‑149 (MQ=255)
                                                                                                                     tcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcg                               <  5:150994/149‑1 (MQ=255)
                                                                                                                     tcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGAGTACTCATGGAGTACGcca                               <  7:112497/149‑2 (MQ=255)
                                                                                                                       cGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcgc                              >  2:233291/1‑148 (MQ=255)
                                                                                                                             gcgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCg                         >  1:192727/1‑148 (MQ=255)
                                                                                                                              cgcCAAAATTAGGCTCGAGCCTCTCGTGCGTTGTTTGGGCCATTCCCTGGGCAACGCGCTGCGCCGCATCCTGTTTTCCTCCATGCCTGTCTGTGCAGTAGTCGAGGCCTAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                      <  5:270962/149‑1 (MQ=255)
                                                                                                                              cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                      <  5:259532/149‑1 (MQ=255)
                                                                                                                              cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                      <  4:6663/149‑1 (MQ=255)
                                                                                                                              cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                      <  2:192727/149‑1 (MQ=255)
                                                                                                                              cgcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAg                      <  7:279810/149‑1 (MQ=255)
                                                                                                                                ccAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTAGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgt                    >  2:113290/1‑149 (MQ=255)
                                                                                                                                 cgaaaccaCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCCTCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtg                   >  7:207931/8‑149 (MQ=255)
                                                                                                                                 cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgt                                                                             <  3:49677/91‑1 (MQ=255)
                                                                                                                                 cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCtgt                                                                             >  4:49677/1‑91 (MQ=255)
                                                                                                                                 cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATc                                                        >  7:123583/1‑112 (MQ=255)
                                                                                                                                 cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATc                                                        <  8:123583/112‑1 (MQ=255)
                                                                                                                                 cAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtg                   <  3:141459/149‑1 (MQ=255)
                                                                                                                                  aaaaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtgt                  <  2:282991/149‑1 (MQ=255)
                                                                                                                                   aaaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGgtgt                  >  4:139607/1‑148 (MQ=255)
                                                                                                                                      ttACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg              >  6:204193/1‑149 (MQ=255)
                                                                                                                                       caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGc                                                  >  2:307949/2‑112 (MQ=255)
                                                                                                                                       caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGc                                                  <  1:307949/111‑1 (MQ=255)
                                                                                                                                        aCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTGCATGAGTACTCCGCGATCGAAGGTGTTCAg              <  6:197816/147‑1 (MQ=255)
                                                                                                                                        aCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCAGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg              >  5:197816/1‑147 (MQ=255)
                                                                                                                                          gCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGa                                       <  6:122846/120‑1 (MQ=255)
                                                                                                                                          gCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGa                                       >  5:122846/1‑120 (MQ=255)
                                                                                                                                             cGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTTGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTACGCGATCGAAGGGGTTCAGGAAGAc        >  8:269999/1‑148 (MQ=255)
                                                                                                                                                   tctcGAGCGTGGTTTCGGCCATACACTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAAtt  >  5:115613/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAATT  >  NC_002947/562860‑563154

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