Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,165,415 Δ1 bp 6.3% coding (2513/3906 nt) hrpA ← ATP‑dependent helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,165,4150T. 6.3% 374.0 / 18.8 127coding (2513/3906 nt)hrpAATP‑dependent helicase HrpA
Reads supporting (aligned to +/- strand):  ref base T (51/67);  new base . (4/4);  total (55/72)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

Only 100 of 128 total aligned reads displayed.

CCGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGCC  >  NC_002947/5165269‑5165564
                                                                                                                                                  |                                                                                                                                                      
ccGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgc                                                                                                                                                     >  2:2768858/1‑150 (MQ=255)
ccGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgc                                                                                                                                                     >  2:611836/1‑150 (MQ=255)
ccGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgc                                                                                                                                                     <  1:3032287/150‑1 (MQ=255)
 cGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACCGCTTCGCGGGCGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgcg                                                                                                                                                    <  2:858977/150‑1 (MQ=255)
 cGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgcg                                                                                                                                                    <  2:3010060/150‑1 (MQ=255)
 cGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGCACTGTGCGGCGGTCACTTCGCTGGCCTcgcg                                                                                                                                                    <  2:3009826/150‑1 (MQ=255)
  ggCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCgg                                                                                                                                                  <  1:350635/150‑1 (MQ=255)
  ggCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCgg                                                                                                                                                  <  1:2567181/150‑1 (MQ=255)
    cAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGc                                                                                                                                                <  2:2393414/150‑1 (MQ=255)
         tCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCa                                                                                                                                           >  2:1491074/1‑150 (MQ=255)
          cGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCAc                                                                                                                                          <  1:193465/150‑1 (MQ=255)
             cAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGtc                                                                                                                                       <  1:3061301/150‑1 (MQ=255)
             cAACAGCGGTGCCGGCACCCGCACCGGCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcgaagt                                                                                                                                       <  2:1477552/150‑6 (MQ=255)
                tgacGGTGCGGGTGCCGGCACCGTGCCACCATCGGGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGTGGGCCAGCACGtcctc                                                                                                                                    >  2:508271/4‑150 (MQ=255)
                 aGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGtcctcc                                                                                                                                   >  1:772473/1‑150 (MQ=255)
                   cGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTACATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCAAGCGGGCCAGACCGtcctccct                                                                                                                                 <  2:821403/150‑3 (MQ=255)
                     gTGCCGGCCCCCGCACCGTCACACCATCACGGGGGTGCTTCGTCGCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTcgc                                                                                                                               <  1:2337582/150‑1 (MQ=255)
                     gTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTcgc                                                                                                                               >  2:200561/1‑150 (MQ=255)
                     gTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGCCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTcgc                                                                                                                               <  2:2396572/150‑1 (MQ=255)
                      tGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTcgcg                                                                                                                              <  2:2529793/150‑1 (MQ=255)
                               ccGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATCCTGTCGACGTACGGTCCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCAGGGTcag                                                                                                                     >  2:576279/1‑150 (MQ=255)
                                 gCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcagca                                                                                                                   <  1:1293615/150‑1 (MQ=255)
                                  ccccGTCACGCCACGACGGGGGTGGCCGGGCATGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcagcag                                                                                                                  <  2:387388/148‑1 (MQ=255)
                                      gTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                              <  1:2117326/150‑1 (MQ=255)
                                      gTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                              <  2:2692103/150‑1 (MQ=255)
                                      gTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                              >  1:1299438/1‑150 (MQ=255)
                                      gTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGGCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                              >  1:1451201/1‑150 (MQ=255)
                                        cacaCCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGc                                                                                                            <  1:1502292/150‑1 (MQ=255)
                                           acCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCTCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTc                                                                                                         >  2:666200/1‑150 (MQ=255)
                                                cACGGGGGTGGCCGGGCTCGACGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCtttt                                                                                                    <  1:200561/150‑1 (MQ=255)
                                                 aCGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGt                                                                                                                >  2:1826540/1‑137 (MQ=255)
                                                 aCGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGt                                                                                                                <  1:1826540/137‑1 (MQ=255)
                                                  cGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTgg                                                                                                  <  2:1299438/150‑1 (MQ=255)
                                                         ggCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCa                                                                                           <  1:2768858/150‑1 (MQ=255)
                                                                 tCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGa                                                                                                                         >  2:392501/1‑112 (MQ=255)
                                                                 tCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGa                                                                                                                         <  1:392501/112‑1 (MQ=255)
                                                                 tCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATca                                                                                                                      >  2:3001088/1‑115 (MQ=255)
                                                                 tCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTAc                                                                                   <  2:776933/150‑1 (MQ=255)
                                                                     aGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAAc                                                                               >  2:1912435/1‑150 (MQ=255)
                                                                     aGTGGCAGCTCATCGGTAGGCTCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTTGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGCTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAAc                                                                               <  2:2090277/150‑1 (MQ=255)
                                                                          tAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCa                                                                          >  2:1326956/1‑150 (MQ=255)
                                                                           aGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCaa                                                                         >  1:340992/1‑150 (MQ=255)
                                                                            gCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACGTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCaaa                                                                        <  1:228540/150‑1 (MQ=255)
                                                                             cTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAg                                                                       >  1:2175448/1‑150 (MQ=255)
                                                                              tCAGCGGCAGGCGCAGATCGCCCACCTGCCTACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAgg                                                                      >  1:2860760/1‑150 (MQ=255)
                                                                              tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGGATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAgg                                                                      >  2:931378/1‑150 (MQ=255)
                                                                              tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCtt                                                                                                      >  1:1082765/1‑118 (MQ=255)
                                                                              tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCtt                                                                                                      <  2:1082765/118‑1 (MQ=255)
                                                                              tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAgg                                                                      <  1:1491074/150‑1 (MQ=255)
                                                                              tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGCACTGTGCGGCGGTCACTTCGCTGGCCTGGCGGGCCAGCACGTCCTCTTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAgg                                                                      >  1:565337/1‑150 (MQ=255)
                                                                                  cGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcg                                                                                                                                           >  1:614018/1‑76 (MQ=255)
                                                                                   ggCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGcc                                                                 <  1:714422/150‑1 (MQ=255)
                                                                                        gcgcAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtg                                                            <  1:2672244/150‑1 (MQ=255)
                                                                                          gcAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtggt                                                          <  2:2364066/150‑1 (MQ=255)
                                                                                           cAGATCGCCCACCTTCATACTGTCGGGGCACCGTGCGGCTGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGTTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtggtg                                                         <  2:565337/150‑1 (MQ=255)
                                                                                             gATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcag                                                                                                                     >  2:2884401/1‑88 (MQ=255)
                                                                                             gATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcag                                                                                                                     <  1:2884401/88‑1 (MQ=255)
                                                                                                  cccACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTTCTTGTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTcc                                                  <  1:2679559/150‑1 (MQ=255)
                                                                                                     cccTGCATACTGTCGGGGTACTGTGAGGCGGTCACTTAGCTAGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTccgcc                                               <  2:2175448/149‑1 (MQ=255)
                                                                                                     aCCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTccgcc                                               >  1:1008720/1‑150 (MQ=255)
                                                                                                     aCCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcgaa                                                                                                                                         <  2:1167964/60‑4 (MQ=25)
                                                                                                     aCCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcgaa                                                                                                                                         >  1:1167964/1‑57 (MQ=25)
                                                                                                         gCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGTTGACCTCGCGGGCCACAACGTCCTTCTAGCTCTTGATCAGCAGGTTGTCGTCCTTTTTGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTTGATTTCCGCCGGCa                                           <  2:1119640/150‑1 (MQ=255)
                                                                                                          cATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTTGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAg                                          >  1:1022719/1‑150 (MQ=255)
                                                                                                            tACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGGATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                        <  1:931378/150‑1 (MQ=255)
                                                                                                            tACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                        >  2:2068118/1‑150 (MQ=255)
                                                                                                            tACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                        <  1:3071883/150‑1 (MQ=255)
                                                                                                              cTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgcgc                                      <  1:615704/150‑1 (MQ=255)
                                                                                                               tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCgg                                             <  2:725257/142‑1 (MQ=255)
                                                                                                               tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCgg                                             >  1:725257/1‑142 (MQ=255)
                                                                                                               tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgcgcg                                     >  1:1234454/1‑150 (MQ=255)
                                                                                                                 tCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCt                                   <  1:2951524/150‑1 (MQ=255)
                                                                                                                  cGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCtt                                  <  2:2528258/150‑1 (MQ=255)
                                                                                                                    gggTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCa                                <  1:666200/150‑1 (MQ=255)
                                                                                                                      ggaCTGTGCGGCGGTCACTTCGCTGTCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATa                              <  2:710782/148‑1 (MQ=255)
                                                                                                                          tgtgCTGCGGTCACTTCTCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTCCCGCCGGCAGGCGCGCTTCATAGaaa                          <  1:576279/150‑1 (MQ=255)
                                                                                                                              cggcggTCACTTCGCTGGCCTCGCGGGCCAGCGCGCGCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGt                      <  2:1022719/150‑1 (MQ=255)
                                                                                                                               ggcggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTg                                                                                                   <  2:679465/72‑1 (MQ=255)
                                                                                                                               ggcggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTg                                                                                                   >  1:679465/1‑72 (MQ=255)
                                                                                                                               ggcggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTa                     >  2:243104/1‑150 (MQ=255)
                                                                                                                                  ggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGt                      <  1:495067/146‑1 (MQ=255)
                                                                                                                                   gTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGc                 >  1:2396835/1‑150 (MQ=255)
                                                                                                                                      gacgtgCTGGCCCGCGAGGCCAGCGCGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCt                                   <  2:934982/124‑1 (MQ=255)
                                                                                                                                      gacgtgCTGGCCCGCGAGGCCAGCGCGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCt                                   >  1:934982/6‑129 (MQ=255)
                                                                                                                                      aCTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATa                              <  2:280930/134‑1 (MQ=255)
                                                                                                                                        ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                              <  1:2986572/52‑1 (MQ=255)
                                                                                                                                        ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc            >  2:1657943/1‑150 (MQ=255)
                                                                                                                                        ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc            >  2:1539397/1‑150 (MQ=255)
                                                                                                                                        ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc            >  2:313094/1‑150 (MQ=255)
                                                                                                                                         gtgCTGGCCCGCGAGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttct           >  1:1831635/3‑150 (MQ=255)
                                                                                                                                          cGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTa                                                                             <  2:47861/83‑1 (MQ=255)
                                                                                                                                          cGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTa                                                                             >  1:47861/1‑83 (MQ=255)
                                                                                                                                           gCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttcttc         <  1:1657943/150‑1 (MQ=255)
                                                                                                                                              ggCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTc      <  1:320813/150‑1 (MQ=255)
                                                                                                                                               gCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                        <  1:911974/115‑1 (MQ=255)
                                                                                                                                               gCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                        >  2:911974/1‑115 (MQ=255)
                                                                                                                                               gCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCa     <  2:1656185/150‑1 (MQ=255)
                                                                                                                                                 cTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttcttc         >  1:2198579/1‑144 (MQ=255)
                                                                                                                                                 cTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttcttc         <  2:2198579/144‑1 (MQ=255)
                                                                                                                                                  tCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGcc  >  2:235805/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CCGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGCC  >  NC_002947/5165269‑5165564

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