Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG43.9% 97.2 / 10.7 57N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (8/24);  new base G (25/0);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.55e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677889‑2678184
                                                                                                                                                       |                                                                                                                                                 
accaGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGt                                                                                                                                                   <  1:7509/151‑1 (MQ=255)
   aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGCAGA‑GGgttt                                                                                                                                                >  1:516961/1‑151 (MQ=255)
       gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTcc                                                                                                                                              <  2:291674/149‑1 (MQ=255)
         ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGGAGA‑GGTTTTTCGcc                                                                                                                                           >  1:76012/1‑150 (MQ=255)
             gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACt                                                                                                                                        <  1:706466/149‑1 (MQ=255)
              cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGGAGA‑GGGTTTCCGCCACTGt                                                                                                                                      >  2:325364/1‑150 (MQ=255)
              cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCCGA‑GGGTTTTCGcc                                                                                                                                           >  1:490393/1‑145 (MQ=255)
               gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCCAGGCAGA‑GGGTTTTCGcccc                                                                                                                                         >  1:734073/1‑144 (MQ=255)
               gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTc                                                                                                                                     <  1:784457/150‑1 (MQ=255)
               gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTc                                                                                                                                     <  1:50276/150‑1 (MQ=255)
               gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCa                                                                                                                                    <  1:292762/151‑1 (MQ=255)
                    gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTCGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCGCGGGAGA‑GGGTTTTCGCCCCTGTCAGccc                                                                                                                                >  1:18903/1‑150 (MQ=255)
                     cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACCGTCAGccc                                                                                                                                >  2:7729/1‑149 (MQ=255)
                       tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCCCTTTCAGCCCCTcc                                                                                                                            >  1:317172/1‑151 (MQ=255)
                             tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCTTccacc                                                                                                                         >  2:229861/1‑148 (MQ=255)
                                 aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCCGACtt                                                                                                                   >  2:275822/1‑150 (MQ=255)
                                  cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTaa                                                                                                                 <  2:202906/151‑1 (MQ=255)
                                   aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTaa                                                                                                                 >  2:458469/1‑150 (MQ=255)
                                       tgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTc                                                                                                              <  1:532389/149‑1 (MQ=255)
                                       agttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCCTCCACCAGACTTAACTcc                                                                                                             >  2:445665/2‑150 (MQ=255)
                                           cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATa                                                                                                        <  1:111850/151‑1 (MQ=255)
                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgcg                                                                                                 <  2:574943/151‑1 (MQ=255)
                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCCTCCACCAGACTTAAATCCTGATAAACcgc                                                                                                  >  2:508744/1‑150 (MQ=255)
                                                       cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTCCCCCCCCGTCAGCCCCTCCCCCCCACCTCACTCCTGATAAAACGCGCCCAc                                                                                            >  1:271128/1‑151 (MQ=255)
                                                           gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCCCCGTCAGCCCTTCCCCCAGACTTAACTCCTGATAAAACCCCCGCAccccc                                                                                        >  1:303535/1‑151 (MQ=255)
                                                              cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCCCCACCGTCAGCCCCTCCACCAGACTTAACTCCTGATAAAACGCGCGCcccccc                                                                                        >  1:218234/1‑148 (MQ=255)
                                                               aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCCCCGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCCCGCACCCCCAGtt                                                                                    >  2:271003/1‑151 (MQ=255)
                                                                gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTc                                                                                   <  2:98586/151‑1 (MQ=255)
                                                                gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTc                                                                                   <  1:845439/151‑1 (MQ=255)
                                                                gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTc                                                                                   <  1:645959/151‑1 (MQ=255)
                                                                 cgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCa                                                                                  <  2:865424/151‑1 (MQ=255)
                                                                    cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCCGACTTAACTCCTGATAAAACGCGCGCACCCCCAGATCAg                                                                                 >  1:11077/1‑149 (MQ=255)
                                                                          gCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCTTCCACCCGACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCAc                                                                           >  2:54564/1‑149 (MQ=255)
                                                                               gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCCCCAGACTTAACTCCTGATAAACCCCGCGCACCCCCAGTTCAGCGGCACGGCGAg                                                                     >  2:150158/1‑150 (MQ=255)
                                                                                ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCg                                                                       >  2:194332/1‑147 (MQ=255)
                                                                                  gggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                                  <  2:629075/150‑1 (MQ=255)
                                                                                     cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAAACGCGCGCCCCCCCCGTTCCGCGGCACGGCGAGGCAGc                                                                >  2:876318/1‑149 (MQ=255)
                                                                                       tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCCCCCCCCGTTCAGCGGCACGGCGAGTCcgcg                                                               >  1:209211/1‑148 (MQ=255)
                                                                                             gcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGt                                                      <  1:124824/151‑1 (MQ=255)
                                                                                               cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGCACGGCGGGTCAGCGGAGCTTCGt                                                      >  1:839031/1‑149 (MQ=255)
                                                                                                aTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGGAGAGGGTTTTCCCCCCCGTTCACCCCTT‑CACCCCACTTAACTCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGCACGGCGGGTCAGCGCAGCTTCGTCAc                                                   >  1:411350/1‑151 (MQ=255)
                                                                                                   tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCCCCAGACCTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCCCGGCGGGTCCGCGGAGCTTCCTCAcc                                                  >  1:334122/1‑149 (MQ=255)
                                                                                                    gCAGGTAATCAATCATGCCGTTGCGGCGGGGGCAACGACCCAGGCAGA‑GGGTTTTCGCCACTGTCAGACCCTCCACCAGAATTAACTCCTGATAAAACCCGCGCACCCCCCGTTCTGCGGCACAGGGGGTCAGCCTAGCCTCGTTACCtt                                                >  1:589086/1‑150 (MQ=255)
                                                                                                          aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTc                                          <  1:760280/150‑1 (MQ=255)
                                                                                                              aatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtc                                     <  1:700490/151‑1 (MQ=255)
                                                                                                                      cGTTGCGTCGGGGGCCCCTTCCGAGGGAGA‑GGGTTTCCCCCCCGTTTAGCCCTTCCACCAAAAGTTAATCCTGATAAAACCCCCGCACCCCCAGTTCAGCGGCACCGCGGGGCAGCGTAGCTTCGTCACCTTCCAGGcgcgggccgc                                 >  1:237889/1‑141 (MQ=255)
                                                                                                                           cggcggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                           >  2:337146/1‑148 (MQ=255)
                                                                                                                            ggcggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                           >  1:99783/1‑147 (MQ=255)
                                                                                                                               ggggCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCTACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                    >  2:197242/1‑151 (MQ=255)
                                                                                                                               ggggCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                    >  1:137980/1‑151 (MQ=255)
                                                                                                                                    cACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc               <  1:458427/151‑1 (MQ=255)
                                                                                                                                        gACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc           <  1:288578/151‑1 (MQ=255)
                                                                                                                                          ccGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg         >  1:704106/1‑151 (MQ=255)
                                                                                                                                           cGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg        <  2:481129/151‑1 (MQ=255)
                                                                                                                                           cGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg        >  1:745150/1‑151 (MQ=255)
                                                                                                                                            gAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCggg       <  1:275794/151‑1 (MQ=255)
                                                                                                                                             aGGCTGA‑GTTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCTCCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCCCCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCggg       <  1:462844/150‑1 (MQ=255)
                                                                                                                                                 aga‑gGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc  <  2:18903/151‑1 (MQ=255)
                                                                                                                                                 aga‑gGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc  <  1:7729/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
ACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677889‑2678184

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