Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG44.7% 125.8 / 15.9 76N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (13/29);  new base G (34/0);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 3.66e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGC  >  NC_000913/2677889‑2678189
                                                                                                                                                        |                                                                                                                                                      
accgtctgctTTTGCTCCCGGCTTAGTGTTGAAACAAGTTGTTCAGGACTAACGACAACGCCCAGCGCCTGTAAGCAGCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCTGCGGGG‑CCACTGACCGAGGCAGAGG‑t                                                                                                                                                        <  2:762570/146‑1 (MQ=255)
  caGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑ttt                                                                                                                                                      <  1:1285449/151‑1 (MQ=255)
       gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGc                                                                                                                                                 <  1:300260/151‑1 (MQ=255)
         ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTTACCAAGGGAGAAG‑GTTtcctcc                                                                                                                                                >  1:191673/1‑150 (MQ=255)
           ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGGACGAGGCAGAGG‑GTTTTCGcc                                                                                                                                                >  1:1065847/1‑148 (MQ=255)
           ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCAc                                                                                                                                              <  1:132053/150‑1 (MQ=255)
              cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CGCCCGAGCCAGGGAGAGG‑GTTTCCGCCAcct                                                                                                                                            >  2:11812/1‑147 (MQ=255)
                ccGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAg                                                                                                                                        <  1:328684/151‑1 (MQ=255)
                   ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGGAGAGG‑GTTTTCGCCCCTGTCAGccc                                                                                                                                     >  1:746957/1‑151 (MQ=255)
                   ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGccc                                                                                                                                     <  2:920440/151‑1 (MQ=255)
                    gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑GCAATGACCGCGGCAGAGG‑GTTTTCGCCACAGTCAGccc                                                                                                                                     >  1:1025507/1‑150 (MQ=255)
                        aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCATCGCCTGTAAGCAGCGTTGGGCATCCTGCGCCATTTGCAGTTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTcc                                                                                                                                 <  2:1089815/150‑1 (MQ=255)
                          cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTGACCCAGGCAGAGG‑GTTTTCGCCACTGCCAGCCCCTtcac                                                                                                                               >  2:732051/1‑150 (MQ=255)
                           gTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTccacca                                                                                                                             <  1:34256/151‑1 (MQ=255)
                             tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTGACCGAGGGGGAGG‑GTTTTCGCCACTCTGAGCCCCTccacc                                                                                                                              >  2:863623/1‑148 (MQ=255)
                               aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGTTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACt                                                                                                                         <  2:624764/151‑1 (MQ=255)
                                 aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTa                                                                                                                       <  1:969950/151‑1 (MQ=255)
                                  cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTGACCCCGGCCGAGG‑GTTTTCGCCACCGTCAGCCCCTCCACCCGAcct                                                                                                                        >  2:1316161/1‑147 (MQ=255)
                                   aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTGACCGAGGCCGAGG‑GCTTCCTCCCCCCTCAGCCCCTCCCCCCGACTTaa                                                                                                                      >  2:696432/1‑150 (MQ=255)
                                    aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACt                                                                                                                    <  2:1034766/151‑1 (MQ=255)
                                     gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACt                                                                                                                    <  2:105341/150‑1 (MQ=255)
                                          tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTGACCCAGGCAGAGG‑GTTTCCCCCCCTGTCAGCCCCTCCCCCAAACTTAAATCCTGAt                                                                                                              >  1:410181/1‑151 (MQ=255)
                                           cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGTG‑CCACTCGCCGGGGCGGAGG‑GTTTCCGCCACCGTCCGGCCTTCCCCCCGACTTTACTcc                                                                                                                  >  2:960326/1‑146 (MQ=255)
                                            aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACCGCCCGAGGGCGAGG‑GTTTCCGCCCCCGTCCGCCCTCCCCCCAgcattacaccc                                                                                                                  >  2:704058/1‑135 (MQ=255)
                                             gggcTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCTTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCTGCCCTTCCACCAGACTTAACTCCTGa                                                                                                               <  2:1322728/144‑1 (MQ=255)
                                             ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                           <  1:43915/151‑1 (MQ=255)
                                               aCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑GCCCTGGCCGAGGCAGAGG‑GTTTTTTACCCCGTTAGCCCCTCCCCCCGACTTAAATCCTCCTaaa                                                                                                           >  2:290345/1‑149 (MQ=255)
                                               aCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCCGACCGAGGCGGCGG‑GTTTTCGCCCCCGTCAGCCCCTCCACCCAACTTTACCCCTGATaaa                                                                                                           >  2:900768/1‑149 (MQ=255)
                                                cTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAcc                                                                                                         <  2:1339164/150‑1 (MQ=255)
                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTTAACCAGGCAGAGG‑TTTTCCGCCCCTGTCCGCCCCTCCACCCCACTTAAATACTGATaa                                                                                                            >  2:1333017/1‑145 (MQ=255)
                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTGACCGAGGCAGAGG‑GTTTTCGCCCCTGTCAGCCCCTCCCCCCGACTTAACTCCTGATaaaac                                                                                                         >  1:40050/1‑146 (MQ=255)
                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgc                                                                                                       <  1:384048/150‑1 (MQ=255)
                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgc                                                                                                       <  1:1339276/150‑1 (MQ=255)
                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACCCACCCAGGCAGAGG‑GTTTTCGCCCCTGTCAGCCCCTCCCCCCGACCTTACTCCTGATaaa                                                                                                           >  1:372266/1‑146 (MQ=255)
                                                       cAACGCCCAGCGCCTGTAAGCAACGTTGGGCTTCCGGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCgcgcgc                                                                                                   <  2:182312/149‑1 (MQ=255)
                                                        aaCGCCCAGCGCCTGTAAGCAGCGTTGGGAATCCTGCGCCGTTTGCAGGAAATCAATCGTGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAc                                                                                                 <  2:726274/150‑1 (MQ=255)
                                                         aCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAccc                                                                                               <  1:406753/151‑1 (MQ=255)
                                                           gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑GGCCGGAGCCGGGCGAAGG‑GTTTTCTCCGCTGTCTACCCTTCCACCCGAATTAAATCCTTATTAACCGCGCGCcccccc                                                                                             >  1:685300/1‑151 (MQ=255)
                                                             ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAccccc                                                                                             <  1:349694/149‑1 (MQ=255)
                                                                gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGT‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGTCTTAACTCCTGATAAACCGCGCGCACCCCCAGTTc                                                                                        <  2:1224631/151‑1 (MQ=255)
                                                                 cgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑GTTTTCGCCCCTGTCAGCCCTTCCACCCGACTTAACTCCTGATAAACCCCGCGCACCCCCCGTTc                                                                                        >  1:212571/1‑150 (MQ=255)
                                                                   ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCa                                                                                       <  2:604254/149‑1 (MQ=255)
                                                                   ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGc                                                                                     <  2:1094391/151‑1 (MQ=255)
                                                                     tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCGCAAACCG‑GTCAGAGG‑GTTTTCGGCCCCGTCAGCCCCTACAACAGCCCTAACTCATGAGGAAACCCCCGCGccccccgctacacgg                                                                                   >  2:1200549/1‑141 (MQ=255)
                                                                        aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑GTTTCCGCCACTGTCAGCCCTTCCACCAAACTTAACTCCTGATAAACCGCGCGCACCCCCCGtt                                                                                         >  1:1142368/1‑142 (MQ=255)
                                                                           cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑GCACTGACCGGGGCCGAGGTTTTTCCCCCCCTCTGA‑CCCTTCCACCCGACCTAAATCCTGATAAAAAGCGCGCAACCCCAGATCATCGGCAc                                                                                >  2:864735/1‑148 (MQ=255)
                                                                            aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑GTTTTCGCCACTGTCAGCCCTTCCACCAGACTTTACTCCTGATAAACCCCGCGCACCCCCCGTTCAGCGGCACg                                                                               >  1:1227610/1‑148 (MQ=255)
                                                                              cGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAg                                                                          <  2:1176195/151‑1 (MQ=255)
                                                                               gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCCCGGCAGAGG‑GTTTTTGCCCCTGTCAGCCCCTCCCCCAGACCTAAATCCTGATAAAACCCGCCCCCCCCCCCGTCAGCGGCACGGCGAg                                                                          >  2:107054/1‑150 (MQ=255)
                                                                                ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCCAGGCAGAGG‑GTTTTCGCCCCTGTCAGCCCTTCCCCCAGACTTTACTCCTGATAAAACGCCCGCACCCCCCGTTCAACg                                                                                    >  1:1197060/1‑139 (MQ=255)
                                                                                 tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACCGAACGCGGCAGAGG‑GTTTTCCCCACTGTCCGCCCTTCCACCCGACTTAACTCCTGCTAAAACGCGCGCACCCCCCCTTCAGCGGCACGGCGAGTc                                                                        >  2:1142473/1‑150 (MQ=255)
                                                                                        ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCCCTGACCGCGGCAGAGG‑GTTTCCGCCCCTGTCAGCCCCTCCACCAAACCTAAATCCCGATAAAACGCCCCCACCCCCCGGTTAGCGGCACGGGGGGGCCGCGTa                                                                  >  1:925258/1‑149 (MQ=255)
                                                                                          tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑GCCCCGACCCCGGCAGAGG‑TTTTCCGCCACCGTCAGCCCCTCCACCAGACTTTAATCCCGATAAAACCCGCCCCCCCCCCGCTCAGCGGCACGGCGGGGCAGCGTAGCtt                                                              >  1:42907/1‑151 (MQ=255)
                                                                                           gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCtt                                                              <  1:752932/150‑1 (MQ=255)
                                                                                            cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑GCACCGACCGAGGCAGAGG‑GTTTTCCCCCCTGTCAGCCCCTCCACCCGACTTAACCCCTGATAAAACCCGCGCCCCCCCCGCTTAGCGGCACGGCGAggaggcgt                                                                   >  2:40746/1‑137 (MQ=255)
                                                                                             gcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGGAGAGG‑GTTTTCGGCACTGTCAGCCCCTCCACCCGACCTAACTCCTGATAAAACCCGCGCCCCCCCCGTTCAGCGGCACGGCggggacgcggagatt                                                              >  2:1257103/1‑134 (MQ=255)
                                                                                             gcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACCGACCGAGGCAGAGG‑GTTTTCCCCACCGTCAGCCCCTCCACCCGACCTAACTCCTCATAAACCCCGCGCACCCCCCCTTCAGCGGCACGGCGAGGCAGCGTAGCtt                                                              >  1:839188/1‑148 (MQ=255)
                                                                                             gcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACCCACCCCGGCAGAGG‑GTTTTCCCCCCCGTCAGCCCCTTCCCCCGACATAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCGGCACGGCGCGTCAGCGGAGCtt                                                              >  2:932167/1‑148 (MQ=255)
                                                                                                       ggTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAg                                                 <  2:73338/151‑1 (MQ=255)
                                                                                                          aatcaatcATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑GTTTTCGCCACTGGCAGCCCTTCCACCCAACTTAAATCCTGATAAAACCCGCGCACCCCCCGCTCAGCGGCACCGCGGGGCAGCGTAGCTTCGTCACCCTCCAGTc                                               >  1:205081/1‑150 (MQ=255)
                                                                                                            tcaatcaTGCCGTTGCGGCGGGG‑CCACCGACCCAGGCAGAGG‑GTTTCCCCCCCTTTCAGCCCTTCCACCAGACTTAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCGGCCCGGCGAGTCAGCGGAGCTTCGTCACCCTCCAGTcgc                                             >  1:71182/1‑150 (MQ=255)
                                                                                                              aatcaTGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtc                                          <  1:8118/151‑1 (MQ=255)
                                                                                                                 caTGCCGTTGCGGCGGGG‑CCACTGACCGCGGCAGAGG‑GTTTTCGCCACCTTCAGCCCCTCCACCCGACCTTACTCCTGATAAACCGCGCGCCCCCCCCGTTCCGCGGCACCGCGAGGCCGCGCAGCTTCTTCCCCTTCCCATCGCGcctc                                        >  2:1357128/1‑150 (MQ=255)
                                                                                                                  aTGCCGTTGCGGCGGGG‑CCACTGACCGAGGCCGAGG‑GTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCCCGCGCACCCCCAGTTCAGCGGCACGGCGAGGCAGCGTAGCTTCGTTACCTTCCCGTCGCGTCCCGc                                      >  1:937933/1‑151 (MQ=255)
                                                                                                                           cggcggGG‑CCACTGACCGAGGCAGAGG‑GTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAAATCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                                >  2:399642/1‑148 (MQ=255)
                                                                                                                            ggcggGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                                >  2:1004020/1‑147 (MQ=255)
                                                                                                                            ggcggGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                                >  2:1222396/1‑147 (MQ=255)
                                                                                                                             gcggGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                            >  1:1239414/1‑150 (MQ=255)
                                                                                                                              cggGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGata                           <  2:325282/150‑1 (MQ=255)
                                                                                                                               gggg‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                          <  2:353439/150‑1 (MQ=255)
                                                                                                                                     cACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc                    >  1:1039563/1‑151 (MQ=255)
                                                                                                                                        tGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTAt                 >  1:1015357/1‑151 (MQ=255)
                                                                                                                                          aCCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATcc               >  2:681040/1‑151 (MQ=255)
                                                                                                                                            cGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg             >  2:1281958/1‑151 (MQ=255)
                                                                                                                                              aGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc           >  2:237681/1‑151 (MQ=255)
                                                                                                                                                  agagG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc       >  1:1189228/1‑151 (MQ=255)
                                                                                                                                                  agagG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc       >  1:811385/1‑151 (MQ=255)
                                                                                                                                                      g‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTg   <  1:549366/151‑1 (MQ=255)
                                                                                                                                                        ttttCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGc  >  2:790632/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
ACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGG‑CCACTGACCGAGGCAGAGG‑TTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGC  >  NC_000913/2677889‑2678189

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