Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	48.4%	69.0 / 42.8	62	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (26/6);  new base A (0/30);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 4.45e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTA  >  NC_000913/436368‑436657                                                                                                                                                        |                                                                                                                                              ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGc                                                                                                                                               <  1:369907/149‑1 (MQ=255)  tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCaa                                                                                                                                             >  1:519664/1‑150 (MQ=255)     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATg                                                                                                                                          <  1:350437/150‑1 (MQ=255)           gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgccg                                                                                                                                   >  2:23545/1‑151 (MQ=255)           gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgcc                                                                                                                                    <  2:248234/150‑1 (MQ=255)             tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGat                                                                                                                                 >  1:501499/1‑151 (MQ=255)              gAACCGGAACAGGCGCTGCACTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGat                                                                                                                                 >  2:814606/1‑150 (MQ=255)                  cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCATACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCg                                                                                                                             >  2:526762/1‑150 (MQ=255)                   ggAACGGGCGCTGCGCGGGGCGCTCTCGGGCGGTGAAGATTACGAGTTGTTTTTCCCTGTGCCGGAACTGAC‑CC‑GGGGCGCGTTGGATGTGGCTCTCGGACACCGGGGCGTACCTTTTCCC‑TGTAGCGGGAAAATGACCGCCGATATCg                                                                                                                             <  2:677863/149‑1 (MQ=255)                      aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAggg                                                                                                                        >  1:205842/1‑151 (MQ=255)                         ggcgcgggggggggCGCTCGCGGGGGGTAAAGATTAGGGGTTGTGTTTCACGGGGCCGGACCAAAC‑CG‑GGGGCGCGAGGGAGGGGGCGCTCGGACACGGGGGGGGCCATTTTACC‑GGAAGGGGGAAAAGGACCGCCGATATCGAAGGGc                                                                                                                       <  2:735521/138‑1 (MQ=255)                              tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                >  1:616973/1‑151 (MQ=255)                                cggtgGGCGCTCTCTGGGGGTGAAGATTGCGGTTTGTGTTTCACTGTGCCGAAACTAAA‑CC‑GGGGCGCGCGGGATGTGGCTCTCGGACCCGGGGGGGTCCCTTTTACC‑TGTATGGGGCAAAGGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                 <  1:635585/145‑1 (MQ=255)                                  ctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                            >  2:9909/1‑151 (MQ=255)                                  cggggCGTTCTCGGGCGGTAAAGATTCGGTTTTGTGTTCCCCTGT‑CCCGAACCGAATCG‑GTGGCGGGTGGGGTTTGTCTCCCGCAAACGGGGGGGGACCTTTTACT‑TTTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                <  2:422974/145‑1 (MQ=255)                                     ggCGCTCTCGGGGGGTGAAGTTTACGGGTTGTGTTTCACTGGGCCGGACCTGAC‑CC‑GGGGCGCGTTGGATGTGGCTCCCGGACCCTGGGGCGTACCGTTTTCT‑TTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                         <  2:729347/151‑1 (MQ=255)                                        gcTCTTTGGGGGGAAAGTTTACGTTTTGTTTTTCAGTGT‑CCCGAACTGAA‑CCGGTGGCGCGTTGGTTGTGGCTCCCGCACACGTGGGTTTCCTTTTTCCC‑TTTAGGGGACAATTAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGa                                                                                                       <  2:728103/150‑1 (MQ=255)                                         ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                     >  2:869434/1‑151 (MQ=255)                                         ctctcgGGCGGTGAAGTTTACGGGTTGTGTTTCCCTGGGCGGGAACTGAC‑CG‑GGGGCGCGCGGGATGTGGCTCTCGGCCACCGGGGCGTCCCTTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGa                                                                                                       <  1:145839/149‑1 (MQ=255)                                            tttGGGGGTGAAGATGACGTGTTGTTTTTCACGGGGCCGGAACTGAA‑CC‑GGGGGGCTGTGGGTGTGGCTGCCGGGCTCCTGGGGGTTCGGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                                      <  1:118515/145‑1 (MQ=255)                                            tctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                   >  2:79998/1‑150 (MQ=255)                                              tGGGGGTGAAGATTCCGATTTGTTTTTCACTGTGCCGGAACTGAA‑CG‑GTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                     <  1:366792/146‑1 (MQ=255)                                               gggggTGAAGATTACGGGTTGTGTTTCACTGTGCCGGAACTGAC‑CG‑GGGGCGCGCTGGATGTGGCTCCCGGCCCCGGGGGCGTACTTTTTCCT‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                    <  2:871324/143‑1 (MQ=255)                                                     gAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                          >  1:188985/1‑150 (MQ=255)                                                      aaGATTGCGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGGTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                         <  1:364967/150‑1 (MQ=255)                                                      aaGATTAGGGGTTGTTTTTCCCTGTGCGGGAACTGAC‑CC‑GGGGCGCGCTGGGTGTGGCTCGCGGACACCGGGGCGTCCCTTTTCCC‑TGTGTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                          <  2:951431/149‑1 (MQ=255)                                                      aaGATTACGGGTTGTTTTTCACTGGGCCGGAACTGAC‑CC‑GGGGCGCGCGGGATGGGGCTCTCGGACCCCGGGGCGTCCTTTTTACC‑TGTGTGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                          <  2:92531/149‑1 (MQ=255)                                                          ttACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                     >  1:473143/1‑150 (MQ=255)                                                            aCGAGTTGTTTTTCACTGTGCGGGAACTGAC‑CC‑GGGGCGCGTTGGATGGGGCTCTCGGCCCCTTGGGCGTCCCTTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                                  <  1:811749/151‑1 (MQ=255)                                                             cGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                                 >  1:844162/1‑151 (MQ=255)                                                                 ttttgtTCCATTGTGCCGAAACGGAC‑GC‑GGGGCGGGTTGGATGTGCCCCGCGGCCCCGGGGGCGTACGTTTTACT‑TGTGTGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                   <  1:94319/142‑1 (MQ=255)                                                                  tgtgTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                             >  1:121164/1‑150 (MQ=255)                                                                   gtgtTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                            >  2:831173/1‑150 (MQ=255)                                                                   gtgtTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                            >  2:343200/1‑150 (MQ=255)                                                                   gtgtTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                            >  1:873247/1‑150 (MQ=255)                                                                      tttCTCGGTGCCGGAACTGAC‑CC‑GGGGCGCGTTGGATGTGGCTCTCGGCCCCCGGGGCGTCCCTTTTACC‑TGTATCGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                         <  2:367653/150‑1 (MQ=255)                                                                        tCACTGTGCCGGAACTAAC‑CC‑GGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTACC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                       <  1:989107/150‑1 (MQ=255)                                                                         cccTGTGCCGGAACGGAA‑CG‑GTGGCGCGTTGGTTGGGGCTCTCGGACCCCTGGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                       <  1:414618/147‑1 (MQ=255)                                                                         cACGGTGCCGGAACTGAC‑CG‑GTGGCGCGCTGGTTGTGGCTCTCGGCCGCGTGGGCGTCCCTTTTCCC‑TGTACCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                       <  1:137047/149‑1 (MQ=255)                                                                          aCTGTGCGGAAAGTGAA‑CG‑GGGGCGCGTTGGATGTGGCTCTCGGACGCTTGGGCGTCCCTTTTCCT‑TGTATCGGACAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                      <  1:415586/149‑1 (MQ=255)                                                                             gtgCCGGAACTGAC‑CG‑GTGGGGCGCTGGTTGGGGCTCCCGGACCCCGGGGCGTCCCTTTTACT‑TGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                  <  1:173337/150‑1 (MQ=255)                                                                               gCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATCCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                               >  2:346260/1‑151 (MQ=255)                                                                                 ccgcactggaccc‑GTGGCGCGGGGGTTGTGTCTCTCGACCTCGGGGGCGTCCTTTTTGCC‑GTTGTGGGGAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                      <  2:279405/132‑1 (MQ=255)                                                                                    aaCTGAC‑CC‑GGGGCGCGTGGGTTGTGCCTCTCGGCCCCCTGGGCGTACTTTTTCCC‑TGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                          <  2:584755/151‑1 (MQ=255)                                                                                    aaCTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                          >  1:349649/1‑151 (MQ=255)                                                                                           accgggggggcgggggATGTGGCCCCCGGACCCGGGGGGGTCCTTTTTCCT‑TGTTCGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                       <  2:574970/135‑1 (MQ=255)                                                                                                      gcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                                          >  1:293910/1‑151 (MQ=255)                                                                                                         gtgtggggCTCTCGGACCCCTGGGCGTCCTTTTTACC‑TGTACCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                           <  2:284771/141‑1 (MQ=255)                                                                                                             gtgCCTCTGGGACCCCTGGGCGTCCCTTTTCCT‑TTTACGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                        <  1:804297/146‑1 (MQ=255)                                                                                                              tgGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                   >  2:290006/1‑150 (MQ=255)                                                                                                                      ggccgctggggCGTCCGTTTTCCC‑TGTAGGGGGAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                              <  2:750234/139‑1 (MQ=255)                                                                                                                      ggcccccggggCGTCCCGTTTCCC‑TGTCGGGGGAAAAGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                              <  1:467576/139‑1 (MQ=255)                                                                                                                       gACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                          >  2:633549/1‑150 (MQ=255)                                                                                                                           ccTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                         >  1:147002/1‑147 (MQ=255)                                                                                                                             ggggctcacggTTTCCC‑GGTATGGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa                      <  1:594151/138‑1 (MQ=255)                                                                                                                             ggggCGT‑CCGTTTTCCTTGTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt                     <  1:23546/148‑1 (MQ=255)                                                                                                                               ggCGTCCCTTTTCCC‑TTTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                 <  1:830855/151‑1 (MQ=255)                                                                                                                                     ccGTTTCCT‑TGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg           <  1:162200/151‑1 (MQ=255)                                                                                                                                     ccGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg           >  2:331772/1‑151 (MQ=255)                                                                                                                                      attttGCC‑TGTTTCGGAAAATTCACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGt          <  1:33293/149‑1 (MQ=255)                                                                                                                                         ttACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc       >  2:551159/1‑151 (MQ=255)                                                                                                                                               tGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt   >  1:969463/1‑150 (MQ=255)                                                                                                                                                 tATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTa  <  1:52749/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                              TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTA  >  NC_000913/436368‑436657

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