Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC30.4% 45.6 / 23.7 56intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base C (2/15);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 3.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                       |                                                                                                                                                     
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:955095/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  2:2051748/150‑1 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  1:426749/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                             >  2:161922/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                             >  1:590371/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:1462527/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: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)
                                                                                                                                                    gccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  1:376195/148‑1 (MQ=255)
                                                                                                                                                      cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  2:826748/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCC‑TCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑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: ‑