Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	49.2%	63.5 / 45.1	61	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (28/3);  new base A (0/30);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 2.57e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGT  >  NC_000913/436366‑436665                                                                                                                                                         |                                                                                                                                                      cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGc                                                                                                                                                       >  2:498659/1‑151 (MQ=255)   ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGc                                                                                                                                                       <  1:857326/149‑1 (MQ=255)   ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCa                                                                                                                                                      >  1:148134/1‑150 (MQ=255)       cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATg                                                                                                                                                  >  2:1230873/1‑150 (MQ=255)        gCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGAc                                                                                                                                                >  2:58767/1‑151 (MQ=255)         ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGAcc                                                                                                                                               >  1:401655/1‑151 (MQ=255)             gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgcc                                                                                                                                            <  1:912056/150‑1 (MQ=255)              ttGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTGCGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgccg                                                                                                                                           <  1:511389/150‑1 (MQ=255)               tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGat                                                                                                                                         >  2:343099/1‑151 (MQ=255)                 aaCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGata                                                                                                                                        >  2:815516/1‑150 (MQ=255)                  aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATc                                                                                                                                      >  1:839828/1‑151 (MQ=255)                  aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATc                                                                                                                                      >  2:590697/1‑151 (MQ=255)                       aaCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAgg                                                                                                                                 >  1:183853/1‑151 (MQ=255)                              gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                          >  1:1172822/1‑151 (MQ=255)                                 gcgcggGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTTTTTCCCTGTGCCGGAACAGAGCCGTGGCGGGTTGGATGTGGCTCCCGGCCACCT‑GGGCGTCCTTTTTCCC‑TGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTg                                                                                                                           <  2:860298/147‑1 (MQ=255)                                   gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGACATCGACGGGCTTTGTTTTa                                                                                                                      >  2:967472/1‑150 (MQ=255)                                      gggCGCTCTCGGGGGGTGAAGTTTAGGAGTTGTGTTTCACGGTGCGGAAACTAACCCGTGGCGCGCTGGGTGTGGCTCGCGGCCACGG‑GGGGGTCCCTTTTCCT‑TGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                                  <  2:691494/151‑1 (MQ=255)                                      gggCGCTCTCGGGGGGTGAAGATTACGAGTTGTTTTTCCCTGGGCCGAAACTGACCCGGGGCGGGTTGGGTGTGGCTCTCGGCCCCCG‑GGGCGTCCTGTTTTCT‑TGTACGGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTa                                                                                                                      <  1:603814/147‑1 (MQ=255)                                      gggCGCTCTCGGGGGGTAAAGATTACGATTTGTTTTTCACTGTGCCGGAACTGACCCGTGGCGCGTTGGATGTGGCTCTCGGACCCGG‑GGGCGTCCCTTTTCCC‑TGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                                  <  1:427778/151‑1 (MQ=255)                                       ggCGTTCTCGGGGGGTAAAAATTGCGAGTTGTTTTTCGCTGTGCGGGATCAGAACGGTGGCGCGTTGGGTGTGCCTCTC‑GACCCCTGGGGCTTCCCTTTTCCT‑TGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                                 <  1:507902/151‑1 (MQ=255)                                       ggCGCTCTCGGGGGGTGAAGATAAGGATTTGTGTTTCCCTGGGCCGGAACGGACCCGGGCCGGGGTGGGGGGGGCTCCGGGCCCCGT‑GGGCGTCCCTTTTCCC‑TTTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                                 <  2:1025593/151‑1 (MQ=255)                                              tctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                          >  2:355479/1‑151 (MQ=255)                                                  gCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                                                      >  2:1054197/1‑151 (MQ=255)                                                   gggTAAAGATTAGGAGTTGTTTTTCCCTGTGCGGGAACTGACCCGTGGCGCGTTGGATGTGGCTCTGGGACCCTG‑GGGCGTACTTTTTCCC‑TGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                     <  2:1132008/150‑1 (MQ=255)                                                    ggTGAAAATTACGATTTGTGTTTCATGGTGCCGGAACTGACCGGTGGGGCGCTGGATGTGGCTCTCGAACACTG‑GGGCGTCCCTTTTCCT‑TGTTTCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                                                      <  1:1010515/149‑1 (MQ=255)                                                    ggTAAAGATTAGGGGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGGGGGTTGTGGCTCTCGGACACGG‑GGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                    <  1:657797/151‑1 (MQ=255)                                                           aTTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                             >  1:1068920/1‑151 (MQ=255)                                                             tACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                             >  2:78433/1‑149 (MQ=255)                                                                    tgttTTTCACTGTGCCGAAACTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACCCCG‑GGGCGTCCCTTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                     <  1:763515/150‑1 (MQ=255)                                                                      tgtTTCCCTGTGCGGGAATTGACCCGGGGCGCGCTGGATGTGGCTCCCGGACCCGG‑GGGCGTCCCTTTTACC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                     <  2:325889/148‑1 (MQ=255)                                                                      tgtTTCACTGTGCCGAAACGGACCCGGGGCGCTTTGGATGTGGCTCGCGGACCCCT‑GGGCGTCCCGTTTCCC‑TTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                    <  2:680741/149‑1 (MQ=255)                                                                         ttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                               >  1:1081465/1‑151 (MQ=255)                                                                          tCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                               >  2:437169/1‑150 (MQ=255)                                                                           ctttgtgCCGGAACGGACCTGGGGCGCGTGGGTTGTGGCTCCGGGCCCCGG‑GGGCGTTCCGTTTCCT‑TTTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                               <  1:146044/146‑1 (MQ=255)                                                                              tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                          >  2:648826/1‑151 (MQ=255)                                                                                 gCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                       >  1:1273603/1‑151 (MQ=255)                                                                                  cggaaaCTGACCGGGGGCGCGTTGGATGGGGCTCGCGGACCCCG‑GGGCGTCCCGTTTACT‑TTTTTCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                         <  1:319054/144‑1 (MQ=255)                                                                                  ccGAAACTGACCCGGGGGGCGTTGGATGTGGCTCTCGGCCCCGT‑GGGCGTCCCTTTTACC‑TGTCGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                       <  2:1104448/150‑1 (MQ=255)                                                                                    ggAACTGACCCGTGGCGCGCTGGTTGTGGCTCTCGGCCCCCG‑GGGCGTCCCTTTTACC‑TGTTGCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                    <  1:496603/151‑1 (MQ=255)                                                                                     gAACTGACCCGTGGCGCGGTGGTTGTGGCTCTCGGACCCTG‑GGGCGTACCTTTTCCT‑TTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                     <  1:694723/149‑1 (MQ=255)                                                                                           acccGTGGCGCGCGGGATGTTCCTCTCGGACCCGG‑GGGGGTCCTTTTTCCC‑TGTACCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                                                <  1:981301/146‑1 (MQ=255)                                                                                           aaCCGGGGCGCGCTGGTTGTGGCTCTCGGCCCCCG‑GGGCGTACCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                             <  2:120108/151‑1 (MQ=255)                                                                                             ccGGGGCGCGTGGGATGTGGCTCTCGGACACGT‑GGGCGTACCTTTTCCC‑TGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                           <  2:69308/151‑1 (MQ=255)                                                                                               ggggCGGGTTGGATGGGGCTCCCGGACCCGT‑GGGCGTACTTTTTCCC‑TTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                         <  2:560056/149‑1 (MQ=255)                                                                                               ggggCGCGGGGGATGTGTCTCCGGGACGCCG‑GGGCGTACTTTTTCCT‑TTTAGGGGAAAAAAGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                         <  2:60093/149‑1 (MQ=255)                                                                                                 ggCGCGGGGGATGGGGCTCCCGGACCCGG‑GGGCGTCCCTTTTCCT‑TTTTGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                        <  1:723981/150‑1 (MQ=255)                                                                                                  gggggttGGTTGTGCCTCCCGGACACCT‑GGTCTTACTTTTTCCC‑TGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                           <  2:148132/140‑1 (MQ=255)                                                                                                    gcgcTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                        >  2:554214/1‑147 (MQ=255)                                                                                                      gcTGGATGTGGCTCTCGGACCCTG‑GGGCGTCCCTTTTCCC‑TGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt                                                    <  1:173648/149‑1 (MQ=255)                                                                                                      gcTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                                                  >  1:960869/1‑151 (MQ=255)                                                                                                              tgGCTCTC‑GACCCCC‑GGGGGTTCCGTTTTCCTTATAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgc                                               <  1:286547/146‑1 (MQ=255)                                                                                                                  tctcGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                      >  2:1224657/1‑151 (MQ=255)                                                                                                                     cGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                                   >  1:990617/1‑151 (MQ=255)                                                                                                                       ggccccgggggCGTCCCTTTTCCT‑TGTACGGGAAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                                  <  2:772284/143‑1 (MQ=255)                                                                                                                             t‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                            >  1:227980/1‑150 (MQ=255)                                                                                                                                       ctttttgcctatacggggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                    <  1:261794/135‑1 (MQ=255)                                                                                                                                        ttttCCC‑TTTATGGGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                    <  2:468798/147‑1 (MQ=255)                                                                                                                                         tttACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg                >  2:5243/1‑151 (MQ=255)                                                                                                                                             cc‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc            >  1:1200138/1‑151 (MQ=255)                                                                                                                                                      gggCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt    >  1:205318/1‑151 (MQ=255)                                                                                                                                                      gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc               <  2:118038/140‑1 (MQ=255)                                                                                                                                                      gggAAAATCACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt       <  1:134056/148‑1 (MQ=255)                                                                                                                                                        gCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGt  >  1:348659/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                      CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGT  >  NC_000913/436366‑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: ‑