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.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                       |                                                                                                                                                     
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                       >  1:669761/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGCCAAATAATGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                       <  2:718291/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAATGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  2:521843/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  2:955095/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  2:2051748/150‑1 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  1:426749/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                             >  2:161922/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                             >  1:590371/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                             >  2:161887/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                             >  1:113416/1‑150 (MQ=255)
           aaaaaGAGATGAAAGTGAGTTATCCCAAACATGACTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGTGCATAAGCGCGCAc                                                                                                                                            >  1:1410752/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCg                                                                                                                                           >  1:310426/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  1:1462527/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  1:2074116/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  1:2074121/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  1:310586/1‑150 (MQ=255)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgggcccacaccc                                                                                                                              >  1:588338/1‑150 (MQ=255)
                              ttATCCCAACTATGTCTGGCCGGACCACCCGGAAAATACTCCCCCGCGGGGACGGACCAAAATTTATTCTTATGTTTGTTTTTCTCCCTCTCCCTTGGGGTGTGGGCGGGGTTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:47374/150‑1 (MQ=255)
                               tATCCCAAACTTGTGGGGCCGCACGACCCGTAAAAAACTGCCCCGACGCCACGAACAAAAATTTATTCGTATGGTGGTTTTTCCCCCCCCCCGGGGGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                                        <  2:252678/150‑1 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCAAAaccacgc                                                                                                                    >  2:405775/1‑146 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCCCCCCAAcccacac                                                                                                                    >  2:1613061/1‑146 (MQ=37)
                                    aaaaCTTGTGGGGCCGTACGACCCGGAAAATACTGCACCGCCGGTACGGCCTAAATTTTTTTCTTATGTTTGTTTTTCTCCCTCTCCGTGTGGGTGGGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  1:1567787/149‑1 (MQ=255)
                                                  cgcaccccccGGAAAAAACTGCCCCGCCGCGCCGGACTAAAAAGTATTCGTACGGTTGGTTTCCACCCCCCCCCTGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                    <  2:1480563/143‑1 (MQ=255)
                                                  ccgcacccccGGCAAAAACCGCCCCCCACGCGACGCACAAAAATTTTTGCTTATGTTTGGTTTCCCCCCCCCCCCGGGGGGTGGGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                     <  1:628169/143‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCGTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                    <  2:1176798/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                    <  2:2142941/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                    <  2:2143259/150‑1 (MQ=255)
                                                     acgccCCGGAAAACACTCCACCGCCGCGACGGACTAAAAAGTATTCGTACGGTTGTTTTACCCCCTCCCCCTGGGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                 <  2:166772/150‑1 (MQ=255)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAAACCACACCCACCAAGGGGGGGGGGccc                                                                                                >  2:376195/1‑150 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCAACCGCACCCGAACCCTCAACCGCCTGGGTGGGGGGACCGATCg                                                                                           >  2:1557834/1‑150 (MQ=25)
                                                                         aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                              <  2:1304957/150‑1 (MQ=255)
                                                                          ccgccgCCAGGGACTAAAATGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:2105829/150‑1 (MQ=255)
                                                                            gacgcgccgGACTAAAAATTATTCGTATGTGTGTTTTTCACCCACTCCCTTGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                           <  2:585944/150‑1 (MQ=255)
                                                                                 gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAACCACAACCCCCAGGGGGAGGGGGCCGCGAGAGCACCAAATTTTAGTGatat                                                                      >  2:269281/1‑150 (MQ=18)
                                                                                           aatgTAATCGAACGGTGGGTTTTCTCCCTCACCGTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                            <  2:1086386/147‑1 (MQ=255)
                                                                                           aaaGTATCCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                            <  2:982381/150‑1 (MQ=255)
                                                                                                atgcGTATGGTTGTTTTTCCCCCTCCCCCTGTGGTTGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                       <  1:1934653/147‑1 (MQ=255)
                                                                                                   cgaaggtttggttTCCCCCCCCCCCGTGGGGTTGGGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  1:1020039/143‑1 (MQ=255)
                                                                                                    gTAACGGTGGTTTTACTCCCACACCCTGGGGGTGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  1:1020049/150‑1 (MQ=255)
                                                                                                      aaGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  2:2069167/150‑1 (MQ=255)
                                                                                                                ttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGTCCCACACCCAAAACCACTCCCGCAACGGGGCGGGGGCCCATCGTGCACAAAATTTTTGGGTTATATTATTATTTGTTGTTACAGAagaaaagaaa                                       >  2:1150242/1‑150 (MQ=12)
                                                                                                                   ctccctcgccctGTGGTTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                     <  2:355570/150‑1 (MQ=255)
                                                                                                                      cctccccgtGTGGGTGTGGGGGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                  <  1:460475/150‑1 (MQ=255)
                                                                                                                           ccctGGGGGAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                             <  2:1366110/150‑1 (MQ=255)
                                                                                                                              tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCCCCCAAAACCACTCCCGCAAGGGGGGGGGGTCCGATAGAGCACCAATTTTTTGAGTTTTATTATTATGTGTTTTTAAAGAATAACAGAAAAAAGGGGtttttg                          >  2:2062183/1‑148 (MQ=255)
                                                                                                                                  ggcgtgggCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      <  1:1996943/145‑1 (MQ=255)
                                                                                                                                    tgagGGCCGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                    <  2:727117/149‑1 (MQ=255)
                                                                                                                                      tgggCGGGGTTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                  <  2:806249/149‑1 (MQ=255)
                                                                                                                                      tgggCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                  <  2:805941/149‑1 (MQ=255)
                                                                                                                                      tgggCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                  <  2:806183/149‑1 (MQ=255)
                                                                                                                                         gcgggggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  1:1263199/147‑1 (MQ=255)
                                                                                                                                          gcGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg              <  1:1411305/149‑1 (MQ=255)
                                                                                                                                          ccGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTACTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg              <  1:1411270/150‑1 (MQ=255)
                                                                                                                                            ggggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:1613061/150‑1 (MQ=255)
                                                                                                                                             gggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg           <  2:2074121/150‑1 (MQ=255)
                                                                                                                                             gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg           <  2:2074116/150‑1 (MQ=255)
                                                                                                                                                    gccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  1:376195/148‑1 (MQ=255)
                                                                                                                                                      cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  2:826748/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  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: