Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5640AT30.4% 175.0 / 14.1 102intergenic (+30/‑166)hrpB/mrcBputative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base T (0/28);  total (33/69)
Fisher's exact test for biased strand distribution p-value = 2.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 113 total aligned reads displayed.

TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAA  >  NC_000913/164415‑164712
                                                                                                                                                       |                                                                                                                                                    
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      <  1:2438292/150‑1 (MQ=255)
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      <  2:2496131/150‑1 (MQ=255)
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      <  1:1962753/150‑1 (MQ=255)
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      <  1:398011/150‑1 (MQ=255)
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGgaga                                                                                                                                                      <  2:1982297/150‑1 (MQ=255)
 tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                                     <  2:357972/150‑1 (MQ=255)
 tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAACGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                                     <  2:1433502/150‑1 (MQ=255)
  cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                                    >  1:2179571/1‑150 (MQ=255)
   tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                                   >  1:70993/1‑150 (MQ=255)
   tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                                   <  1:85434/150‑1 (MQ=255)
    ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggga                                                                                                                                                  <  1:962653/150‑2 (MQ=255)
    ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                                  <  2:463342/150‑1 (MQ=255)
    ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                                  >  2:2146090/1‑150 (MQ=255)
     gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                 >  2:42459/1‑150 (MQ=255)
     gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                 >  2:169509/1‑150 (MQ=255)
     gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                 >  1:1090423/1‑150 (MQ=255)
      ataGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                <  1:2610995/148‑1 (MQ=255)
      aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                <  1:2610948/150‑1 (MQ=255)
       aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                               <  2:1737567/150‑1 (MQ=255)
       aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                               >  2:563238/1‑150 (MQ=255)
        aGGTGCGTACCGTTAGGTGAAAAAAGAGATGAAAGGGCGTTATCCCAATCTTTTCTGGCCGGACGACCCTGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCggg                                                                                                                                              <  1:482042/150‑1 (MQ=255)
        aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                              <  2:1598257/150‑1 (MQ=255)
        aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                              <  2:1098707/150‑1 (MQ=255)
         ggAGCGTACCGAGAGGTGCAAAAAGCGATGAACGGGCGGTATCCCAAACCTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                             <  2:1979399/150‑1 (MQ=255)
          gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                            >  2:288594/1‑150 (MQ=255)
          gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                            >  2:220272/1‑150 (MQ=255)
            gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                          <  1:1407445/150‑1 (MQ=255)
            gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                          <  2:1894531/150‑1 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        <  2:1998895/150‑1 (MQ=255)
               tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                       <  2:1929066/150‑1 (MQ=255)
                  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                    <  1:1057705/150‑1 (MQ=255)
                    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                  <  2:2720332/150‑1 (MQ=255)
                       ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                               >  1:622564/1‑150 (MQ=255)
                       ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                               >  1:622628/1‑150 (MQ=255)
                        gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                              >  2:1644934/1‑150 (MQ=255)
                        gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                              >  2:1644030/1‑150 (MQ=255)
                          gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                            >  2:839405/1‑150 (MQ=255)
                          gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                            >  1:1166002/1‑150 (MQ=255)
                           cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                           >  1:847580/1‑150 (MQ=255)
                            aaaaaGAGATTAATGGGCTTTATCCCAACATTGTCTGGCCGGACGACCCGGAAAATACTGAACCGACGGGACGGACTAAAAATTATTCTTATGGTTGGTTTTCGCCATATCCCGGTGGTAGTGGGCGGGGTTTTGGGCATCAGCGCGCAc                                                                                                                          <  2:2552139/150‑1 (MQ=255)
                            aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTAGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                          <  1:803029/150‑1 (MQ=255)
                             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGACAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGAGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                         >  2:1360360/1‑150 (MQ=255)
                              aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                        >  1:1129346/1‑150 (MQ=255)
                               aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGtg                                                                                                                       >  2:41325/1‑149 (MQ=255)
                                agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggc                                                                                                                      >  1:1457954/1‑147 (MQ=255)
                                    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttgtccc                                                                                                                  >  2:1841751/1‑145 (MQ=255)
                                     tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccg                                                                                                                 >  1:2451197/1‑143 (MQ=255)
                                     tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccc                                                                                                                 >  1:1972691/1‑143 (MQ=255)
                                        acgggCTTTAGCCAAAACATGTCGGGCCGCCCGACCGGGCAAATACTGCCCCGACGCGAGGGACTAAAAAGTTTTCGTAAGGTTGGTTTTCCCCCTCCCCCGGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                              <  2:1793015/148‑1 (MQ=25)
                                            gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccacaaccaa                                                                                                          >  1:1858371/1‑148 (MQ=255)
                                             ctttATCCAAAACATGTCTGGCCGGACCACCCGCAAAATACTGCACCGACGCGACGCACTAAAAATTTTTCTTATGTTTGGTTTTCTCCCACCCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                         <  2:874755/148‑1 (MQ=255)
                                                    ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCGCCCCCCAAcacccac                                                                                                  >  2:1968101/1‑144 (MQ=37)
                                                     cAACCATGTCGGGCCGGACGACCCGGAAATTACTGCCCCGACGCGAGGGACGAAAATTTATTCTTATGGTTGGTTTTCTCCCTCCCCGGGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                 <  2:1434115/150‑1 (MQ=255)
                                                     cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTGGCCCGCAACCCAAcacacacc                                                                                                 >  1:2435407/1‑143 (MQ=37)
                                                     aaaaCTTGTCGGGCCGGACGACCCGGAAATAACTGCACCTAGGCGCCGGAGTAAAAAGTTTTCTTAAGGTTGTTTTTCCCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                 <  2:1434105/149‑1 (MQ=255)
                                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCAAACCCACAccc                                                                                                >  1:2433430/1‑150 (MQ=37)
                                                       aaCTTGTGTGGCGGGACGCCCCGCAAAATCCGGCCCCGACGCTCCGTACTAAAATGTTTTGTTATGGTTGTTTTTCGCCCCCGCCGTGGGGGTGTGGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                               <  1:2182071/150‑1 (MQ=255)
                                                       aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCGAAACCACAccca                                                                                               >  2:2568444/1‑149 (MQ=37)
                                                           tttggggCCGGACTCCCCTGAAAATACTGCCCCGACGCCACGTACTAAAATTTATTCTTATGTTTGTTTTTCTCCACCGCCCTGTGGTTGTGGGCGGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                           <  2:48062/145‑1 (MQ=255)
                                                           gttCTGGCCGTACTACCCGGAAAATATGGCACCGACGCTACGTAGGAAAATTTTTTCGTACGGTTGTTTTTCTCCATCACCCTGTGGGTTTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                           <  2:2090411/148‑1 (MQ=255)
                                                            gTCGGGCCGGACGACCCGGCAAATCCTGCACCGCCGCGCCGGACGAAAACGTATACGTATGGTTGGTTTTCGCCCCCCCCCTGTGGGTGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                          <  2:1965762/150‑1 (MQ=255)
                                                                 gccggccgACCCGCCAAATACGGCACCGACGCGACGTACTAAAAAGTTTGCGTACGGTTGTTTTTCTCCCGCTCCCTGGGGGAGCGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                     <  2:1059190/150‑1 (MQ=255)
                                                                 gCCGTACGCCCCGCAAAATACGCCCCCGCCGCGACGGACTAAAAATTATCCTTATGGTTGGTTTTCTCCCACCCCCTGGGGGAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                     <  2:70993/150‑1 (MQ=255)
                                                                  ccGGCCGCCCCGCAAATTACTCCCCCGCCGCCACGCACTAAAATTTTTTCTAATGTTGGTTTTTCACCCCCCCCCTGGGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                    <  1:2568444/150‑1 (MQ=255)
                                                                  ccGGACGACCCGGCAAATCCGGCACCGACGCGACCGACGAAAAAGTATTCGCAAGGTTGGTTTTATCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                    <  2:1446060/150‑1 (MQ=255)
                                                                     gcccacCCGCAAAATACTGCCCCCCCGCGACGGACTAAAAATTATTCTTATGTTGGTTTTACCCCCACCCCGTGGGGGCGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                 <  1:2486083/146‑1 (MQ=255)
                                                                       cgccccGGCAATTCCTGCACCGACGCTACGGACTAAAAAGTTTTCTTAAGGTTGGTTTGCCCCATCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                               <  2:2482548/147‑5 (MQ=255)
                                                                       cgccccGGCAAATACGCCACCGCCGCGACGTACTAAAAAGTTTTGTTACGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                               <  2:2482569/147‑5 (MQ=25)
                                                                          cccGCCAAATCCTCCCCCGCCGCCCCGGCCCAAAAAGTATTCGTAAGGTTGGTTTTCCCCCTCGCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                            <  2:1726542/150‑1 (MQ=37)
                                                                           ccgtaaaaTACTCCACCGACGCTACGTACTAAAATTTTTTCTTATGGTTGTTTTACCCCCTCCCCCTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                           <  1:1671395/145‑1 (MQ=255)
                                                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGTTTTTCTCCCTCCCCCGGGGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                           <  2:1087100/150‑1 (MQ=39)
                                                                           ccGGAAAATACTGACCCCACGCCACGGCCTAAAAAGTATTCTAAGGGTTGGTTTTCCCCCGCTCCCTGGGGGCGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                           <  2:2386502/150‑1 (MQ=255)
                                                                             ggaaaaTACGCCCCCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCGCCCCCTCCCTGGGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                         <  2:2112108/147‑1 (MQ=255)
                                                                             gcaaaaTACTCCCCCGCCGCGACGGACGAAAAATTATTGGTATGGTTGTTTTTCACCCTCCCCCTGTGGGCGTGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                         <  1:2650151/147‑1 (MQ=255)
                                                                              taaaaTACTGCACCGAGGCGAGGGACTAAAATTTTTTGTTATGGTTGGTTTTCCCCCTCGCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                        <  1:1757154/148‑1 (MQ=255)
                                                                               aaaaTCCTCCCCCGCCGCCCCGGACCAAAAAGTATCCGAATGGTTGTTTTTCCCCCCCCCCCTGGGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                       <  2:847580/149‑1 (MQ=255)
                                                                                        ccccgccgCTACGCACAAAAATTTATCGTTAGGGTGGTTTTTCCCCCCCCCCGGGGGGGTGCGGGCGGGGGGGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                              <  1:81668/148‑1 (MQ=255)
                                                                                        ccccgccgCGCCGGACTAAAATGTTTTCTTATGTTGGTTTTTCACCGAACCCCTGGGGTAGTGGGGGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                              <  2:1193014/148‑1 (MQ=255)
                                                                                        cACCGCCGCGACGGACTAAAATGTTTTCGTATGGTTGGTTTTCTCCCACGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                              <  2:2179571/150‑1 (MQ=38)
                                                                                           cgacgcgccgGACTAAAAAGTTTTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                           <  2:2724096/150‑1 (MQ=38)
                                                                                               gcgacgGACTAAAAATTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                       <  1:2721698/150‑1 (MQ=38)
                                                                                               gccacgtactAAAAATTTTTCTTATGGTTGGTTTACGCCCTCCCCGGGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                       <  1:2471749/147‑1 (MQ=255)
                                                                                                 gccggactAAAAATTTTTCGTACGGTTGGTTTTCTCCCCCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                     <  2:2599962/148‑1 (MQ=38)
                                                                                                  ccggactAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                    <  2:1090423/149‑1 (MQ=38)
                                                                                                    ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTGGCCCCCACCCCAACCCCCACCCCCCCGGGGGCGGGGGCCGCTGGCGCCCAAAAATTTGGAGATAtatt                                                  >  1:2479045/1‑150 (MQ=16)
                                                                                                              ttATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                        <  2:1825242/149‑1 (MQ=255)
                                                                                                                tttctaaaGGTTGTTTTACCCCCGCCCCCGGTGGGAGCGGGCCGGGGTTAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                      <  2:554555/144‑1 (MQ=255)
                                                                                                                tttCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                      <  2:2451197/149‑1 (MQ=255)
                                                                                                                  tcttttggttgtttTTCGCCCCCTCCGTGTGGGTTTGGGCCGGGTTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                    <  1:41325/144‑1 (MQ=255)
                                                                                                                  tccttcggttgtttTCCGCCCCCCCCCTGTGGTCTCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                    <  1:1208422/144‑1 (MQ=255)
                                                                                                                  tCGTAAGTTTGGTTTTCTCCCACTCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                    <  2:870392/150‑1 (MQ=255)
                                                                                                                  tCGTAAGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                    <  2:1822037/150‑1 (MQ=255)
                                                                                                                   cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                   <  2:2581350/150‑1 (MQ=255)
                                                                                                                   cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                   <  2:2581552/150‑1 (MQ=255)
                                                                                                                           ttttttttCTCCATCTCCCTGTGGTTTTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                            <  1:563238/146‑1 (MQ=255)
                                                                                                                                       cctctcccgGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTTTGTCTTGTAACAGAAGAACAGAAAATCgg                 <  2:2128600/150‑1 (MQ=255)
                                                                                                                                        ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAAACCACTCCCTCAAGGGGGGGGGGGCCGAGAGAGCCCCAAATTTGTGGGTTATATTATTATGTGTTTTTTCAGAAGAACAGGAAATAggg                >  2:1496340/1‑150 (MQ=12)
                                                                                                                                             cgtgtgggtgtgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt           <  1:169509/148‑1 (MQ=255)
                                                                                                                                                   ggAGAGAGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt     <  1:1360360/150‑1 (MQ=255)
                                                                                                                                                      gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa  <  2:2134763/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAA  >  NC_000913/164415‑164712

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: