Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,165,419:1 +A 6.7% coding (2509/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,4191.A 6.7% 473.3 / 27.1 119coding (2509/3906 nt)hrpAATP‑dependent helicase HrpA
Reads supporting (aligned to +/- strand):  ref base . (50/61);  new base A (4/4);  total (54/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 126 total aligned reads displayed.

CGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGCCAGGA  >  NC_002947/5165270‑5165568
                                                                                                                                                      |                                                                                                                                                     
cGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACCGCTTCGCGGGCGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgcg                                                                                                                                                        <  2:858977/150‑1 (MQ=255)
cGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgcg                                                                                                                                                        <  2:3010060/150‑1 (MQ=255)
 ggCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCgg                                                                                                                                                      <  1:2567181/150‑1 (MQ=255)
 ggCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCgg                                                                                                                                                      <  1:350635/150‑1 (MQ=255)
   cAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGc                                                                                                                                                    <  2:2393414/150‑1 (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)
                   ggTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTcg                                                                                                                                    <  2:2419800/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)
                                gCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcagca                                                                                                                       <  1:1293615/150‑1 (MQ=255)
                                 ccccGTCACGCCACGACGGGGGTGGCCGGGCATGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcagcag                                                                                                                      <  2:387388/148‑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)
                                                aCGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTg                                                                                                       <  2:1506416/150‑1 (MQ=255)
                                                 cGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTgg                                                                                                      <  2:1299438/150‑1 (MQ=255)
                                                        ggCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCa                                                                                               <  1:2768858/150‑1 (MQ=255)
                                                               cTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTa                                                                                        <  1:611836/150‑1 (MQ=255)
                                                                tCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATAGTGTCGGGGTACTGTGCGGCGGGCACTTGGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGGTCAGCAGGTTGGCGTCCTTTTGGCTTCCCATGCGGTAc                                                                                       >  1:323449/1‑150 (MQ=255)
                                                                tCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATca                                                                                                                          >  2:3001088/1‑115 (MQ=255)
                                                                tCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATca                                                                                                                          <  1:3001088/115‑1 (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                                                                                                          <  2:1082765/118‑1 (MQ=255)
                                                                             tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCtt                                                                                                          >  1:1082765/1‑118 (MQ=255)
                                                                             tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAgg                                                                          <  1:1491074/150‑1 (MQ=255)
                                                                             tCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGCACTGTGCGGCGGTCACTTCGCTGGCCTGGCGGGCCAGCACGTCCTCTTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAgg                                                                          >  1:565337/1‑150 (MQ=255)
                                                                                 cGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcg                                                                                                                                               <  2:614018/77‑2 (MQ=255)
                                                                                 cGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcg                                                                                                                                               >  1:614018/1‑76 (MQ=255)
                                                                                  ggCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGcc                                                                     <  1:714422/150‑1 (MQ=255)
                                                                                       gcgcAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtg                                                                <  1:2672244/150‑1 (MQ=255)
                                                                                        cgcAGATCGCCCACCTGCATACTGTCGGGGTACTGGGAGGAGGTCACTTCGATGGCCTCGCGGGCCAGCACGTCATCATCGCGCATGATGAGAAGGTTGGCGTAATTTTGGCTGAATAAGGGATAAAAACTATAAAAGGTGGAGGTCtgg                                                               >  2:1227602/1‑150 (MQ=255)
                                                                                          cAGATCGCCCACCTTCATACTGTCGGGGCACCGTGCGGCTGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGTTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtggtg                                                             <  2:565337/150‑1 (MQ=255)
                                                                                            gATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcag                                                                                                                         >  2:2884401/1‑88 (MQ=255)
                                                                                                    cccTGCATACTGTCGGGGTACTGTGAGGCGGTCACTTAGCTAGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTccgcc                                                   <  2:2175448/149‑1 (MQ=255)
                                                                                                    aCCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTccgcc                                                   >  1:1008720/1‑150 (MQ=255)
                                                                                                    aCCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcgaa                                                                                                                                             >  1:1167964/1‑57 (MQ=25)
                                                                                                        gCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGTTGACCTCGCGGGCCACAACGTCCTTCTAGCTCTTGATCAGCAGGTTGTCGTCCTTTTTGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTTGATTTCCGCCGGCa                                               <  2:1119640/150‑1 (MQ=255)
                                                                                                        gCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCa                                               >  1:431956/1‑150 (MQ=255)
                                                                                                         cATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTTGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAg                                              >  1:1022719/1‑150 (MQ=255)
                                                                                                           tACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGGATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                            <  1:931378/150‑1 (MQ=255)
                                                                                                           tACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                            >  2:2068118/1‑150 (MQ=255)
                                                                                                             cTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgcgc                                          <  1:615704/150‑1 (MQ=255)
                                                                                                              tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCgg                                                 >  1:725257/1‑142 (MQ=255)
                                                                                                              tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCgg                                                 <  2:725257/142‑1 (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)
                                                                                                                     gTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATa                                  <  1:2087202/150‑1 (MQ=255)
                                                                                                                           tgCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCAt                                   <  2:2658115/143‑1 (MQ=255)
                                                                                                                           tgCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCAt                                   >  1:2658115/1‑143 (MQ=255)
                                                                                                                             cggcggTCACTTCGCTGGCCTCGCGGGCCAGCGCGCGCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGt                          <  2:1022719/150‑1 (MQ=255)
                                                                                                                              ggcggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTg                                                                                                       >  1:679465/1‑72 (MQ=255)
                                                                                                                              ggcggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTg                                                                                                       <  2:679465/72‑1 (MQ=255)
                                                                                                                                 ggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGt                          >  2:495067/1‑146 (MQ=255)
                                                                                                                                 ggTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGt                          <  1:495067/146‑1 (MQ=255)
                                                                                                                                  gTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGc                     >  1:2396835/1‑150 (MQ=255)
                                                                                                                                  gTCACTTCGCTGGCCTCGCGGGCCAGCACGTCATCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGc                     >  1:2254679/1‑150 (MQ=255)
                                                                                                                                     gacgtgCTGGCCCGCGAGGCCAGCGCGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCt                                       <  2:934982/124‑1 (MQ=255)
                                                                                                                                     gacgtgCTGGCCCGCGAGGCCAGCGCGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCt                                       >  1:934982/6‑129 (MQ=255)
                                                                                                                                     aCTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATa                                  >  1:280930/1‑134 (MQ=255)
                                                                                                                                     aCTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATa                                  <  2:280930/134‑1 (MQ=255)
                                                                                                                                       ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                                  <  1:2986572/52‑1 (MQ=255)
                                                                                                                                       ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                                  >  2:2986572/1‑52 (MQ=255)
                                                                                                                                       ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc                >  2:313094/1‑150 (MQ=255)
                                                                                                                                       ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc                >  2:320813/1‑150 (MQ=255)
                                                                                                                                       ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc                >  2:1539397/1‑150 (MQ=255)
                                                                                                                                       ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc                >  2:1657943/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)
                                                                                                                                                cTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCg                                                                                           >  1:1791931/1‑66 (MQ=255)
                                                                                                                                                cTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttcttc             >  1:2198579/1‑144 (MQ=255)
                                                                                                                                                cTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttcttc             <  2:2198579/144‑1 (MQ=255)
                                                                                                                                                 tCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGcc      >  2:235805/1‑150 (MQ=255)
                                                                                                                                                  cgcgGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGCCa     <  2:2396835/150‑1 (MQ=255)
                                                                                                                                                     gggCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGCCAGGa  >  2:1128614/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTCGTCAGCCAGGA  >  NC_002947/5165270‑5165568

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: