Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5690CG30.1% 149.3 / 34.3 83intergenic (+35/‑161)hrpB/mrcBputative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (0/25);  total (29/54)
Fisher's exact test for biased strand distribution p-value = 2.06e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTG  >  NC_000913/164421‑164710
                                                                                                                                                      |                                                                                                                                              
aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                >  1:450068/1‑149 (MQ=255)
aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                <  1:698027/149‑1 (MQ=255)
aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                >  2:1050182/1‑149 (MQ=255)
   ggAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                             <  2:273188/149‑1 (MQ=255)
    gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                            <  1:354821/149‑1 (MQ=255)
     aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                           >  1:1043256/1‑149 (MQ=255)
      gcttACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGTACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTATCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          <  2:1033978/146‑1 (MQ=255)
      gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          >  1:1513989/1‑149 (MQ=255)
      gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          <  1:813383/149‑1 (MQ=255)
      gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          >  2:28612/1‑149 (MQ=255)
       cGTACCGTGAGGTGCAAAAAGAGATGCAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAC‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                         <  2:1265746/149‑1 (MQ=255)
        gTACGGTGTGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:550717/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATCCTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  2:1975523/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:1505776/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  2:331198/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:213792/149‑1 (MQ=255)
        gTACCGTGAGGTGAAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGGA‑AAATACTGCACCGACGCGACGGACTAAAAAGTATTCTTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:550730/149‑1 (MQ=39)
         tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                       >  2:1842447/1‑149 (MQ=255)
           ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                     >  2:973531/1‑149 (MQ=255)
           ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGAAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                     >  1:1468256/1‑149 (MQ=255)
                    gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                            >  2:1413109/1‑149 (MQ=255)
                    gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGCCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                            >  2:2014494/1‑149 (MQ=255)
                      aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATTCTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                          >  2:320965/1‑149 (MQ=255)
                      aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                          >  1:1364871/1‑149 (MQ=255)
                        aacgacaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGAGGCGAGGGACTAAAAAGTATCCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        <  2:1249123/143‑1 (MQ=38)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                        >  2:669674/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:166398/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:921891/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:1306677/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCAAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:719629/1‑149 (MQ=255)
                         aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                       >  1:47066/1‑149 (MQ=255)
                          agagATCCAAGGGCGTTAACCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATAGTTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGt‑t                                                                                                                     <  2:2083344/149‑1 (MQ=38)
                            ataTGAATGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGA‑AAATACTGCACCGCCGCGCCGGACGAAAACGTTTTCTTAT‑GGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑tca                                                                                                                   <  2:166398/147‑1 (MQ=25)
                            agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGtgga                                                                                                                    >  2:1134151/1‑146 (MQ=255)
                               tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggcccc                                                                                                                 >  2:1785939/1‑143 (MQ=255)
                                gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tgacccg                                                                                                               >  2:1475562/1‑148 (MQ=255)
                                 aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tgacccac                                                                                                              >  1:2130270/1‑147 (MQ=255)
                                  atgggCTTTATCCAAAACTTGTCTGGCGGGACGCCCGGTC‑AAATACTCCACCGACGCGACGGACTAAAATTTATT‑GTAATGGTTGTTTTTCTCCCTCTCCCTGTGGTTTTGGGCGGGGTTGTGGGCATCAGCGCGCACGT‑tcaccctca                                                                                                             <  2:1865866/147‑1 (MQ=255)
                                    gggCGTTATCCCAAACATGTCTGGCCGGGCGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑gcacccgcacc                                                                                                           >  1:1837934/1‑149 (MQ=255)
                                        gTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACCCACCCCCaaa                                                                                                       >  2:1184534/1‑149 (MQ=38)
                                           aTCCCAAACATTTCGGGCCGGACG‑CCCGGCAAAATCCTGCACCGACGCGACGGCCGAAAACTTTTACGTAA‑GGTTGTTTTTCCCCCCCTCCGTGTGGGAGCGGGCGGGGTTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAAccc                                                                                                    <  2:1382279/149‑1 (MQ=255)
                                           aTCCCAAACATGTCTGGCGGGACGACCCGTA‑AATTACTCCACCGCGGCGCCTGACTAAAAATTTTTCTTTT‑TTTTGTTTTTCCCCCCCCCGGTGGGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAAccc                                                                                                    <  2:544305/149‑1 (MQ=255)
                                             cccAAACATGTCGGGCGGGACGACCCGGA‑AAATACTGCACCGACGCGACGGACTAAAATTTATTCTTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCctc                                                                                                  <  2:1937136/149‑1 (MQ=37)
                                               cAAACTTGTGGGGCGGGCCGACCCGCA‑AATTAGTGCACCGACGCGCGGGACTAAAAATTATTGGTAT‑GTTTGTTTTTCCCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCtctc                                                                                                <  2:355189/149‑1 (MQ=255)
                                               cAAACAGGTCGGGCCGCACCCCCCGAA‑AAAAACTCCCCCGACGCGACGGACTAAAAATTATTCTTAT‑GTTTGGTTTTCCCCCTCCCCGTGGGGTAGCGGGCGGGGGGGTGGCCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCtctc                                                                                                <  2:244395/149‑1 (MQ=255)
                                                aaccATGTCTGGCCGTACCCCCCGGC‑AAATACTGCACCGCCGCCCCGTACTAAAAATTATTCGTAA‑GTTGGGTTTTCGCCCACTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTcc                                                                                               <  2:942493/146‑1 (MQ=255)
                                                 aaCATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGT‑GGGCCCccacccccacccacaacc                                                                                              >  1:1862891/1‑143 (MQ=37)
                                                     ggcggggCGGGCCGACCCGCA‑AACTACGCCACCGACGCGCCGCACTAAAATTTATTCTTAC‑GGTTGGTTTCCCCCCTCCCCCTGGGGGCGGGGGCGGGGGTGGGGCCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                          <  1:2003351/144‑1 (MQ=255)
                                                       tCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGG‑TCACCCGCCCCCGGAGCCGCTCCCTCACgg                                                                                        >  2:1312122/1‑149 (MQ=37)
                                                        cTGGCCGGCCCCCCCGGC‑AAATACGGCCCCGCCGCGCCGGACTAAAAAGTATACGTAT‑GGTTGGGTTGCTCCCTCCCCCTGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAggg                                                                                       <  2:1084543/149‑1 (MQ=255)
                                                           gCCGGACGCCCCGGA‑AACTACTGCACCGAGGCGCCGGACTAAAAAGTATTCTAAC‑GGTTGGTTGACCCCCACCCCCGGGGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                    <  2:236371/149‑1 (MQ=255)
                                                           gCCGGACGCCCCGCA‑AAATACGGCCCCGCCGCGAGGCACTAAAAATTATTCGAAA‑GTTTGGTTTTCCCCCCCCCCCTGTGGGTGTGGGCGGGGTTGCGGCCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                    <  1:1506663/149‑1 (MQ=255)
                                                            ccGGACGACCCCGA‑AACTACTGCACCGACGCGCCGCACGAAAAAGTATCCGTAA‑GTTTGGTTTCCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                   <  2:1413724/149‑1 (MQ=37)
                                                                accacCCGGC‑AAATCCTCCCCCGACGCGACGCACGAAAAATTATTCTTAC‑GGTTGTTTTACCCCCCCACCGTGTGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                               <  1:1568941/146‑1 (MQ=255)
                                                                   acCCGCA‑AAATACTCCCCCCCCGCGCCGTACAAAAAAGTATTCTTAT‑GGTGGTTTTCCCCCCCCCCCCGGGGGGCGCGGGCCGGGGGGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                            <  2:562866/149‑1 (MQ=255)
                                                                     cccgcaaaaTACGCCACCGCCGCGCCGGCCTAAAATTTTTTGGTAT‑GGTTGTTTTTCCCCCTCTCCGTGGGGGTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                           <  1:1732796/143‑1 (MQ=255)
                                                                           aaaGACTGCACCGACGCGACGGACCAAAAAGTCTTCGTAA‑GGTTGGTTTTCTCCCACTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCACACCCTAACCCTCTACCTCAAGGGATATGGGAGCGATCGTGc                                                                     <  2:2056087/149‑1 (MQ=37)
                                                                              tACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGCCCCCCAACCAAACCCCCTCCCGCACGGGGGAGGGGGCCGCTCGAGcgca                                                                  >  1:2092573/1‑149 (MQ=25)
                                                                                    ccctccgcgacgTAATAAAATTTTTTCTTAT‑GGTTGGTTTTCCCCCTCCCCGTGGGGGTGTGGGCGGGGTTGTGGCCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                            <  2:638020/144‑1 (MQ=255)
                                                                                    ccccacgcgccgGACTAAAATGTTTTCGTAA‑GTTTGTTTTTCTCCCTCCCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                            <  2:526290/145‑1 (MQ=255)
                                                                                    aCCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGG‑GCCCCCCAACCCAAACCCACACCCGCAAGGGGGAGGGGCCCGCTCGCGCCCCAAtgtt                                                            >  2:2054880/1‑149 (MQ=21)
                                                                                    aCCGACGCGACGGAAGAAAATGTTTTCTTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                            <  2:1864433/149‑1 (MQ=255)
                                                                                     ccGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑GCCCCCTCCCCCGAACCCACTCCCGCGAGGGGGCGGGGCCCCAGCCAGCGCACAATTTg                                                           >  2:1358849/1‑149 (MQ=21)
                                                                                      cgccgcgaggGACTAAAAATTTTTGTTAT‑GGTTGGTTTTCCCCCTCCCCGGGGGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                          <  2:548003/146‑1 (MQ=255)
                                                                                          gcgacgGACTAAAAAGTTTTCGTTA‑GGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                      <  2:1653948/149‑1 (MQ=255)
                                                                                             acggactAAAATGTTTTCTTAT‑GTTTGGTTTTCGCCCACTCCGGGTGGGAGTGGGCCGGGGGGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                   <  2:1364871/149‑1 (MQ=255)
                                                                                                    aaattttATTCGTAT‑GGTTGGTTTTCACCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                            <  1:2014494/143‑1 (MQ=38)
                                                                                                    aaaaaGTTTTCTTAT‑GGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                            <  2:450068/149‑1 (MQ=38)
                                                                                                       aaGTATTGTTAT‑GGTTGTTTTTCTCCCTCCCCGTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGG‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                         <  2:1306677/149‑1 (MQ=255)
                                                                                                           aTTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                     <  2:1169680/149‑1 (MQ=255)
                                                                                                            atCGAAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGAGGGCATCAGCCCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                    <  1:2115819/148‑1 (MQ=255)
                                                                                                                  atggttggttTTCTCACACGCCCTGGGGCAGTGGGCGGGGGGGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTGTGAGAGATATCTTCTACTGTCTTGTAACag                               <  2:1406442/147‑1 (MQ=255)
                                                                                                                   tggttgtttTTCTCCCACTCCCTTTGGTAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                              <  2:112222/148‑1 (MQ=255)
                                                                                                                   tggttggttATCCCCCTCACCCTGTGGGATCGGGCCGGGGTGAGGCCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                              <  1:2085287/148‑1 (MQ=255)
                                                                                                                   agttttgttTTCCCCCTCGCCGTGGGGGTGGGGGCGGGGGTGTGGCCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                              <  1:1031121/143‑1 (MQ=255)
                                                                                                                    ggttgtttTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                             <  1:1511225/149‑1 (MQ=255)
                                                                                                                    ggttggttTACCCCCTCCCCGTGTGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                             <  2:1043256/149‑1 (MQ=255)
                                                                                                                           tttCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                      <  2:1779887/149‑1 (MQ=255)
                                                                                                                                ccctctccctGTGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                 <  1:1109585/149‑1 (MQ=255)
                                                                                                                                  ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACCCCCAACCAAAACCCCACCCGCAAGGGCGTGGGGTCCCCTCGAGCACCAATTTTGAGAGATATATTATTCAGTCTTTTTACAGAAAAAAAAAAAATAgg               >  2:1554526/1‑149 (MQ=16)
                                                                                                                                     tccctGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt            <  2:1989918/149‑1 (MQ=255)
                                                                                                                                       ccggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt          <  2:1468256/146‑1 (MQ=255)
                                                                                                                                         tgtgggtgtgggCGGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg        <  1:1475562/149‑1 (MQ=255)
                                                                                                                                             ggtgtgggCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc    <  2:1797774/144‑1 (MQ=255)
                                                                                                                                             ggAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc    <  2:213210/149‑1 (MQ=255)
                                                                                                                                               agagGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCCACCGCCACCC‑AACCCACACCCCCAAGGGAGGGGGGGCCCATCGGGCACCAATTTTGGGAGATAAATTGTTATGGTGTGTAAGCGAAGAACAGAAAAAAGGGCgttttggcccg  >  2:684559/1‑142 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
AAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTG  >  NC_000913/164421‑164710

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