Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	45.8%	77.4 / 36.4	59	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (26/6);  new base A (0/27);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 2.32e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGT  >  NC_000913/436368‑436665                                                                                                                                                        |                                                                                                                                                      ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGc                                                                                                                                                       <  2:736376/149‑1 (MQ=255)   ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCaaa                                                                                                                                                    >  1:38546/1‑150 (MQ=255)     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATg                                                                                                                                                  <  2:360976/150‑1 (MQ=255)     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATg                                                                                                                                                  <  2:389676/150‑1 (MQ=255)     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGa                                                                                                                                                 <  1:993324/151‑1 (MQ=255)        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGAccg                                                                                                                                              >  1:585847/1‑151 (MQ=255)                     aaCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAgg                                                                                                                                 >  2:385822/1‑151 (MQ=255)                        aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGc                                                                                                                               >  1:522194/1‑150 (MQ=255)                          gcgctgcgctGGGGTCTCTCGGGCGGTGAATTTTACGGTTTGTGTTTCCCTGGGCCGAAACTAACCCGTGGCGCGATGGGTGTGGCTCCGGCACACGGGGG‑GGTCCTTTTTCCC‑TGTAGC‑GGAAAAATGACCGCCGATATCGAAGGGCttt                                                                                                                            <  1:931814/151‑1 (MQ=255)                              tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                        >  2:723746/1‑151 (MQ=255)                              tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                        >  2:557610/1‑151 (MQ=255)                               gcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                       >  1:993561/1‑151 (MQ=255)                                         ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                             >  2:1125419/1‑151 (MQ=255)                                             ctgggggaaaaGATTAGGGGTTTTTTTTCGCTGCCCCGGAACTGACCGGGGGCGCCGGGGGTGGGGC‑CCCCGACCGCTGGGCCGTCCTTTTTCCC‑TGTTTG‑GGGAAAATGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                             <  1:704262/138‑1 (MQ=255)                                               gggggTGAAGTTTCCGAGTTGTGTTTCACTGTGCCGGAACGGACCGGTGGCGCGCTGGATGTGGCTCTCGGACCCTGGGG‑CGTCCCTTTTCCC‑TGTTCG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                        <  2:88668/147‑1 (MQ=255)                                               gggggTGAAGATTAGGGGTTGTTTTTCCCTGTGCCGGAACTAACCCGTGGCGCGCGGGAGGTGCCTCTCGGACACTGGGG‑CGTACCTTTTCCC‑TTTACC‑GGGAAAAGGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                        <  1:418993/147‑1 (MQ=255)                                               gggggTAAAGATTGCGTTTTTTTTTGCACTGCGCGGGAACTGACCCGGGGCGCGTTGGGGGGGCCTCCCGGACCCGGGGG‑CTTACCTTTTCCC‑TTTACC‑GGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                          <  1:1055157/145‑1 (MQ=255)                                               agCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                                       >  1:136883/2‑151 (MQ=255)                                                ggggTGAAGTTTACGAGTTGTGTTTCCCTGGGCCGAAACTGACCCGGGGCGCGCTGGATGTGGCTCCCGGACCCCGGGG‑CTTCCCTTTTCCT‑TGTACG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                                                      <  1:777153/149‑1 (MQ=255)                                                ggggTGAAGATTACGGTTTGTGTTTCCCTGTGCCGGAACTGAACCGTGGCGCGGTGGGTGGGGCTCTCGGACACGGGGG‑GGTCCCTTTTCCC‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                                              <  2:151058/141‑1 (MQ=255)                                                  gggaaagttaaCGATTTGTTTTTCACTGCGCGGAACCGGACCCGGGGCGCGGTGGATGTGGCTCCCGAACCCGGGGG‑CGTCCTTTTTCCT‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                    <  2:465208/141‑1 (MQ=255)                                                    tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                  >  2:677171/1‑151 (MQ=255)                                                      aagttgaggggttgttTTTCCCTGGGCCGGACCTGACCGGGGGCGCGTTGGATGTGGCTCCCGGACCCCTGGG‑CGTCCCGTTTCCT‑TGTTTG‑GGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                                 <  2:963728/140‑1 (MQ=255)                                                           tccgtgttttttttCACTGTGCCGAACCTGACCCGTGGCGCGTTGGTTGTGCCCCCCGGACCCGGGGG‑CTTCCTTTTTCCC‑TGTATG‑GGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                            <  1:637546/139‑1 (MQ=255)                                                             cGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                                         >  1:147737/1‑151 (MQ=255)                                                                  tgtgTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                     >  1:72690/1‑150 (MQ=255)                                                                   gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                    >  2:988751/1‑150 (MQ=255)                                                                      tttCCCTGAGCCGGAATGGATCCGCGGCGCGTTGGTTGTGCCTCTGGGACCCTGGGG‑CGTCCTTTTTCCT‑TGTACC‑GGGCAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                 <  1:974211/150‑1 (MQ=255)                                                                      tttCACTGTGCCGGAACGGACCGGGGGCGCGTGGGATGGGCCTCCGGGCCCCGTGGG‑CGTCCCTTTTCCC‑TGTAGC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                    <  1:119682/147‑1 (MQ=255)                                                                          aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                            >  2:172411/1‑151 (MQ=255)                                                                          aCTGGGCCGGAACGGACCCGTGGCGCGCTGGTTGTGGCTCTCGGACACCGGGG‑GGTCCCGTTTCCT‑TTTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                            <  2:55678/151‑1 (MQ=255)                                                                             gggCCGGAACTGACCGGGGGCGCGTTGGATGTGTCTCTCGGCCCCTTGGG‑CGTCCTGTTTCCT‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                          <  2:470269/148‑1 (MQ=255)                                                                                ccGAAACTGACCGGGGGCGCGTTGGTTGTGCCTCTCGGACCCGTGGG‑CGTGCCTTTTCCT‑TGTTTG‑GGGAAAAAGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa                                                                           <  1:115565/146‑1 (MQ=255)                                                                                ccGAAACTGACCCGGGGCGCGCTGGATGTGGCCCCCGGACCCTGGGG‑CGTACCTTTTACC‑TGTAGC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                            <  1:1128558/145‑1 (MQ=255)                                                                                  ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGCAAGGATATGACCAtttt                                                                    >  2:421005/1‑151 (MQ=255)                                                                                  gaaacgaaaCGGTGGCGGGCTGGATGGGGCCCCCGCCCACGGGGG‑CGTCCTTTTTCCT‑TTTATG‑GGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                      <  1:796478/142‑1 (MQ=255)                                                                                    aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                  >  1:969536/1‑151 (MQ=255)                                                                                       tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                               >  1:489578/1‑151 (MQ=255)                                                                                       ggaccggggccgcgcgGGAGGTGCCTCTCGAAC‑CCGGGG‑GCTTCCGTTTTCCTTGTATG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                               <  2:176743/141‑1 (MQ=255)                                                                                         accgggggCGCGTTGGATGTGGCTCTCGGACCCGTGGG‑CGTACCGTTTCCT‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                             <  2:711905/145‑1 (MQ=255)                                                                                           ccGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                            >  1:590044/1‑150 (MQ=255)                                                                                            cGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                          >  1:386494/1‑151 (MQ=255)                                                                                             ggggggcgggggttgtgTCTCTGGGACGCTGGGG‑CGTCCCGTTTACT‑TGGAGG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                          <  2:1195049/137‑1 (MQ=255)                                                                                                      tGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                                >  2:273253/1‑151 (MQ=255)                                                                                                       ggttgtgtcccggagACATCTGGG‑AGTCCTGTTT‑CC‑TGTTTTGGGGAAAATTACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                                <  1:838338/136‑1 (MQ=255)                                                                                                       ggttgggcctctCGGCCCCCGGGG‑CGTCCCTTTTCCT‑TGTAGG‑GGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                   <  2:1005465/139‑1 (MQ=255)                                                                                                              gCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGt                                        >  1:472820/1‑151 (MQ=255)                                                                                                                  tcGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc                                    >  2:370127/1‑151 (MQ=255)                                                                                                                    ggaccgggggg‑CGTCCTTTTTCCC‑TTTACG‑GGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                       <  1:1160889/138‑1 (MQ=255)                                                                                                                    gacccccgggg‑CGTACTTTTTCCT‑TGTATG‑GGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                                  <  2:326044/143‑1 (MQ=255)                                                                                                                       cactgggg‑CGTCCCTTTTCCC‑TGTATG‑GGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                                   <  1:651639/142‑1 (MQ=255)                                                                                                                           tgggggttcgcGTTTCCT‑TGTATG‑GGGCAAGGACCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa                              <  2:311245/139‑1 (MQ=255)                                                                                                                            ggg‑CGTCCCTTTTCCC‑TGTAGG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                            <  2:217607/149‑1 (MQ=255)                                                                                                                             gg‑CGTACTTTTTCCT‑TGTATC‑GGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                         <  2:38546/151‑1 (MQ=255)                                                                                                                                 gTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt                      >  1:836559/1‑151 (MQ=255)                                                                                                                                   tcccttttccttGTATG‑GGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa                        <  2:1145337/137‑1 (MQ=255)                                                                                                                                   tcccttttccttGTAGC‑GGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa                        <  2:368905/137‑1 (MQ=255)                                                                                                                                       tttACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg                 >  2:88279/1‑150 (MQ=255)                                                                                                                                            c‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc            >  1:885062/1‑150 (MQ=255)                                                                                                                                                       gCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGt  >  2:938113/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGT  >  NC_000913/436368‑436665

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