Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 19.9% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC19.9% 136.5 / 13.6 87intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base C (2/17);  total (35/52)
Fisher's exact test for biased strand distribution p-value = 3.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728
                                                                                                                                                     |                                                                                                                                                    
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  1:369194/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGTCGGGGTGAGGGCAt                                                                                                                                                     >  2:473573/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCag                                                                                                                                                     >  1:2581194/1‑149 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  2:1040223/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     <  2:697964/150‑1 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCAGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                               >  2:3920518/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgcg                                                                                                                                              >  1:2005029/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           >  1:3336334/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           >  2:3537717/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           >  1:1367068/1‑150 (MQ=255)
                gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                      >  2:1605340/1‑150 (MQ=255)
                      aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccccc                                                                                                                                >  1:3419221/1‑147 (MQ=255)
                        aGGGCGTTATCCCAAACATGTCGGGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                              <  2:2625144/150‑1 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCCAAc                                                                                                                         >  2:184491/1‑150 (MQ=38)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCaaaa                                                                                                                         >  2:3217745/1‑149 (MQ=255)
                             gTTATCCCAAACATGTCGGCCCGCACGACCCGGCAAATACGGCACCGACGCGACGCACGAAAAATTTTGCGTACGGTTGGTTTTCCCCCTCCCCCTGTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                         <  2:2143591/150‑1 (MQ=255)
                             gTTATCCCAAAAAGGTCTGGCCGGACGACCCGCCAAATACTGCCCCGACGCGACGCACTAAAAAGTATTCGTATGGTTGGTTTTCACCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                         <  2:2143610/150‑1 (MQ=37)
                             gTTAGCCCAAACCGGCCTGGCCGGACCACCCGGCAAATCCTGCACCGACGCGACGGACGAAAATTTATCCGAAAGGTTGTTTTTCTCCCTCTCCCCGTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                         <  2:2143623/150‑1 (MQ=255)
                                  cccAAACCTGTGGGGCCGGACCCCCCGGAAAACCCGCCCCCGCCGCGCCGGACAAAAATGTTTTGTTTTGTTGGGTTTACCCCCCCGCCGGGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                    <  1:1537820/150‑1 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAAaccaca                                                                                                                    >  2:2202277/1‑147 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGGCCGCACCCCAAccccca                                                                                                                    >  1:726945/1‑147 (MQ=37)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAAaccaca                                                                                                                    >  1:574477/1‑147 (MQ=38)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGTGGACCCGCAACCCAAaccaca                                                                                                                    >  1:574526/1‑147 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCAAAcccacac                                                                                                                   >  2:1999880/1‑146 (MQ=255)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAACCCACAccc                                                                                                                 >  1:534611/1‑150 (MQ=37)
                                         ctggcgggccggccGCCCGGGCAAATCCGGCCCCGCCGCGACGGAGGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGCGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:3383516/144‑1 (MQ=25)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:2420494/150‑1 (MQ=255)
                                         aTGTCGGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAATTTATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:2003609/150‑1 (MQ=38)
                                           gTCGGGCCGCACGACCCGGAAAATACGGCCCCGACGCGAGGGACTAAAATTTTTTCGTACGGTTGTTTTTCGCCCACGCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:392282/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCGACCCCACCCCAACCCTAACCCACACgggcggg                                                                                                      >  2:2057916/1‑146 (MQ=37)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  2:3905575/150‑1 (MQ=255)
                                                 ccGGACCACCCGCAAAATACTCCACCGAGGCGCCTGACGAAAAATTATTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  1:868897/150‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGTCCCACACCCAAAACCACACCCGCAAAGGAGAgggg                                                                                                   >  1:3138580/1‑150 (MQ=37)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCAAAACCCCACCCCCAAaggggagggg                                                                                                   >  2:1809060/1‑150 (MQ=37)
                                                    gacgccCCGCAAAAGACCGCCCCCACGCGCCGCCCGAAAACGTATTCGTACGGTTGTTTTTCCCCCGCCCCGGGGGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  2:3660711/150‑1 (MQ=255)
                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAACCACACCCCCAAggggggggggc                                                                                                  >  2:3582522/1‑149 (MQ=37)
                                                    ccccacCCGCAAAATACTCCCCCGCCGCGCCGGACCAAAATGTTTTCTTATGGTTGTTTTTCCCCCCCCCCCGGGGGTTGCGGGCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  2:3336334/146‑1 (MQ=255)
                                                         cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACACCCCAACCCACACCCACAACGGGGAGGGGGccccg                                                                                             >  1:3347825/1‑147 (MQ=25)
                                                               aaagccTGCCCCGACGCGAGGGCCGAAACCGTTGGCGTCGGGTTGGTTTTCGCCCTCCCCCTGGGGGGGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                      <  2:52656/145‑1 (MQ=255)
                                                               aaaTACTGCACCGACGCGACGGACGAAAATGTTTTCGTTTGTTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGGGGGGGTGAGGTCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                      <  2:52697/150‑1 (MQ=37)
                                                                         ccGACGCGCGGGACTAAAATGTATTCTTATGGTGGGTTTTCGCCCCCTCCGTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  2:3814320/150‑1 (MQ=255)
                                                                         ccGACGCGCCGCACGAAAATTTTTTCGTAGGGTTGGTTTCCCCCCTCGCCGTGGGGGAGGGGGCGGGGGTGGGGCCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  1:2029620/150‑1 (MQ=255)
                                                                          cgacgcgccgAACTAAAATTTTTTGGTTTGTTTGGTTTACCCCCCCCCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                           <  1:1076144/150‑1 (MQ=255)
                                                                          cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCACAACCACACCCACCAAGGCGACGGGGCCCATCGAGCACaaaaaaaatg                                                                           >  1:3029145/1‑143 (MQ=21)
                                                                          cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCGCGCCCACAAGGGAGAGGGGGCCCAGCGAGCCCAAAATTTgag                                                                           >  1:1624718/1‑150 (MQ=25)
                                                                          cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCCAACCCCCGACCCCAACGGCGGGGGGCCCGGTGGTCCACAAATATTtag                                                                           >  1:1917904/1‑150 (MQ=21)
                                                                              gcgacgGATTAAAATGTTTGCGTATGGTTGTTTTGCCCCCTCGCCGTGTGGTAGTGGGCGGGGTTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                       <  2:3686796/150‑1 (MQ=255)
                                                                               cgccgtacgAAAATTTTTTGGTATGGTTGTTTTGCCCCCTCTCCCTGGGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                      <  2:3422178/144‑1 (MQ=255)
                                                                                  cggacgAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                   <  2:3241967/150‑1 (MQ=255)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCAAACCCAATCCCGCAACGGCGAGGGGACCCATCCAGCCCAAAATTTTAGAGaaataat                                                                  >  2:2874189/1‑147 (MQ=21)
                                                                                              ttttCGTATGGTTGTTTTTCGCCCTCCCCGGGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  2:1106477/148‑1 (MQ=255)
                                                                                              tATTCGTAAGGTTGGTTTTCCCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  1:3464733/150‑1 (MQ=255)
                                                                                                 tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                    <  2:2095910/150‑1 (MQ=255)
                                                                                                   gtttgggtggttTTCCCCCTCTCCCTGTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                  <  1:255338/143‑1 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAGGGAGAGGGGACCCAACCAGCACAAAATTTTAGAGATATCTTCTTCTGTCTTGTaaaac                                                 >  2:3885063/1‑147 (MQ=25)
                                                                                                     atggttggttTTCGCCCACGCCCTGTGGTAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                <  1:2353412/148‑1 (MQ=255)
                                                                                                     aaGGTTGGTTTTCCCCCACACCCTGTGGGATAGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                <  1:3385089/150‑1 (MQ=255)
                                                                                                      tgtttggttTTCTCCCCCGCCGTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                               <  1:1216689/147‑1 (MQ=255)
                                                                                                      cggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                               <  1:2965978/149‑1 (MQ=255)
                                                                                                       ggttggttTTCTCCATCTGACGGTGGGAGTGTGCGGGGGTTCGGGCCTCATCGTGCACGTTCACCCTCACGCTAACGGTCTCGCTGTTGGGTGTTGGGTCGGATCGAGCACTATTTTTGTTATATATCTTCTTCTGTCTTGTAACagtag                                              <  2:2605181/150‑1 (MQ=255)
                                                                                                        gttggttCTCCCCCTCCCCGTGGGGGTGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                             <  1:1112359/150‑1 (MQ=255)
                                                                                                         ttgtttTCCGCCATCTCGGTGTGGTTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  2:3638093/150‑1 (MQ=255)
                                                                                                         ttggttGCCCCCCTCGCCCGGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  2:3382571/150‑1 (MQ=255)
                                                                                                                accccctctccctGTGGTCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                     <  2:280916/147‑1 (MQ=255)
                                                                                                                 ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCGAAACCTCTCCCACAAGGGAGAGGGGGCCCCGCGCGCACAAAATTTGAGAGATATCATCTTATGTCTTGTAACAGAAGAACACaaaaa                                    >  2:2921547/1‑149 (MQ=25)
                                                                                                                 ctccctcgccctGGGGTTGGGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                    <  1:489889/150‑1 (MQ=255)
                                                                                                                     ctctcccgGTGGGAGGGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                <  1:1093371/150‑1 (MQ=255)
                                                                                                                          cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAAACCACAACCAAAAGGGGGTGGGGGCCCCTCGCGCACAAAATTTTAGAGAAAAATTCTTATGTGTTGTTAAATAAGAACACAAAAAAGGGCgttt                           >  2:3549957/1‑150 (MQ=12)
                                                                                                                          ccacagggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCACACCCGCAAGGGGGAGGGGGCCGGTCGAGCACAAAATTTGATGGCTATGTTATTAAGTGGTTTAACAGAATAACAGAAAAACGGGCgttt                           >  2:2429892/6‑150 (MQ=255)
                                                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCTAAGCCTCTGCCTCAAGGGTGAGGGGGCCGATCCAGGACAAAATTTGTGATTTATCTTCTTCTGTGTTGGGACGGACGAACCGAAAATggggggttttg                         >  2:2893174/1‑143 (MQ=19)
                                                                                                                              tgGTAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  1:3140986/150‑1 (MQ=255)
                                                                                                                              tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  2:2956589/150‑1 (MQ=255)
                                                                                                                              ggggAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  1:3140604/149‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                     <  1:2466451/150‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                     <  2:292784/150‑1 (MQ=255)
                                                                                                                                gctgtgggCGGGGTTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                     <  1:334852/145‑1 (MQ=255)
                                                                                                                                    agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                 <  2:261392/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAAACCACACACACAAAGGGGAGGGGGCCCCTAGAGCCCCAAATTTTGGATATATATTCATCAGTCTTGTAAACCAAGAACAGAAAAACCGGCGTTTGCGccgggctat               >  2:1266197/1‑143 (MQ=255)
                                                                                                                                        ccGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  1:3755616/150‑1 (MQ=255)
                                                                                                                                         gggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  2:1796223/149‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  2:1796233/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:305702/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTGCAACCACACCCAAAACCACAACCGCAAAGGAGAGGGGCCCCCTCGAGCCCAAATTTTGAGTGATATATTCTAAATTCTTGTGAAACAAGAACAGAAAATAGGGCTTTTGCGCCTGAAAAtttgggg         >  2:2408090/1‑145 (MQ=255)
                                                                                                                                               gAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa      <  2:1531354/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  2:574477/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCAACCACAACCCAAACCAAACCCACAAGGGAGAGGGGACCCAAACAGCACCAAATTTTATAAATAAAATATTAAGTCATTTAAAAGAAGAAAAGAAAAAAGGGCATTTTCGCCGGAATATTGCGGCGaaaaa   >  2:2523939/1‑150 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTAACCCACACCCGAAACCACTCCCACAAGGGAGGGGGGACCGATCGAGCACAAAATTTTAGAGATATAATATTCAGTCTTTTTAAAGAAGAAAAGAAAAACGGGCTTTTTAGCCTGAATATTGGGGGGaaaaa   >  2:970448/1‑150 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  2:574526/150‑1 (MQ=255)
                                                                                                                                                   gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg  <  1:463976/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728

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