Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5790GC22.4% 130.9 / 12.1 58intergenic (+45/‑151)hrpB/mrcBputative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base C (0/13);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 8.81e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGG  >  NC_000913/164431‑164720
                                                                                                                                                      |                                                                                                                                              
aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                               >  1:2143259/1‑150 (MQ=255)
aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                               >  1:2142941/1‑150 (MQ=255)
 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                              >  1:669761/1‑150 (MQ=255)
 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGCCAAATAATG‑CCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                              <  2:718291/150‑1 (MQ=255)
    tGAGGTGCAAAAAGAGATGAATGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                           <  2:521843/150‑1 (MQ=255)
    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                           <  2:2051748/150‑1 (MQ=255)
    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                           <  2:955095/150‑1 (MQ=255)
        gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                       >  1:426749/1‑150 (MQ=255)
           cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                    >  1:590371/1‑150 (MQ=255)
           cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                    >  2:161922/1‑150 (MQ=255)
           cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                    >  2:161887/1‑150 (MQ=255)
           cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                    >  1:113416/1‑150 (MQ=255)
            aaaaaGAGATGAAAGTGAGTTATCCCAAACATGACTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGTGCATAAGCGCGCAc                                                                                                                                   >  1:1410752/1‑150 (MQ=255)
             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCg                                                                                                                                  >  1:310426/1‑150 (MQ=255)
             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                  >  1:2074116/1‑150 (MQ=255)
             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                  >  1:2074121/1‑150 (MQ=255)
             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                  >  1:1462527/1‑150 (MQ=255)
             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                  >  1:310586/1‑150 (MQ=255)
                          gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgggcccacaccc                                                                                                                     >  1:588338/1‑150 (MQ=255)
                               ttATCCCAACTATGTCTGGCCGGACCACCCGGAAAATACTC‑CCCCGCGGGGACGGACCAAAATTTATTCTTATGTTTGTTTTT‑CTCCCTCTCCCTTGGGGTGTGGGCGGGGTTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                <  2:47374/150‑1 (MQ=255)
                                tATCCCAAACTTGTGGGGCCGCACGACCCGTAAAAAACTG‑CCCCGACGCCACGAACAAAAATTTATTCGTATGGTGGTTTTT‑CCCCCCCCCCGGGGGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                               <  2:252678/150‑1 (MQ=255)
                                    ccAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCAAAaccacgc                                                                                                           >  2:405775/1‑146 (MQ=37)
                                    ccAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCCCCCCAAcccacac                                                                                                           >  2:1613061/1‑146 (MQ=37)
                                     aaaaCTTGTGGGGCCGTACGACCCGGAAAATACTG‑CACCGCCGGTACGGCCTAAATTTTTTTCTTATGTTTGTTTTT‑CTCCCTCTCCGTGTGGGTGGGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                          <  1:1567787/149‑1 (MQ=255)
                                                   cgcaccccccGGAAAAAACTG‑CCCCGCCGCGCCGGACTAAAAAGTATTCGTACGGTTGGTTTC‑CACCCCCCCCCTGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAggg                                                                                           <  2:1480563/143‑1 (MQ=255)
                                                   ccgcacccccGGCAAAAACCGCCCCCCACGCGACGCACAAAAATTTTTGCTTATGTTTGGTTTC‑CCCCCCCCCCCGGGGGGTGGGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAgg                                                                                            <  1:628169/143‑1 (MQ=255)
                                                   cGGACGACCCGGCAAATACTG‑CACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCCCCGTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAggg                                                                                           <  2:1176798/150‑1 (MQ=255)
                                                   cGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TTGAGGGAGAggg                                                                                           <  2:2142941/150‑1 (MQ=255)
                                                   cGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TTGAGGGAGAggg                                                                                           <  2:2143259/150‑1 (MQ=255)
                                                      acgccCCGGAAAACACTC‑CACCGCCGCGACGGACTAAAAAGTATTCGTACGGTTGTTTTA‑CCCCCTCCCCCTGGGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGAc                                                                                        <  2:166772/150‑1 (MQ=255)
                                                        gacCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAAACCACACCCACCAAGGGGGGGGGGccc                                                                                       >  2:376195/1‑150 (MQ=25)
                                                            cGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCAACCGCACCCGAACCCTCAACC‑GCCTGGGTGGGGGGACCGATCg                                                                                  >  2:1557834/1‑150 (MQ=25)
                                                                          aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                     <  2:1304957/150‑1 (MQ=255)
                                                                           ccgccgCCAGGGACTAAAATGTTTTCGTATGGTTGGTTTT‑CTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                    <  1:2105829/150‑1 (MQ=255)
                                                                             gacgcgccgGACTAAAAATTATTCGTATGTGTGTTTTT‑CACCCACTCCCTTGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                  <  2:585944/150‑1 (MQ=255)
                                                                                  gacggacgAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAACCACAACC‑CCCAGGGGGAGGGGGCCGCGAGAGCACCAAATTTTAGTGatat                                                             >  2:269281/1‑150 (MQ=18)
                                                                                            aatgTAATCGAACGGTGGGTTTT‑CTCCCTCACCGTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                   <  2:1086386/147‑1 (MQ=255)
                                                                                            aaaGTATCCGTAAGGTTGGTTTT‑CTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                   <  2:982381/150‑1 (MQ=255)
                                                                                                 atgcGTATGGTTGTTTTT‑CCCCCTCCCCCTGTGGTTGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                              <  1:1934653/147‑1 (MQ=255)
                                                                                                    cgaaggtttggttTC‑CCCCCCCCCCGTGGGGTTGGGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                           <  1:1020039/143‑1 (MQ=255)
                                                                                                     gTAACGGTGGTTTTACTCCCACACCCTGGGGGTGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                           <  1:1020049/150‑1 (MQ=255)
                                                                                                       aaGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                        <  2:2069167/150‑1 (MQ=255)
                                                                                                                 tt‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGTCCCACACCCAAAACCACTCCC‑GCAACGGGGCGGGGGCCCATCGTGCACAAAATTTTTGGGTTATATTATTATTTGTTGTTACAGAagaaaagaaa                              >  2:1150242/1‑150 (MQ=12)
                                                                                                                    ctccctcgccctGTGGTTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                            <  2:355570/150‑1 (MQ=255)
                                                                                                                       cctccccgtGTGGGTGTGGGGGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                         <  1:460475/150‑1 (MQ=255)
                                                                                                                            ccctGGGGGAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                    <  2:1366110/150‑1 (MQ=255)
                                                                                                                               tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCCCCCAAAACCACTCCC‑GCAAGGGGGGGGGGTCCGATAGAGCACCAATTTTTTGAGTTTTATTATTATGTGTTTTTAAAGAATAACAGAAAAAAGGGGtttttg                 >  2:2062183/1‑148 (MQ=255)
                                                                                                                                   ggcgtgggCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt             <  1:1996943/145‑1 (MQ=255)
                                                                                                                                     tgagGGCCGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa           <  2:727117/149‑1 (MQ=255)
                                                                                                                                       tgggCGGGGTTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat         <  2:806249/149‑1 (MQ=255)
                                                                                                                                       tgggCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat         <  2:805941/149‑1 (MQ=255)
                                                                                                                                       tgggCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat         <  2:806183/149‑1 (MQ=255)
                                                                                                                                          gcgggggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt      <  1:1263199/147‑1 (MQ=255)
                                                                                                                                           gcGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg     <  1:1411305/149‑1 (MQ=255)
                                                                                                                                           ccGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTACTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg     <  1:1411270/150‑1 (MQ=255)
                                                                                                                                             ggggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg   <  1:1613061/150‑1 (MQ=255)
                                                                                                                                              gggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg  <  2:2074121/150‑1 (MQ=255)
                                                                                                                                              gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg  <  2:2074116/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
ACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGG  >  NC_000913/164431‑164720

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