Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	606,619	0	T	C	40.3%	84.7 / 21.9	62	G255G (GGA→GGG)	ybdK	weak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base T (30/7);  new base C (0/25);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 2.38e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAA‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCG  >  NC_000913/606470‑606759                                                                                                                                                            |                                                                                                                                             cgacgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAA‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAAt                                                                                                                                              <  2:1207804/150‑1 (MQ=255)     gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAA‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGc                                                                                                                                         >  1:119316/1‑151 (MQ=255)                   ggATCGGTGATGACGCCTTCAA‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGt                                                                                                                           >  1:206937/1‑151 (MQ=255)                      tggtttatgcCGCCTCCAA‑GCCCATGGCGCCGGGCGGGAAAA‑GGTTTGAATTTAAACACAAGATAATTTTTTTCTT‑TTTGTTAAAG‑G‑GGGCCTCCGGGAAGAACCG‑AGTGGGGGGTAGCATAAATTCACCCCGCCATATTTACTGCGTg                                                                                                                          <  1:334916/143‑1 (MQ=255)                       cGGTGATGACGCCTTCAA‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCt                                                                                                                       >  1:457828/1‑151 (MQ=255)                            aggccGCCTTCAA‑CCCCATAGCGACAGGCTTGGAAA‑GGGTTAATTTTATACAGGAGTTATTTTTTTTCTT‑GATTTTTAAG‑C‑GGGCTTTCCGTCAAAACCG‑GGTGGGCGGTAGCTTAATTTAATCCCGCCATATTTACTGCGTGGCTAAgg                                                                                                                   <  2:1043963/146‑1 (MQ=255)                                 cccTTCAC‑CCCAATGGCGAC‑GGGCTGGAAAGGGGTTAATTTTATACAGCAGTTAATCTTTTTCTT‑G‑TTTTTAAA‑CGGGGCTTTCCTTCAGTACCC‑GGGGGGGGGTAGCTTAATTTAACCCCGCCATATTTACTGCGTGGCTAAGGGt                                                                                                                 <  1:110578/146‑1 (MQ=255)                                  ccTTCAA‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACg                                                                                                            >  2:131486/1‑151 (MQ=255)                                   cTTCAACGCCCAGAGCGAAGGGCCTGAAAT‑CGGTTAATTTTAAAAAGAAGGTAATTTTTTCCCT‑GATGTAAAAG‑C‑GGGCTTGCTGTCATTAACG‑GGGGGGGGGGACCTTCAATAAATCCCGCCATATTTACTGCGTGGCTAAGGGt                                                                                                                 <  1:406566/146‑1 (MQ=255)                                        a‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTAt                                                                                                      >  1:1071184/1‑151 (MQ=255)                                           cccagggCGACGGGCCGGGAAA‑GGTTTAATTTTATACGGCA‑GTTATTTTTTTCCTGGATTTTTAAG‑G‑GGGGTTTCCGTCAGTACCC‑GGTGGGGGGTACCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATcc                                                                                                    <  2:961275/145‑1 (MQ=255)                                              aaaGCGACGGGCCGGAAAA‑GGTTTAAATTTATACAGCAGTTAATTTTTTCCTT‑GATGTTAAAC‑G‑GGGCGTTCCGTAAGAACCC‑GGGGGGCGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCAt                                                                                                  <  2:116582/148‑1 (MQ=255)                                              aTAGCGACAGGCCTGGAAA‑GGTTTAAATTTATACAGCAGTTATTTTTTTTCCG‑TATTTTTAAC‑C‑GGGCGTTCCGTCAGTACCC‑GGGGGGCGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATc                                                                                                 <  1:628820/151‑1 (MQ=255)                                                   gacgggctgggAAA‑CGGTTGATTTTATACAGCAGGTAATTTTTTTCTT‑TTTGTTTAAA‑C‑GGGCGTCCGGTCAGTAACG‑AGGGGGCGTTAGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCAc                                                                                               <  2:225214/139‑1 (MQ=255)                                                       ggCGTGAAAG‑CGTTTAAATTTATACAGCAGGTATTTTTTTTCCT‑GATTTTAAAA‑G‑GGGCTTCCCTTCAGTACCC‑GGGGGGGGGTAGCCTGAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCAc                                                                                               <  1:535714/144‑1 (MQ=255)                                                       ggCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAAcc                                                                                        >  2:59773/1‑151 (MQ=255)                                                       ggCCAGGAAA‑CGGTTAAATTAATCCACAAGGTATTTTTTTTCCT‑TTTGTTAAAT‑C‑GGGCGTTCCGTAAGAAACG‑GGTGGGGGTTAGCTTGATTAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAAcc                                                                                        <  2:866182/151‑1 (MQ=255)                                                        gCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCt                                                                                       >  1:7509/1‑151 (MQ=255)                                                           tGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTcc                                                                                     >  1:196751/1‑150 (MQ=255)                                                            ggaaaggggTTGAATTTATACAGCAGTTATTTTTTTTCCT‑GATGTTTAAA‑G‑GGGCGTCCCGTCAGTACCG‑GGGGGGGGGTAGCTTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTc                                                                                      <  2:1096972/142‑1 (MQ=255)                                                            ggAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCTATCACCCGAACCTCCAc                                                                                   >  2:1209833/1‑151 (MQ=255)                                                             gAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCAcc                                                                                  >  2:523734/1‑151 (MQ=255)                                                              aaa‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCg                                                                                 >  2:529337/1‑151 (MQ=255)                                                               aagggtttaattttATACAGCAGTTATTTTTTTTCCT‑GTTTTTTAAC‑G‑GGGCGTTCCGTAAGTAACG‑GGGGGGGGGTAGCCTGAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCAcc                                                                                  <  2:1187933/139‑1 (MQ=255)                                                                       gAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCaaa                                                                          >  2:865737/1‑150 (MQ=255)                                                                           ttATACAGCGGGTAATTTTTTCCGT‑TTTGTTTAAC‑C‑GGTCTTCCCGCCGGAACCC‑AGGGGGCGGTCCCTTAAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCaaa                                                                          <  1:1256851/146‑1 (MQ=255)                                                                                  gcagTTATTTTTTTCCCT‑G‑TTTTTAAA‑CGGGGCTTCCCTTCAGAACCC‑GGGGGGGGGTACCCTAAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGAc                                                                  <  1:1150081/147‑1 (MQ=255)                                                                                  gcagTTAATTTTTTTCCT‑GATTTTTAAC‑G‑GGGCTTTCCGTCAGTAACC‑GGTGGGGGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAgg                                                                    <  2:1211337/145‑1 (MQ=255)                                                                                    agatattttttttCCT‑G‑TTTTTAAA‑CGGGGCTTCCCTTAAGAACCC‑GGGGGGGGGTACCCTAAATAACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACt                                                                 <  1:350738/138‑1 (MQ=255)                                                                                    agGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGc                                                            >  2:50122/1‑151 (MQ=255)                                                                                       tattttttttCCT‑GATGTTAAAC‑C‑GGGCGTCCCGTCAGTAACG‑AGGGGGGGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTgg                                                               <  1:789990/140‑1 (MQ=255)                                                                                        attttttttCCT‑GTTTTTAAAG‑C‑GGGCTTTCCGTCAGAACCG‑GGGGGGGGGTAGCCTAATTTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGa                                                          <  2:128170/145‑1 (MQ=255)                                                                                          tttttttCGT‑TATGTTAAAG‑G‑GGGCTTCCCGTCAGAAGCC‑GGGGGGCGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                           <  1:860250/144‑1 (MQ=255)                                                                                          tttttttCGT‑GTTTTTCAAA‑C‑GCGCTTCCCGAAAGAACCC‑GGGGGGCGGTACCCTGATTTATCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAata                                                       <  1:1077125/148‑1 (MQ=255)                                                                                            ttttcCCT‑TATGTTTAAA‑C‑GGGTTTTCCG‑AAATAACCGGGGGGGCGGTACCTTAATTTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTgggg                                                             <  2:813380/142‑1 (MQ=255)                                                                                             ttttCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATccc                                                   >  1:650997/1‑151 (MQ=255)                                                                                               ttCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATccc                                                   >  2:559323/1‑149 (MQ=255)                                                                                                  cT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              >  1:622104/1‑151 (MQ=255)                                                                                                   t‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              >  1:685564/1‑150 (MQ=255)                                                                                                      aTGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              >  1:802640/1‑148 (MQ=255)                                                                                                       tGTTTAAC‑C‑GGGCGTCCCGTAAATACCG‑GGGGGGAGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTg                                               <  1:1210757/146‑1 (MQ=255)                                                                                                       tGTTAAAC‑C‑GGGCTTTCCGT‑AATAACCTGGGGGGGGGTTGCCTAATTAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAg                                                 <  2:989297/144‑1 (MQ=255)                                                                                                        ttttAAAGC‑GGGCGTCCCGTCAGAACCG‑GGGGGGCGGTAGCCAAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAg                                                 <  2:779781/143‑1 (MQ=255)                                                                                                        gTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              >  1:786307/1‑146 (MQ=255)                                                                                                        gTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              >  2:1031843/1‑146 (MQ=255)                                                                                                         ttttaaagggggCTTTCCGTCAGAAACG‑GGGGGGCGGTAGCTTGAATTATTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              <  2:504190/138‑1 (MQ=255)                                                                                                         tttAAA‑G‑GGGCGTTCCGTCAGTAATC‑GGGGGGCGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAg                                                 <  2:344620/142‑1 (MQ=255)                                                                                                                   ggCGTTCCGTCAGTAACC‑GGGGGGGGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAAt                                   <  2:119316/148‑1 (MQ=255)                                                                                                                     cGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGt                              >  1:779351/1‑151 (MQ=255)                                                                                                                      gTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTc                             >  2:902398/1‑151 (MQ=255)                                                                                                                       ttCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCg                            >  2:1255081/1‑151 (MQ=255)                                                                                                                         ccGTCAGAACCC‑GGGGGGCGGTACCATAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGat                          <  1:250645/151‑1 (MQ=255)                                                                                                                          cggtaagaacccggggggCGGTGGCTTAATTTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGa                           <  2:321692/135‑1 (MQ=255)                                                                                                                          cGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatc                         >  1:67591/1‑151 (MQ=255)                                                                                                                           gTCAGTAACC‑AGGGGGGGGTAGCCTGAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGat                          <  2:1066982/149‑1 (MQ=255)                                                                                                                             cAGTACCC‑GGTGGGGGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatcatc                      <  2:1088370/151‑1 (MQ=255)                                                                                                                             cAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTTCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatcatc                      >  1:504914/1‑151 (MQ=255)                                                                                                                              aGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCg                     >  2:482446/1‑151 (MQ=255)                                                                                                                                  aCC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgtggt                 >  2:1142195/1‑151 (MQ=255)                                                                                                                                         gggCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACt           <  2:100864/151‑1 (MQ=255)                                                                                                                                                 gCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGc   >  2:83888/1‑151 (MQ=255)                                                                                                                                                 gCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGc   >  1:279722/1‑151 (MQ=255)                                                                                                                                                   cTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCg  >  1:1224536/1‑150 (MQ=255)                                                                                                                                                            |                                                                                                                                             CGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAA‑GCCCATAGCGACAGGCCTGGAAA‑CGGTTGAATTTATACAGCAGGTAATCTTTTTCCT‑GATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACC‑AGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCG  >  NC_000913/606470‑606759

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