Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,4640AT44.1% 81.1 / 29.9 59R9R (CGT→CGAcitEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (6/27);  new base T (26/0);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 2.32e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGT  >  NC_000913/650318‑650614
                                                                                                                                                  |                                                                                                                                                       
aTGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACgtt                                                                                                                                                     <  2:946531/151‑1 (MQ=255)
aTGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACgtt                                                                                                                                                     <  2:315668/151‑1 (MQ=255)
     gcgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGc                                                                                                                                                <  1:554810/151‑1 (MQ=255)
     gcgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCgtttttt                                                                                                                                                 >  2:109623/1‑149 (MQ=255)
     gcgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCgttgtt                                                                                                                                                  >  1:602703/1‑149 (MQ=255)
                 tttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGAGCGAGGTTTTTGTTTGTGCGGCGAAGCgg                                                                                                                                      >  1:421224/1‑149 (MQ=255)
                  ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGGTTTTCGTTGTTTTAGCGCAGcag                                                                                                                                      >  1:686012/1‑146 (MQ=255)
                  ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTTTTGCAGCGAAACgg                                                                                                                                      >  1:910405/1‑148 (MQ=255)
                  ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTTTTTCAGCGAAACGGaa                                                                                                                                    >  1:972756/1‑150 (MQ=255)
                  ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCCCGTTTTTCGGTGTTGGACCGAAGCGGaa                                                                                                                                    >  1:110477/1‑150 (MQ=255)
                    ttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTTGTTGTTGCAGCGAAACGGaaa                                                                                                                                   >  1:788639/1‑149 (MQ=255)
                    ttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAAt                                                                                                                                  <  1:156997/150‑1 (MQ=255)
                      cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCa                                                                                                                                <  2:857530/150‑1 (MQ=255)
                      cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGCGGGTTTATGGTTTTTCAGCGAAGCGGGAAACa                                                                                                                                >  2:692046/1‑150 (MQ=255)
                           aaTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTTCGTTGTTGCAGCGAAGCGGAAATCAtt                                                                                                                              >  1:270628/1‑147 (MQ=255)
                           aaTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGAAGAGGGAAACa                                                                                                                                >  1:560850/1‑145 (MQ=255)
                              gCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAtc                                                                                                                       <  1:843016/151‑1 (MQ=255)
                                tACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGGCGCGAAGCGGGAAACATTGGGAAtcc                                                                                                                      >  1:441993/1‑150 (MQ=255)
                                 aCGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTTTGTTTGTGGAACGAAACGGAAATCATTTGCaa                                                                                                                         >  1:761444/1‑146 (MQ=255)
                                 aCGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAtcctc                                                                                                                    <  2:602724/151‑1 (MQ=255)
                                     aGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTGTTGCAGCGAAGCGGAAAACATTGGCAATCCTccc                                                                                                                  >  1:530677/1‑149 (MQ=255)
                                         ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGGGGCGGAGCGGGAAACAATGGGCAACCTccc                                                                                                                  >  2:166612/1‑145 (MQ=255)
                                          tCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGCGGTTTTCGTTGTTGCAGCGAAGCGGAAATCATTTGCAATCCCCCCCTgg                                                                                                              >  2:393797/1‑148 (MQ=255)
                                             aGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTTGTTTTTGCAGCGAAACGGAAAACAATGGCAATCCTCCCATGGGag                                                                                                           >  1:920608/1‑148 (MQ=255)
                                               gTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTACGTTGTTTGGGAGAAACGGGCAACATTTGCAATACTCCCCTGGCCGAGcc                                                                                                      >  2:226788/1‑150 (MQ=255)
                                                 cAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCgg                                                                                                    <  1:1052723/150‑1 (MQ=255)
                                                  aaaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTTGTTGTTGCAGCGAAGCGGAAAACATTGGCAATCCTCCCCTGGCAGAGCCggg                                                                                                   >  1:260425/1‑150 (MQ=255)
                                                   aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTTGTTGTTGCAGCGAAAAGGGAAACATTTGCAATCCTCCCCTGGCCGAGCCggg                                                                                                   >  1:461493/1‑149 (MQ=255)
                                                       tCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTTCGGCGAAGAGGGAAAAATTGGGAAACCCCCCCTGGGAGAGGCCgggga                                                                                                 >  2:1113631/1‑146 (MQ=255)
                                                        cAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCg                                                                                             <  2:972792/150‑1 (MQ=255)
                                                                cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGCGTTTTTTGGTGTTTGCGCGGAGAGGAAAAAATTTGGAATCCACCCATGGCAGAACCGGGGTAcccccgg                                                                                         >  1:848206/1‑141 (MQ=255)
                                                                      ggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                              <  2:990778/151‑1 (MQ=255)
                                                                       gTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                                                                             <  2:867015/151‑1 (MQ=255)
                                                                                   gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                                 <  2:1070986/151‑1 (MQ=255)
                                                                                   gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                                 <  2:1085887/151‑1 (MQ=255)
                                                                                   gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                                 <  2:788663/151‑1 (MQ=255)
                                                                                      gCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGt                                                               <  2:173967/150‑1 (MQ=255)
                                                                                       cTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTGTCGCAGCGAAGCGGAAAACATTGGGAATCCTCCCCTGGCAGAGCCGGGGTACCGCCGGCGCGGGCCAGCAGGGCTTCCCGTc                                                                 >  1:1030538/1‑147 (MQ=255)
                                                                                       cTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGAAGCGGAAAAAATTGGCAATCCCCCCCTGGCCGAGCCGGGAGACCGCCGGCGCGGGCCCGCAGGGCTTcc                                                                     >  2:50296/1‑143 (MQ=255)
                                                                                        tGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGAGGTTTTCGTTTTTTCAGCGGAGCGGAAAACAATTGCAATCCCCCCCTTGCAGAGGCGGGGAACCGCCGGGGCGGGCCCGCAGGGGTTcc                                                                     >  2:481958/1‑142 (MQ=255)
                                                                                         gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc                                                           <  2:842926/151‑1 (MQ=255)
                                                                                          aCCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACg                                                          <  1:244677/151‑1 (MQ=255)
                                                                                               cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTaaaa                                                     <  1:734053/151‑1 (MQ=255)
                                                                                               cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTaaa                                                      <  1:810602/150‑1 (MQ=255)
                                                                                                gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAAt                                                    <  2:202560/151‑1 (MQ=255)
                                                                                                             cAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGGTTTTTCAGCGAAACGGGAATCATTGGCAAACCCCCCCTGGCAGAGCCGGGGTACCCCTGGCGCGGGGCAGCAGGGCCTCCCGTCGTGCACGGAAAAATCAGTCCAGTg                                         >  2:139279/1‑149 (MQ=255)
                                                                                                                gCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCt                                    <  2:530487/151‑1 (MQ=255)
                                                                                                                  acaaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                                  <  1:783356/151‑1 (MQ=255)
                                                                                                                  acaaacaaCATGCTGCGGCGGGGGGGGGGTTTTCGTTGTTTCAGCGAAACGGGAAACAATTGCAAACCTCCCCTGGCCGAGCCGGGGGACCCCCGGGGCGTGGCCGCCGGGCTTCCAGGCGGGCCCGGAAAAAGCAGTcccgggggccct                                   >  2:18558/1‑140 (MQ=255)
                                                                                                                    aaacaaCATGCTGCGGCGGGGGCGGGGTTTTTTTTGTTTCCGCGAAAAGGAAAACAATGGGAAACCCCCCCTGGCAGAGCCGGGGTACCCCCGGCGCGTGCCCGCAGGGGTTCCAGTCGGGCACGTAAAAAACAGGCCAGAGGGCCctt                                  >  2:478748/1‑149 (MQ=255)
                                                                                                                                gcggcgGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACAtt                      <  2:734528/149‑1 (MQ=255)
                                                                                                                                   gcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCt                 <  2:685336/151‑1 (MQ=255)
                                                                                                                                   gcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCt                 >  1:183517/1‑151 (MQ=255)
                                                                                                                                     gggTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGt               <  2:761464/151‑1 (MQ=255)
                                                                                                                                       gTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc             <  1:299804/151‑1 (MQ=255)
                                                                                                                                           gAGTTTTACGTTGTCGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgccg         <  2:1030575/151‑1 (MQ=255)
                                                                                                                                                ttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCAc    >  1:824935/1‑151 (MQ=255)
                                                                                                                                                ttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCAc    >  1:1092303/1‑151 (MQ=255)
                                                                                                                                                 tACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACg   >  1:24143/1‑151 (MQ=255)
                                                                                                                                                  aCGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt  >  2:585674/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
ATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGT  >  NC_000913/650318‑650614

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: