Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	75,079	0	C	T	46.7%	54.0 / 49.8	60	L134L (CTG→CTA)	thiB	thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (4/28);  new base T (28/0);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 3.55e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTG  >  NC_000913/74931‑75229                                                                                                                                                      |                                                                                                                                                         gCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAg                                                                                                                                                        <  1:1113990/151‑1 (MQ=255)      gCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGttt                                                                                                                                                     >  2:573909/1‑149 (MQ=255)          gtgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTtttttttt                                                                                                                                               >  2:331881/1‑151 (MQ=255)            ggggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTg                                                                                                                                              <  2:671255/150‑1 (MQ=255)             gggCGTCATCGCCATAGCCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCATGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTc                                                                                                                                            <  2:1220590/151‑1 (MQ=255)              ggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTATTTGTtttt                                                                                                                                               >  2:841075/1‑147 (MQ=255)                cGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTtttttttttt                                                                                                                                             >  1:245032/1‑145 (MQ=255)                    aTCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTAGTttttttttttttta                                                                                                                                         >  1:869197/1‑142 (MQ=255)                    aTCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAAc                                                                                                                                     <  1:67074/150‑1 (MQ=255)                      cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCGGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGaa                                                                                                                                  <  1:1106416/151‑1 (MQ=255)                             gACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTTTTTTTTTTTTTTAT‑AAAAGAAgggg                                                                                                                              >  2:850445/1‑146 (MQ=255)                                tttttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTGTTCTTGTCTT‑AAAAGAAGGCGGAg                                                                                                                           >  1:699927/1‑148 (MQ=255)                                  tttGCACCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGTTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGc                                                                                                                       <  2:560240/150‑1 (MQ=255)                                      cATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTTAT‑AAAAGAAAGGGAAGGAggc                                                                                                                      >  1:241516/1‑145 (MQ=255)                                      cATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTGTCAT‑AAACaaaaggggaggaagca                                                                                                                     >  1:140278/1‑135 (MQ=255)                                         ccATAGCAACAGACCCAGCCCCGGTTTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCa                                                                                                              <  2:867874/151‑1 (MQ=255)                                                aaCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTCT‑AAAACAAAGGGAAGGAGCCA‑TAATCaaa                                                                                                            >  1:728639/1‑146 (MQ=255)                                                aaCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTCTTGTCCT‑AAACGAAGGGGAAGGAGGCACTATTCaaaag                                                                                                          >  1:839925/1‑147 (MQ=255)                                                 aCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGGTTTTTTTTTTTAT‑AAAAGAAGGCGGAGGAGACA‑TAATCaaa                                                                                                            >  1:527434/1‑145 (MQ=255)                                                   aGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCAT‑AAACGAAGGGGAAGAAGCCC‑TAATTAAAAGGTACg                                                                                                     >  1:553558/1‑150 (MQ=255)                                                      cccAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTCGTTTTTTTTTTTCAA‑AAAAGAAAGGGAAGTAGCCA‑Taaataaaacgttcgcaa                                                                                                  >  1:1270301/1‑135 (MQ=255)                                                       ccAGCCCCGGTGTACTGTTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGaa                                                                                                   <  2:230072/148‑1 (MQ=255)                                                        cAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTTCT‑AAAACAAAGCGAAGTAGCCA‑TAATCAAACCGGAAGaaa                                                                                                  >  2:879664/1‑148 (MQ=255)                                                           ccctGGTGTACGGGTGCGCGGATCCTGAAAAACCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTc                                                                                             <  2:630623/150‑1 (MQ=255)                                                            cccGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCAT‑TAAAGAAAGCGGAGTAACCC‑caaacaaaacggaagaaagggt                                                                                              >  1:146140/1‑134 (MQ=255)                                                                 tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTGTTTTTTTCAT‑AAAAGAAGGGGGAGAAACCA‑AAATCAAAAAGGACGAAAGGGTGattatt                                                                                       >  1:1234387/1‑150 (MQ=255)                                                                  gtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTc                                                                                      <  2:66489/150‑1 (MQ=255)                                                                    aCTGGTGCGGGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCCAACGGTACGAAAGTGTCATTATTCCa                                                                                    <  1:287625/150‑1 (MQ=255)                                                                      tGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTTAT‑AAAAGAAAGCGAAGTAACCA‑AAAACAAACCGGACGAAAgggtgcatatttttgcc                                                                                 >  1:805952/1‑138 (MQ=255)                                                                       ggTGCGGGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAg                                                                                   <  2:812054/148‑1 (MQ=255)                                                                       ggTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAgc                                                                                  <  2:484667/149‑1 (MQ=255)                                                                       ggTGCGCGGATCCTAAAAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAgc                                                                                  <  2:1010136/149‑1 (MQ=255)                                                                           cgcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGccgccg                                                                             <  1:123032/150‑1 (MQ=255)                                                                            gcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTTTTTTTTCTTGTTCT‑AAAATAAAGCGAAGTGGGCA‑TAATTAAACAGGAAGAAAGTGTCATTATTTCAgcaccc                                                                              >  2:762752/1‑144 (MQ=255)                                                                              ggATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGa                                                                          <  1:251679/150‑1 (MQ=255)                                                                                   cTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTTT‑AAAAAAAGGCGGAGGAACCA‑TAATTAAAAAGGACGAAAAGGGCATTATTTCCGCCCCCCGGGACggt                                                                     >  1:404308/1‑148 (MQ=255)                                                                                    tGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCAT‑AAAAAAAAGCGAAGTAGGCA‑TAAACAAAAAGGACGAAAGGGTCTTTATTTCAGCCGCCggggaaggt                                                                     >  2:717761/1‑143 (MQ=255)                                                                                    tGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTGTTTTTTTCATAAAACAAAGGGGAGGAGGCCA‑‑AAACAAAAGGGACGGAAAAGGCATTATTTCAGCCGCCggg                                                                           >  1:40357/1‑143 (MQ=255)                                                                                     gATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTa                                                                    <  1:276156/149‑1 (MQ=255)                                                                                      aTAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAAc                                                                  <  1:380651/150‑1 (MQ=255)                                                                                      aTAAACCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAAc                                                                  <  2:910950/150‑1 (MQ=255)                                                                                        aaaTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACg                                                                 <  2:147610/149‑1 (MQ=255)                                                                                         aaTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTTTTAT‑AAAAGAAGGCGAAGGAGGCA‑TAATTAAAAGGGACGAAAGAGGCATTATTTCAGCCGCCGGGGACGTTAAAGGc                                                               >  2:140330/1‑150 (MQ=255)                                                                                              cccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATcc                                                           <  2:1238459/149‑1 (MQ=255)                                                                                               ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTCTTTTCAT‑AAAAGAAGGCGAAGGAGCCA‑TAATCAAACGGTAAGAAAGAGTCATTTTTTCAGCCCCCGGGGACGTTAACGGCATcccc                                                         >  1:572666/1‑148 (MQ=255)                                                                                                      ttttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTTTCAT‑AAAAGAAGGGGAAGGAGCCA‑TAATCAAAAAGGACGGAAGTGTCATTATTTCCGCCCCCCGGGACGTTAAAGGCATCCCCCGcccc                                                   >  2:932900/1‑147 (MQ=255)                                                                                                             cGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCtttt                                          <  1:498600/151‑1 (MQ=255)                                                                                                             cGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCAT‑AAAAGAAGGGGAAGAAGCCC‑TAATTAAAAGGTACGAAAATGTCATTATTTCAGCCGCCGGGGAAGATTACGGCATCCGCCGCCACACCCCcttt                                          >  1:1118868/1‑151 (MQ=255)                                                                                                                 ctcAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTTTCAT‑AAAAGAAAGGGGAGTAGCCA‑TAATCAAACGGGACGAAAGTGTCATTATTTCAGCCGCCGGGGACGGTAACGGCATCCGCTGCCACACCCCctttt                                         >  1:851164/1‑147 (MQ=255)                                                                                                                     aCCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTGTCCT‑AAAACAAAGCGAGGAAGCCC‑TAATCAAAACGTACGAAAATGTCATTATTTCAGCCGCCCGGAAAGTTAACGGCATCCGCTGCCCCACCCCCTTTTGCaaa                                    >  1:263678/1‑149 (MQ=255)                                                                                                                       cAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGt                                <  1:140332/151‑1 (MQ=255)                                                                                                                               tCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGttt                         <  2:947259/150‑1 (MQ=255)                                                                                                                                        gtgGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcgg                <  1:922360/150‑1 (MQ=255)                                                                                                                                         tgGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGgcggcg              >  2:283775/1‑151 (MQ=255)                                                                                                                                          ggC‑GGTTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGgcggcg              <  2:1010835/150‑1 (MQ=255)                                                                                                                                          ggC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGt             <  2:699918/151‑1 (MQ=255)                                                                                                                                            c‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt           >  1:1046592/1‑151 (MQ=255)                                                                                                                                                gTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAc        >  2:625234/1‑151 (MQ=255)                                                                                                                                                   tttCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAgg     >  1:677727/1‑151 (MQ=255)                                                                                                                                                   tttCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAg      <  2:630702/150‑1 (MQ=255)                                                                                                                                                      cAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGgttg  <  1:690096/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         GCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCA‑TAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTG  >  NC_000913/74931‑75229

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