Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,133,311	G→A	P232P (CCG→CCA)	metF →	5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,133,311	0	G	A	93.8%	146.0 / ‑0.6	48	P232P (CCG→CCA)	metF	5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (15/30);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATCACACAGTGATCACAAAT  >  NC_000913/4133070‑4133546                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             aTCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATg                                                                                                                                                                                                                                     <  1:274151/251‑1 (MQ=255)             aCCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTc                                                                                                                                                                                                                         <  1:680383/251‑1 (MQ=255)                            gcgcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGgatgat                                                                                                                                                                                                          <  2:688369/251‑1 (MQ=255)                                  aGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGcc                                                                                                                                                                                                    <  1:28799/251‑1 (MQ=255)                                  aGGCGGATTTGCTTAATCTGAAACGAAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATGGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGTCCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGcc                                                                                                                                                                                                    <  1:614390/251‑1 (MQ=255)                                        aTTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAAcc                                                                                                                                                                                              <  1:40360/251‑1 (MQ=255)                                          ttGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCTACGGTCTGGATGATGATGCCGAAAccc                                                                                                                                                                                             <  1:220502/250‑1 (MQ=255)                                          ttGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAAccc                                                                                                                                                                                             <  1:597134/250‑1 (MQ=255)                                                aaTCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACt                                                                                                                                                                                      <  2:268917/251‑1 (MQ=255)                                                          gCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGcgcg                                                                                                                                                                            <  2:597197/251‑1 (MQ=255)                                                                tGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCTATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTTTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGTATTTTGCCGGTATCTAACTTTAAACAGGCGAAGGAAT‑TTGCCGATATCACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAatat                                                                                                                                                                       <  1:361935/250‑1 (MQ=255)                                                                 ggATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCGTTGTTGTGGG‑AATTTTT‑TCTGGGGTTTTGTCGGTATTTTACTTTTAAAAAGAGAAGAAAT‑TTTCTGGTTTTACCAACGTGCGTATTTCAGCGTGGGTGGCGCAAATGTTTTACGGTGTGGATGGTGTTGCCGAAAcccccgcaccggttggggcggcgattt                                                                                                                                                                     >  2:342631/1‑227 (MQ=255)                                                                    tGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGGCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAatat                                                                                                                                                                       <  1:368058/246‑1 (MQ=255)                                                                        ggAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGa                                                                                                                                                              <  1:375125/251‑1 (MQ=255)                                                                                      gATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGAttt                                                                                                                                                <  1:342171/251‑1 (MQ=255)                                                                                             cAGTTCCTCTTCGATGCTGAAAGCTACCTGCGTTTTCGTGACCGCTGTGGATCGGCGGGCATTGATGGGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTGGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCg                                                                                                                                        <  2:454338/251‑1 (MQ=255)                                                                                                              cGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTTGTGTGGG‑AATTATT‑TCGGGGATTTTGCCGGGATCTAAATTTTAAAAGGCGAAGAAAT‑TTTCCGGGATGACCAACGTGCGTATTTCAGCGGGGATGGGGCAAAAGTTCGACGGGCTGGGTGATTATGCCGCAAACCACAAACTGGGTGGGGCGGATATTGTCATGGGTATGGTGAAGATTTTTAAGC‑GAGAAGGAGGGAAAGAtt                                                                                                                       >  1:424958/1‑251 (MQ=255)                                                                                                              cGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAgttttgaacgatt                                                                                                                       <  1:325519/251‑13 (MQ=255)                                                                                                                  aGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCa                                                                                                                   <  1:710466/251‑1 (MQ=255)                                                                                                                    cTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTGATGTGGA‑AATTATTTCCGGGA‑TTTTTCCGGTGTCTTACTTTTAACAGGCGGAGAAAT‑TTTCCGGCATGACCCACGTGCGTATTCCAGCGGGGATGGGGCCAATGTTCGACGGTCTGGGTTATGATGCCGAAACCCGCCAACCGGTTGGGGCGGATATTGCCATGGGGTTGGTGAAGGTTTT‑AAGCCGTGAAGGAGGGAAAGAttttccct                                                                                                                 >  1:404643/1‑246 (MQ=255)                                                                                                                           cGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTGATGTGGGGATTTTTT‑CGGGGATTTTTCCCGGT‑TCTTACTTTTAAAAGGCGAAGAAATATTTCCGA‑ATGACCAAAGGGAGGATTTCGGCGTGGGTGGGGGAAATGTTCGAAGGTCTGGGTGATGATGCCGACACCCCCACAAAGGTTGGGGCGCAGATTTCCATGGATATGGTGAAGATTTT‑TAACCGTGGAGGGGAGAAAAATTTTCAATTCtata                                                                                                           >  2:246408/1‑250 (MQ=255)                                                                                                                                tCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCtt                                                                                                      <  1:41869/250‑1 (MQ=255)                                                                                                                                    gACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATg                                                                                                                                                                                                                            <  1:89850/129‑1 (MQ=255)                                                                                                                                             gtATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGAGTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGaaa                                                                                         <  1:479093/250‑1 (MQ=255)                                                                                                                                                 cggcggGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGt                                                                                    <  1:605076/251‑1 (MQ=255)                                                                                                                                                 cggcggGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGt                                                                                    >  1:267706/1‑251 (MQ=255)                                                                                                                                                    cggGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAc                                                                                 >  1:345219/1‑251 (MQ=255)                                                                                                                                                                       ttatt‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTTCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGGTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGTCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTc                                                                                                                     <  1:117824/197‑1 (MQ=255)                                                                                                                                                                       ttatt‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGCCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATCTGGTGCAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTc                                                                                                                     >  2:117824/1‑197 (MQ=255)                                                                                                                                                                         att‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATAAGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTggg                                                             >  2:282366/1‑251 (MQ=255)                                                                                                                                                                              cGGGAATTTTGCCGGTATCTAACTTTAAACAGTCGAAGAAAT‑TTGACGCTATGACCAACGTGCGTATTCCAGCGTGGATGACGAAAATGTTCGACGGTCTGGATGATGATGCCGAATCCCGCAAACTGGTTGGCGCGAATCTTGCAATGGATATGGTGAAACTTTT‑AAGCCGAGAAGGAGTGAAAGATTTCCACTTCTATCCGCTTAACCGGGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGtt                                                         >  2:486262/1‑251 (MQ=255)                                                                                                                                                                                ggAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCg                                                       <  1:589799/251‑1 (MQ=255)                                                                                                                                                                                   aTTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGAcc                                                    <  2:639614/251‑1 (MQ=255)                                                                                                                                                                                   aTTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGAcc                                                    >  1:5207/1‑251 (MQ=255)                                                                                                                                                                                     tttGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCt                                                   <  1:616567/250‑1 (MQ=255)                                                                                                                                                                                      ttGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TGGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGACAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGAGATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTg                                                  <  2:438374/250‑1 (MQ=255)                                                                                                                                                                                                        aaaCAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGCTTTCCACTTCTATACGCTTAACCGTGCTGACATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAACTAGTGTGGCt                               >  2:547152/1‑251 (MQ=255)                                                                                                                                                                                                         aaCAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCtt                              <  1:100565/251‑1 (MQ=255)                                                                                                                                                                                                            aGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTtg                           <  1:750361/251‑1 (MQ=255)                                                                                                                                                                                                            aGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCAGCGAGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAGACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGCAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATACATAGTGTGGCTTTtg                           >  1:386855/1‑251 (MQ=255)                                                                                                                                                                                                                gaagaaAT‑TTGCCGATATGACCAACGTGCGTATTCAAGCGTGGATGGAGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGCGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGa                        >  2:288242/1‑250 (MQ=255)                                                                                                                                                                                                                           gCCGATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCCAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑ACGCCGTGAAGGCGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATaacac               >  2:433962/1‑249 (MQ=255)                                                                                                                                                                                                                              gATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATCACACAGTg           >  1:648456/1‑250 (MQ=255)                                                                                                                                                                                                                              gATATGACCAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATCACACAGTGa          <  1:288242/251‑1 (MQ=255)                                                                                                                                                                                                                                tatGACCAACGTGCGCATTCCAGCGTGGATGTCGCAAATGTTCGACGGTCTGGATGATGATGCCGAATCCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATGTT‑AAGCCGTGAAGGAGTGAAAGATTTCCGCTTCTATACGCTTAACCGTGCGGAAATGAGTTACGCGATTTGCCATACGCTGTGTGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAAGCACACATTGAt         <  1:81465/250‑1 (MQ=255)                                                                                                                                                                                                                                      cAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATCACACAGTGAt         <  2:382515/244‑1 (MQ=255)                                                                                                                                                                                                                                      cAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATCACACAGTGAt         >  1:382515/1‑244 (MQ=255)                                                                                                                                                                                                                                      cAACGTGCGTATTCCAGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATCACACAGTGATCACAAAt  >  2:81770/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             ATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGA‑AATTATT‑CCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAT‑TTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTT‑AAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGACCTGGTTTATAAATAGTGTGGCTTTTGTGAAAATCACACAGTGATCACAAAT  >  NC_000913/4133070‑4133546

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