Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC34.1% 71.4 / 27.7 91intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base A (21/39);  new base C (3/28);  total (24/67)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164433‑164729
                                                                                                                                                       |                                                                                                                                                    
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     <  1:1057705/150‑1 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   <  2:2720332/150‑1 (MQ=255)
     ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                                >  1:622564/1‑150 (MQ=255)
     ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                                >  1:622628/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                               >  2:1644934/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                               >  2:1644030/1‑150 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                             >  1:1166002/1‑150 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                             >  2:839405/1‑150 (MQ=255)
         cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                            >  1:847580/1‑150 (MQ=255)
          aaaaaGAGATTAATGGGCTTTATCCCAACATTGTCTGGCCGGACGACCCGGAAAATACTGAACCGACGGGACGGACTAAAAATTATTCTTATGGTTGG‑TTTT‑CGCCATATCCC‑GGTGGTAGTGGGCGGGGTTTTGGGCATCAGCGCGCAc                                                                                                                                           <  2:2552139/150‑1 (MQ=255)
          aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTAGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           <  1:803029/150‑1 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGACAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGAGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                          >  2:1360360/1‑150 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                         >  1:1129346/1‑150 (MQ=255)
             aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGtg                                                                                                                                        >  2:41325/1‑149 (MQ=255)
              agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggc                                                                                                                                       >  1:1457954/1‑147 (MQ=255)
                 tattaagggggGTTATCCCAAAAATGTGGGGCGGTACGACCCGGAAAATACTGCACGGACGCGACGGACTAAAATTTTTTCTTATGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                    <  2:2124625/143‑1 (MQ=21)
                  aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttgtccc                                                                                                                                   >  2:1841751/1‑145 (MQ=255)
                   tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccg                                                                                                                                  >  1:2451197/1‑143 (MQ=255)
                   tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccc                                                                                                                                  >  1:1972691/1‑143 (MQ=255)
                      acggggtttATCCCAAACTGGTCGGGCCGGACGACCCGCAAAATACTCCACCTCCGCTACGTACTAAAATTTATTCTTATGTTTGT‑TTTT‑CTCCCCCTCCC‑TGGGGGTGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                               <  2:888418/143‑1 (MQ=255)
                      acgggCTTTAGCCAAAACATGTCGGGCCGCCCGACCGGGCAAATACTGCCCCGACGCGAGGGACTAAAAAGTTTTCGTAAGGTTGG‑TTTT‑CCCCCTCCCCC‑GGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                               <  2:1793015/148‑1 (MQ=25)
                      aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                               >  1:874755/1‑150 (MQ=255)
                        gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccgccccc                                                                                                                             >  2:1913935/1‑150 (MQ=255)
                          gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccacaaccaa                                                                                                                           >  1:1858371/1‑148 (MQ=255)
                           ctttATCCAAAACATGTCTGGCCGGACCACCCGCAAAATACTGCACCGACGCGACGCACTAAAAATTTTTCTTATGTTTGG‑TTTT‑CTCCCACCCCC‑TGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                          <  2:874755/148‑1 (MQ=255)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCGCCCCCCAAcacccac                                                                                                                   >  2:1968101/1‑144 (MQ=37)
                                   cAACCATGTCGGGCCGGACGACCCGGAAATTACTGCCCCGACGCGAGGGACGAAAATTTATTCTTATGGTTGG‑TTTT‑CTCCCTCCCCG‑GGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                  <  2:1434115/150‑1 (MQ=255)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTGGCCCGCAACCCAAcacacacc                                                                                                                  >  1:2435407/1‑143 (MQ=37)
                                   aaaaCTTGTCGGGCCGGACGACCCGGAAATAACTGCACCTAGGCGCCGGAGTAAAAAGTTTTCTTAAGGTTGT‑TTTT‑CCCCCTCTCCC‑TGTGGTAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                  <  2:1434105/149‑1 (MQ=255)
                                    aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCAAACCCACAccc                                                                                                                 >  1:2433430/1‑150 (MQ=37)
                                     aaCTTGTGTGGCGGGACGCCCCGCAAAATCCGGCCCCGACGCTCCGTACTAAAATGTTTTGTTATGGTTGT‑TTTT‑CGCCCCCGCCG‑TGGGGGTGTGGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                <  1:2182071/150‑1 (MQ=255)
                                     aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCGAAACCACAccca                                                                                                                >  2:2568444/1‑149 (MQ=37)
                                        aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCGCCGCCCCCAAACCCACTcccacca                                                                                                             >  2:2619317/1‑146 (MQ=37)
                                         tttggggCCGGACTCCCCTGAAAATACTGCCCCGACGCCACGTACTAAAATTTATTCTTATGTTTGT‑TTTT‑CTCCACCGCCC‑TGTGGTTGTGGGCGGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:48062/145‑1 (MQ=255)
                                         tgacgggCGGGACCACCCGCAAAATCCTCCCCCGCCGCGAGGGACTAAAATGTATTCGAATGGTTGG‑TTTA‑CCCCCACGCCC‑GGGGGGAGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:48111/145‑1 (MQ=255)
                                         gttCTGGCCGTACTACCCGGAAAATATGGCACCGACGCTACGTAGGAAAATTTTTTCGTACGGTTGT‑TTTT‑CTCCATCACCC‑TGTGGGTTTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:2090411/148‑1 (MQ=255)
                                          gTCGGGCCGGACGACCCGGCAAATCCTGCACCGCCGCGCCGGACGAAAACGTATACGTATGGTTGG‑TTTT‑CGCCCCCCCCC‑TGTGGGTGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:1965762/150‑1 (MQ=255)
                                               gccggccgACCCGCCAAATACGGCACCGACGCGACGTACTAAAAAGTTTGCGTACGGTTGT‑TTTT‑CTCCCGCTCCC‑TGGGGGAGCGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  2:1059190/150‑1 (MQ=255)
                                               gCCGTACGCCCCGCAAAATACGCCCCCGCCGCGACGGACTAAAAATTATCCTTATGGTTGG‑TTTT‑CTCCCACCCCC‑TGGGGGAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                      <  2:70993/150‑1 (MQ=255)
                                                ccGGCCGCCCCGCAAATTACTCCCCCGCCGCCACGCACTAAAATTTTTTC‑TAATGTTGGTTTTT‑CACCCCCCCCC‑TGGGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                     <  1:2568444/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATCCGGCACCGACGCGACCGACGAAAAAGTATTCGCAAGGTTGG‑TTTT‑ATCCCTCTCCC‑TGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                     <  2:1446060/150‑1 (MQ=255)
                                                   gcccacCCGCAAAATACTGCCCCCCCGCGACGGACTAAAAATTATTC‑TTATGTTGG‑TTTTACCCCCACCCCG‑TGGGGGCGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  1:2486083/146‑1 (MQ=255)
                                                     cgccccGGCAATTCCTGCACCGACGCTACGGACTAAAAAGTTTTCTTAAGGTTGG‑TTTG‑CCCCATCTCCC‑TGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                                                <  2:2482548/147‑5 (MQ=255)
                                                     cgccccGGCAAATACGCCACCGCCGCGACGTACTAAAAAGTTTTGTTACGGTTGT‑TTTT‑CTCCCTCTCCC‑TGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                                                <  2:2482569/147‑5 (MQ=25)
                                                        cccGCCAAATCCTCCCCCGCCGCCCCGGCCCAAAAAGTATTCGTAAGGTTGG‑TTTT‑CCCCCTCGCCC‑GGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                             <  2:1726542/150‑1 (MQ=37)
                                                         ccgtaaaaTACTCCACCGACGCTACGTACTAAAATTTTTTCTTATGGTTGT‑TTTA‑CCCCCTCCCCC‑TGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  1:1671395/145‑1 (MQ=255)
                                                         ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGT‑TTTT‑CTCCCTCCCCC‑GGGGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:1087100/150‑1 (MQ=39)
                                                         ccGGAAAATACTGACCCCACGCCACGGCCTAAAAAGTATTCTAAGGGTTGG‑TTTT‑CCCCCGCTCCC‑TGGGGGCGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:2386502/150‑1 (MQ=255)
                                                           ggaaaaTACGCCCCCGACGCGACGGACGAAAAAGTATTCGTATGGTTGG‑TTTT‑CGCCCCCTCCC‑TGGGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:2112108/147‑1 (MQ=255)
                                                           gcaaaaTACTGCCCCGACGCGACGGACTAAAATGTATTCGTATGGTTGG‑TTTT‑CCCCCTCTCCC‑TGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:1664797/147‑1 (MQ=255)
                                                           gcaaaaTACTCCCCCGCCGCGACGGACGAAAAATTATTGGTATGGTTGT‑TTTT‑CACCCTCCCCC‑TGTGGGCGTGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  1:2650151/147‑1 (MQ=255)
                                                            taaaaTCCTCCACCGACGCCTCGGACAAAATTTTATTCTTATGTTGTG‑TTTT‑CGCCCGCCCCC‑TGGGGGCGGGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  1:1503584/148‑1 (MQ=255)
                                                            taaaaTACTGCACCGAGGCGAGGGACTAAAATTTTTTGTTATGGTTGG‑TTTT‑CCCCCTCGCCG‑TGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  1:1757154/148‑1 (MQ=255)
                                                             aaaaTCCTCCCCCGCCGCCCCGGACCAAAAAGTATCCGAATGGTTGT‑TTTT‑CCCCCCCCCCC‑TGGGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                        <  2:847580/149‑1 (MQ=255)
                                                                      ccccgccgCTACGCACAAAAATTTA‑TCGTTAGGGTGGTTTTT‑CCCCCCCCCCG‑GGGGGGTGCGGGCGGGGGGGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:81668/148‑1 (MQ=255)
                                                                      ccccgccgCGCCGGACTAAAATGTTTTCTTATGTTGGT‑TTTT‑CACCGAACCCC‑TGGGGTAGTGGGGGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  2:1193014/148‑1 (MQ=255)
                                                                      cACCGCCGCGACGGACTAAAATGTTTTCGTATGGTTGG‑TTTT‑CTCCCACGCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  2:2179571/150‑1 (MQ=38)
                                                                        ccgccgCGACGGACGAAAAAGTATTCTTAAGGTGGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:2466183/150‑1 (MQ=255)
                                                                         cgacgcgccgGACTAAAAAGTTTTCGTATGGTTGG‑TTTT‑CCCCCTCTCCC‑TGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                            <  2:2724096/150‑1 (MQ=38)
                                                                             gcgacgGACTAAAAATTATTCGTATGGTTGT‑TTTT‑CTCCCTCTCCC‑TGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                        <  1:2721698/150‑1 (MQ=38)
                                                                             gccacgtactAAAAATTTTTCTTATGGTTGG‑TTTA‑CGCCCTCCCCG‑GGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                        <  1:2471749/147‑1 (MQ=255)
                                                                               gccggactAAAAATTTTTCGTACGGTTGG‑TTTT‑CTCCCCCGCCC‑TGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                      <  2:2599962/148‑1 (MQ=38)
                                                                                ccggactAAAAAGTTTTCGTATGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                     <  2:1090423/149‑1 (MQ=38)
                                                                                  ggacgAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTGGCCCCCACCCCAACCCCCACCCCCCCGGGGGCGGGGGCCGCTGGCGCCCAAAAATTTGGAGATAtatt                                                                   >  1:2479045/1‑150 (MQ=16)
                                                                                            ttATTCGTATGGTTGG‑TTTT‑CTCCCTCTCCG‑TGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:1825242/149‑1 (MQ=255)
                                                                                              tttctaaaGGTTGT‑TTTA‑CCCCCGCCCCC‑GGTGGGAGCGGGCCGGGGTTAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                       <  2:554555/144‑1 (MQ=255)
                                                                                              tttCGTATGGTTGG‑TTTT‑CTCCCTCTCCG‑TGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                       <  2:2451197/149‑1 (MQ=255)
                                                                                                tcttttggttgt‑ttTT‑CGCCCCCTCCG‑TGTGGGTTTGGGCCGGGTTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  1:41325/144‑1 (MQ=255)
                                                                                                tccttcggttgt‑ttTC‑CGCCCCCCCCC‑TGTGGTCTCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  1:1208422/144‑1 (MQ=255)
                                                                                                tCGTAAGTTTGG‑TTTT‑CTCCCACTCCC‑TGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  2:870392/150‑1 (MQ=255)
                                                                                                tCGTAAGGTTGG‑TTTT‑CGCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  2:1822037/150‑1 (MQ=255)
                                                                                                 cGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  2:2581350/150‑1 (MQ=255)
                                                                                                 cGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  2:2581552/150‑1 (MQ=255)
                                                                                                     ttgttgg‑ttTT‑CCCCCCCCCCCGGGGGGGCGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  1:220272/148‑1 (MQ=255)
                                                                                                         tttttttt‑CTCCATCTCCC‑TGTGGTTTTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                             <  1:563238/146‑1 (MQ=255)
                                                                                                                     cctctccc‑gGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTTTGTCTTGTAACAGAAGAACAGAAAATCgg                                  <  2:2128600/150‑1 (MQ=255)
                                                                                                                      ctctccc‑tGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAAACCACTCCCTCAAGGGGGGGGGGGCCGAGAGAGCCCCAAATTTGTGGGTTATATTATTATGTGTTTTTTCAGAAGAACAGGAAATAggg                                 >  2:1496340/1‑150 (MQ=12)
                                                                                                                            cgtgtgggtgtgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                            <  1:169509/148‑1 (MQ=255)
                                                                                                                                  ggAGAGAGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      <  1:1360360/150‑1 (MQ=255)
                                                                                                                                     gggggCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                   <  1:703025/148‑1 (MQ=255)
                                                                                                                                     gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                   <  2:2134763/150‑1 (MQ=255)
                                                                                                                                         gCCGGGGTGCGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  1:2146090/150‑1 (MQ=255)
                                                                                                                                           cGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:847278/150‑1 (MQ=255)
                                                                                                                                            ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:2498431/150‑1 (MQ=255)
                                                                                                                                            ggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:1139936/150‑1 (MQ=255)
                                                                                                                                            ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:1457954/150‑1 (MQ=255)
                                                                                                                                              ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  2:1972691/150‑1 (MQ=255)
                                                                                                                                                    ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:52240/150‑1 (MQ=255)
                                                                                                                                                    ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:1248755/150‑1 (MQ=255)
                                                                                                                                                     ccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg   <  1:823464/149‑1 (MQ=255)
                                                                                                                                                      cATCAGCGCGCACGATCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  2:2361018/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164433‑164729

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: ‑