Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,748	G→A	S404S (AGC→AGT)	yaaJ ←	putative transporter YaaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,748	0	G	A	96.9%	180.7 / NA	66	S404S (AGC→AGT)	yaaJ	putative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (29/33);  minor base C (1/1);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAAT‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCC  >  NC_000913/6511‑6990                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    gTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTgc                                                                                                                                                                                                                                        >  1:60162/1‑251 (MQ=255)   tGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTgcc                                                                                                                                                                                                                                       <  2:680376/250‑1 (MQ=255)          tAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGATGCCGCCGATg                                                                                                                                                                                                                                >  1:189973/1‑250 (MQ=255)            gTTGTGATTCCTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGAcc                                                                                                                                                                                                                             <  2:182919/251‑1 (MQ=255)            gTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGAcc                                                                                                                                                                                                                             <  1:214282/251‑1 (MQ=255)                     ccacattcGAAACATCACCCCCCGCGTCCGGATAAAGCTGGCGGTCGATACCCGTATACCGCAAC‑GGAACAAACACGGGGCGCACGCCGCGTTTACGCTGGCGTAGATAATCTCTGGCTATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGTTAAGAGCCAGGCAGGCCA‑TGATTATATCTTCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTCGCGaa                                                                                                                                                                                                                     <  2:506470/243‑1 (MQ=255)                     aCTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAAT‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCCCATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGa                                                                                                                                                                                                                      >  1:498801/1‑249 (MQ=255)                          tGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGACAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGAAAAATTGGTAATAGCCAAGCAGGCCA‑TTATGATATCTGCCAGTTGCCACAACAACGGGAGTCTTAGGAAGGGGACGCCGATGGCCGTTGCGGAggcgcc                                                                                                                                                                                                               >  1:11563/1‑247 (MQ=255)                              aaaCATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAAAACCATGAAGGCCA‑TTATTATATCTGCCAGTTGCCACAACAGCGGGAGACTTAGAAAGGTGCCGCAGCTGACCGTTGCGCAgggggggaa                                                                                                                                                                                                            >  2:612871/1‑243 (MQ=255)                                 catcatCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATTTGCCAGTTGCCACATCAGACGAAGACTTAGCAAGGTGCCGCCGATGACCGTTTCGAAGGGGCAGA‑TCCGc                                                                                                                                                                                                       >  1:214076/1‑251 (MQ=255)                                  atcaccCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCa                                                                                                                                                                                                      <  1:347417/251‑1 (MQ=255)                                         cACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGGTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCa                                                                                                                                                                                               <  1:418313/251‑1 (MQ=255)                                             cgTCCGAAGAAAGCTGGCGGCCGAATTCCGGATAACG‑AATCGGATCAAACACCGGTCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAAGGGTATGAACCTCAGTCGAGAGCAGTAAAATGTCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGACGCCGATGACCGTTGCGAAGGTGCAGT‑TCCGCAAACACCAGATc                                                                                                                                                                                           <  2:615531/251‑1 (MQ=255)                                              gTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCC‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCCAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGA‑TCCGCAAACAcc                                                                                                                                                                                                >  1:89096/1‑245 (MQ=255)                                                      aaaGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCTCACGCAGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCGTTAggg                                                                                                                                                                                  <  2:569042/251‑1 (MQ=255)                                                        aGCTGGCGGCCGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGgtt                                                                                                                                                                                <  1:333223/251‑1 (MQ=255)                                                             gtgtgcgttatCCGGATAACGAAAA‑TGATCAAACACCGGGGGCTCGCCGAGTTTACGCTGGCGTAGATAATCATTGGAAAGGTTATGAACCACAGGTGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGCCTTAGCAAGGTGCCGGCGATGACCGTTGCGAAGTTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGTTTGTTCa                                                                                                                                                                           <  2:443623/243‑1 (MQ=255)                                                                 ccGATATCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCCAATTGGTAATAACCCTGCAGGCCA‑ATATTATATCTGCCAGTTGCCAACTACGCGGGAGACTTAGACAGGTGCCGCCCATGAACGTTGCGAAGGGGGAGA‑GCCGCCAACAACAGACCGACTTAGGGGTGTTCAGGCg                                                                                                                                                                       >  2:32345/1‑251 (MQ=255)                                                                      atCCGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGAc                                                                                                                                                                                                                              <  1:21308/192‑1 (MQ=255)                                                                        ccGGATAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGGCA‑TTATGATATCTGCCAGTTGCCACATCAGCAGAAGACTTAGCCAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGA‑TACGCAAACACCACATCGCCTTAGGggtgttttggcgtgcaaa                                                                                                                                                                 >  2:306393/1‑238 (MQ=255)                                                                          ggATAACGCAAC‑GTATCGAACCCCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGATCCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATCGCCATGCAGGGCA‑TTAAGATATCTGCCAGGTGCCACAGAAGCGGGAGACTTAGCAAGGGTCCGCCGGTGGCCGTTGCGTCGGCGCCGAACCCGCAA‑CACCAGGTCAGTTTAGGGTTGTTCGGGCggaactaga                                                                                                                                                               >  2:91786/1‑242 (MQ=255)                                                                            aaaaggCAAC‑GGAGCAAACCCCGGGGGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAAACACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAATACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCTCTTGAGGGTTGTGCAGGCGTAAAAagaaga                                                                                                                                                            <  2:80741/246‑1 (MQ=255)                                                                            aTAACGCAAC‑GGATCAAACCCCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAagaaga                                                                                                                                                            <  1:501103/251‑1 (MQ=255)                                                                            aTAACGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAagaaga                                                                                                                                                            <  1:248842/251‑1 (MQ=255)                                                                               aCGCAAC‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGAt                                                                                                                                                         <  1:612871/251‑1 (MQ=255)                                                                                 cgcaacggtccAAACTCCGGGCGAACGTCGAGTTTACGCTGGCGGAGATAATCCCTGGCAATGGTATGAACCACAGCCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATATCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCCCATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGTGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGAtt                                                                                                                                                        <  2:452925/241‑1 (MQ=255)                                                                                   aaacgacttaaacaccgtgtccACGACGAGTTATCGCGGGCGTAGATAAACACTGGAAATGTGATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTTCCAGTTGCCACAGCAGCGGAAGACTTAGCAAGGTGCCGCTGATGACCGTTGCGAAGGTTCAGA‑TCCGCAAACTCCAGATCGCTTTAGGGTTGTTCAGTCGTAAAAAGAAGAGATTGt                                                                                                                                                      <  2:396470/230‑1 (MQ=255)                                                                                      cggATCAAACCCCGGGCCCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGA‑TCCGCAAACTCCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGtttt                                                                                                                                                   <  2:89096/250‑1 (MQ=255)                                                                                        gATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCg                                                                                                                                                 <  1:596861/251‑1 (MQ=255)                                                                                              acacCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATAGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCATGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGGGCCGCCGATGCCCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGGTGTTCAGGCGTAAAAAGACGAGATTTTTTTCGGcgta                                                                                                                                            >  1:236007/1‑247 (MQ=255)                                                                                              acacAGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATATGTAATAGCCATGCAGGCCAATTAGGA‑ATCTGCCAGTTGACACATCAGCGGAAGAATTAGCAAGGTGACGCCGATAAACGGTGCGAAGGGGCAGA‑ACCGCAAAAACCCGCACGCATTAGGGTAGTTCAAGAGTAAAAAGAAGAGGTTGTTTTctggg                                                                                                                                              >  2:189124/1‑244 (MQ=255)                                                                                                   gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATATGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGGTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGGTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGt                                                                                                                                       >  1:92124/1‑250 (MQ=255)                                                                                                                        tGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCCCATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGGAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGa                                                                                                                  >  1:496598/1‑250 (MQ=255)                                                                                                                        tGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAACCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑ACCGCAAACACCAGATCGCTTTAGGGGTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTgc                                                                                                                  >  1:538214/1‑249 (MQ=255)                                                                                                                          gCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTTTGATATCTGCCCGTTTCCACATCAGCGGAAGCCTTAACAAGGTGCCGCCGCTGCCCGTTGCGCAGGTGCAGA‑TCCGCAAACACCAGATCCCTTTAGGGGTGTTTAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTgccgg                                                                                                               >  1:398087/1‑247 (MQ=255)                                                                                                                                   aTCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGaa                                                                                                      <  1:99354/251‑1 (MQ=255)                                                                                                                                     ccttGGCAATGGTAGGAACCACAGGAGAGAGCAGTAAAATGGCGGGCAAATTTGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAAGCAAACAGAATa                                                                                                    <  2:109994/248‑1 (MQ=255)                                                                                                                                         ggCAATGGTATGAACCACAGGCGATAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATTCAGGCCA‑TTATGATATCTTCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTCCAGTCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTGGGCAACGATGGAGCTGAAGGCAAACAGAATAAc                                                                                                  <  1:561601/249‑1 (MQ=255)                                                                                                                                                      aCCACAGGCGAGAGCAGGAAAATGGCGGACCAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCAGCCAGTTGCCACATAAGCGGAAGACTTAGCAAGGGGCCGCCGATGACCGTTGGGAAGGGGCAGG‑TCCGCAAAAACCAGATCGCTTTAGGGTTGTTCAGGCGTGAAAAGAAGAGATTGTTTTCGGCGTACATGTAGTTGGCAAAGCTGGAGCTGAAGGAAAAAAGCATAAACACGAGGGTAACaa                                                                                    >  2:109751/1‑250 (MQ=255)                                                                                                                                                        cacaGGCGAGATCAGTAAAATGGCGGTTAAAGTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGGTGCCACATCAGCGGAAGACTTAGCAAGGTGTCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTTGTTGGCAACGATGGAGCTGAAGGCAAACATAATAACCACGAGGGTAACAAACt                                                                                 <  2:392602/251‑1 (MQ=255)                                                                                                                                                         acaGGAGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGGCATGCAGGCCA‑GTATGATATCTGCCAGTTGCCACAGCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTc                                                                                <  2:39106/251‑1 (MQ=255)                                                                                                                                                              cGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGgtt                                                                                                                                                                                <  1:664421/150‑1 (MQ=255)                                                                                                                                                               gagaGCAGTAAAATGGCGGGCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGAGGAAGACTTAGCAAGGGGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACCCCAGATCGCTTTAGGGGTGTTCAGGGGTAAAAAGAAGAGAATGTTTTCGGCGGAAATGTAGTTGGGAACGATGGAGCTGAAGGCAAACAGAAAAACCCCGAGGGGAACAAAATCAGCaac                                                                          >  2:474523/1‑249 (MQ=255)                                                                                                                                                                  agCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCa                                                                       <  1:306393/251‑1 (MQ=255)                                                                                                                                                                    cAGTAAAATGGCGGACAAATTGGTAATAGACATGCAGGCCA‑TTATGATATCTGCCAGTTGCAAAATCAGAGGAAGAATTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCCCAAACACCAGATCGCTTTAGGGATATTCAGGCGTAAAAAGAAAAGCATGTTTTCGGCGCGAATGTAGTTGGCAAAGAAGGAGCTGAAAGCAAACAAAATAACCACGAGGGTAACAAAcactgaacaccagg                                                                     >  2:514613/1‑238 (MQ=255)                                                                                                                                                                      gTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGaa                                                                   <  1:585904/251‑1 (MQ=255)                                                                                                                                                                           attgCGGTCAAAGTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTTCCGCCTATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCAt                                                              <  1:652653/248‑1 (MQ=255)                                                                                                                                                                              cccgatctatggggaaTAGCCATGCAGGCCA‑TTATGATATCTGCCAGTGGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTaa                                                           <  2:157409/237‑2 (MQ=255)                                                                                                                                                                                ggTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTaaca                                                         <  1:297440/251‑4 (MQ=255)                                                                                                                                                                                ggTCAAATTGGTAATAGCCAAGCAGGCCA‑TTAGGATATCTGCCAGTTGCCACAACAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTaaca                                                         <  2:60162/251‑4 (MQ=255)                                                                                                                                                                                    aaaTTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGACGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCg                                                     >  1:472681/1‑251 (MQ=255)                                                                                                                                                                                      aTTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTaa                                                           >  1:157409/1‑242 (MQ=255)                                                                                                                                                                                                gCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGCAGACTTGGCACGGTGCCGCCGGTGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAAAACCCGCATCGCCTTc                                           >  2:532575/1‑249 (MQ=255)                                                                                                                                                                                                 ccATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTgg                                        <  1:410701/251‑1 (MQ=255)                                                                                                                                                                                                 ccATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTgg                                        >  1:410889/1‑251 (MQ=255)                                                                                                                                                                                                 ccATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTgg                                        >  1:556621/1‑251 (MQ=255)                                                                                                                                                                                                 ccATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGACCCCCAGGAACCCATTAACACCCGCATCGCCTTCTgg                                        >  2:56726/1‑251 (MQ=255)                                                                                                                                                                                                  cATGCAGGCCA‑TTATGAGATCTGCCAGGTGCCACATCAGCGGAATACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGa                                       <  2:62207/251‑1 (MQ=255)                                                                                                                                                                                                           cA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGAt                                      <  2:230314/243‑1 (MQ=255)                                                                                                                                                                                                           cA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGAt                                      >  1:230314/1‑243 (MQ=255)                                                                                                                                                                                                            a‑TTATGATATCTGCCAGATGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGTAGGCAAACAGAATAGCCACGAGGGTAACCAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAAt                             >  2:560600/1‑251 (MQ=255)                                                                                                                                                                                                               tATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATAc                           >  1:183429/1‑251 (MQ=255)                                                                                                                                                                                                                                  cacaTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGt        >  1:38219/1‑251 (MQ=255)                                                                                                                                                                                                                                  cacaTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGAAGGt        >  1:328806/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTtgt     >  1:632483/1‑251 (MQ=255)                                                                                                                                                                                                                                        aGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGACAAACAGAATAACCACGCGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGcctggagtgtgtgcc  >  2:289806/1‑237 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    GTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAAT‑GGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCA‑TTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGA‑TCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCC  >  NC_000913/6511‑6990

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