Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,637,436	G→A	intergenic (‑26/‑206)	yhiN ← / → pitA	putative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,637,436	0	G	A	98.8%	305.3 / NA	80	intergenic (‑26/‑206)	yhiN/pitA	putative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/56);  minor base C (0/1);  total (23/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATAC  >  NC_000913/3637189‑3637676                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 atgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagggcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTg                                                                                                                                                                                                                                              <  2:518821‑M2/29‑1 (MQ=255)   gtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGa                                                                                                                                                                                                                                             <  1:529261‑M2/30‑1 (MQ=255)   gtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGa                                                                                                                                                                                                                                             >  2:283460‑M2/221‑250 (MQ=255)     gcacatcaaccacgactgctgtctggaaagcgccgtgttgaaaggggacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaattccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAtt                                                                                                                                                                                                                                          <  2:447613‑M2/33‑1 (MQ=255)     gcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAtt                                                                                                                                                                                                                                          <  1:343363‑M2/33‑1 (MQ=255)       acatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGt                                                                                                                                                                                                                                        >  1:29106‑M2/217‑251 (MQ=255)           caaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAg                                                                                                                                                                                                                                    <  2:542187‑M2/39‑1 (MQ=255)               cacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgttcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCGCTTA‑AAAAAACGGACGTGAATTGTAAAGATTc                                                                                                                                                                                                                                <  1:512855‑M2/43‑1 (MQ=255)                 cgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCt                                                                                                                                                                                                                              <  1:386319‑M2/45‑1 (MQ=255)                  gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaattccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                             <  2:523249‑M2/46‑1 (MQ=255)                  gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                             >  1:384433‑M2/206‑251 (MQ=255)                      actgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTg                                                                                                                                                                                                                         <  1:32323‑M2/50‑1 (MQ=255)                      aatgtctggaaatcgccgtgttgaaaggtgacatgggagacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacacagggccactagcacacctgattcctctctaaatactctatccgaGTTGCTCTCCTCTTA‑AACAAACGGACGTGAATTGTAAAGTTTCCTCAGTg                                                                                                                                                                                                                         >  2:630636‑M2/202‑251 (MQ=255)                                aatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAtt                                                                                                                                                                                                                                          >  2:613838‑M2/192‑224 (MQ=255)                                aatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAtt                                                                                                                                                                                                                                          <  1:613838‑M2/33‑1 (MQ=255)                                   cgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCAt                                                                                                                                                                                                            <  1:230402‑M2/63‑1 (MQ=255)                                    gccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                           >  1:169703‑M2/188‑251 (MQ=255)                                    gccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgagtcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                           <  2:201775‑M2/64‑1 (MQ=255)                                     ccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATc                                                                                                                                                                                                           >  1:566127‑M2/187‑250 (MQ=255)                                     ccgtgttgaaaggtgacatgggtgacgtgcagcatgttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCa                                                                                                                                                                                                          <  1:292543‑M2/65‑1 (MQ=255)                                       gtgttgaaagttgtcatgggtgacgtgcagcattttgccgatgacgtgatcgcgcagctcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctcgctaaataccctatccgaGTTGCGCTCCTCGTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGACCCAGCATCAat                                                                                                                                                                                                        <  2:321397‑M2/67‑1 (MQ=255)                                       gtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAat                                                                                                                                                                                                        <  1:22628‑M2/67‑1 (MQ=255)                                       gtgttgaaaggtgacatgggtcacgtgcagcattttgccgatgacgttatcgccccgcgcggcgtgcggcaccgtcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctgtggtgaacgatcaacaaggggccactagcacgcctgattcctctctaattaccctagccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCATTGGTCGCGCCAGCATCAat                                                                                                                                                                                                        <  2:541373‑M2/67‑1 (MQ=255)                                         gttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAatat                                                                                                                                                                                                      <  1:448970‑M2/69‑1 (MQ=255)                                         gttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAata                                                                                                                                                                                                       <  2:614787‑M2/68‑1 (MQ=255)                                              aaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgctgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                 <  1:262043‑M2/74‑1 (MQ=255)                                              aaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcagcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTa                                                                                                                                                                                                 <  1:268268‑M2/74‑1 (MQ=255)                                               aggtgagatgggtgacgtgcagcattttgccgattacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcccttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGGTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaa                                                                                                                                                                                                <  2:105930‑M2/75‑1 (MQ=255)                                               aggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattccgctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaa                                                                                                                                                                                                <  2:189952‑M2/75‑1 (MQ=255)                                                   gacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAg                                                                                                                                                                                             <  1:248180‑M2/78‑1 (MQ=255)                                                    acatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGa                                                                                                                                                                                           <  1:567610‑M2/80‑1 (MQ=255)                                                      atgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAg                                                                                                                                                                                         <  1:150547‑M2/82‑1 (MQ=255)                                                       tgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGt                                                                                                                                                                                        >  1:566747‑M2/169‑251 (MQ=255)                                                        gggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGt                                                                                                                                                                                        <  1:466497‑M2/83‑1 (MQ=255)                                                         ggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcattagcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAt                                                                                                                                                                                      <  1:225558‑M2/85‑1 (MQ=255)                                                                   agcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAAt                                                                                                                                                                            <  1:50629‑M2/95‑1 (MQ=255)                                                                   agcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAAt                                                                                                                                                                            <  2:398491‑M2/95‑1 (MQ=255)                                                                     cattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATta                                                                                                                                                                          <  1:559410‑M2/97‑1 (MQ=255)                                                                      attttgccgatgacgttatcgcccagcgcggagtgcggcacgggcatttgcagtgcgtgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaatacccgatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtat                                                                                                                                                                         <  2:528194‑M2/98‑1 (MQ=255)                                                                         ttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtatat                                                                                                                                                                       <  2:240949‑M2/100‑1 (MQ=255)                                                                          tgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtatata                                                                                                                                                                      <  1:22590‑M2/101‑1 (MQ=255)                                                                           gccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAg                                                                                                                                                                    <  1:463787‑M2/103‑1 (MQ=255)                                                                           gccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAg                                                                                                                                                                    <  2:73055‑M2/103‑1 (MQ=255)                                                                             cgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGAt                                                                                                                                                                  <  1:138044‑M2/105‑1 (MQ=255)                                                                             cgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGAt                                                                                                                                                                  <  1:51872‑M2/105‑1 (MQ=255)                                                                                              agcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTg                                                                                                                                                  >  1:208651‑M2/130‑250 (MQ=255)                                                                                               gcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGAGAATTATGTTTTGAGTGaa                                                                                                                                                <  2:236436‑M2/123‑1 (MQ=255)                                                                                                         cggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAg                                                                                                                                      <  1:283460‑M2/133‑1 (MQ=255)                                                                                                         cggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGACAATTATTTTTTGAGTGAAATCCATAc                                                                                                                                        >  1:568812‑M2/119‑249 (MQ=255)                                                                                                             acgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAggggg                                                                                                                                  >  1:551378‑M2/115‑251 (MQ=255)                                                                                                              cgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGc                                                                                                                                 <  2:474299‑M2/138‑1 (MQ=255)                                                                                                               gggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCa                                                                                                                                <  1:565527‑M2/139‑1 (MQ=255)                                                                                                                  cagttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAATATAATTATTTTTTGAGTGAAAGCCATACAGGGGGCAAAt                                                                                                                             <  2:535925‑M2/142‑1 (MQ=255)                                                                                                                             cttgacgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAATAaagac                                                                                                                  >  2:618085‑M2/99‑249 (MQ=255)                                                                                                                                     aggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACTGGGGGCAAATCAAAAAAAGTCTATATTTc                                                                                                          >  2:347874‑M2/91‑251 (MQ=255)                                                                                                                                     aggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTc                                                                                                          <  1:296257‑M2/161‑1 (MQ=255)                                                                                                                                        aagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACt                                                                                                       <  2:384433‑M2/164‑1 (MQ=255)                                                                                                                                           attgagtctattcttgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTg                                                                                                    >  2:448481‑M2/85‑251 (MQ=255)                                                                                                                                           attgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTg                                                                                                                                                                               >  1:534654‑M2/85‑176 (MQ=255)                                                                                                                                           attgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAGTTGCCTg                                                                                                                                                                               <  2:534654‑M2/92‑1 (MQ=255)                                                                                                                                              gagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCGCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAAGCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGccc                                                                                                 <  2:169703‑M2/170‑1 (MQ=255)                                                                                                                                                  ctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgc                                                                                             >  1:210178‑M2/78‑251 (MQ=255)                                                                                                                                                      tattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAGTGTAAAGATTCCTCAGTGGTCGCTCCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGGGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgcgccgc                                                                                          <  1:477813‑M2/177‑1 (MQ=255)                                                                                                                                                            gtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTc                                                                                   >  2:520868‑M2/68‑251 (MQ=255)                                                                                                                                                                      gaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAAt                                                                          <  1:520868‑M2/193‑1 (MQ=255)                                                                                                                                                                       aacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATgc                                                                        <  1:120253‑M2/195‑1 (MQ=255)                                                                                                                                                                          gatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGcgcc                                                                     >  1:494115‑M2/54‑251 (MQ=255)                                                                                                                                                                           atcaacaaagggccactagctcacctgattccactctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTGTTGAGTTAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAACGcgccg                                                                    <  1:3909‑M2/199‑1 (MQ=255)                                                                                                                                                                                        cactagcacacctaactcctctctaaataccctatccgaGTTGCACTCCTCTTA‑AAAACACGGCCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGGAAATCAAAAAAAGTCTATATTTCACTGTGCCCGCGCCGCGAATGTCACTGCTGATGCGCTGCGTTCATGTCCTc                                                       <  2:340268‑M2/212‑1 (MQ=255)                                                                                                                                                                                                 acctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCg                                              <  2:29106‑M2/221‑1 (MQ=255)                                                                                                                                                                                                                   taccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAt                            >  2:97252‑M2/13‑251 (MQ=255)                                                                                                                                                                                                                            cgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCt                   <  1:448481‑M2/248‑1 (MQ=255)                                                                                                                                                                                                                             gaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAt                            <  2:528001‑M2/239‑1 (MQ=255)                                                                                                                                                                                                                             gaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAt                            >  1:528001‑M2/3‑241 (MQ=255)                                                                                                                                                                                                                              aGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTgg                 >  1:481418‑M2/2‑251 (MQ=255)                                                                                                                                                                                                                                       cctcTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAttt        <  1:123685/251‑1 (MQ=255)                                                                                                                                                                                                                                          tcttaaaaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGc      >  1:424268/6‑251 (MQ=255)                                                                                                                                                                                                                                             taaaaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCAt    <  1:354420/248‑1 (MQ=255)                                                                                                                                                                                                                                              taaaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTCTATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGTATTTGCATAc  >  2:511104/2‑251 (MQ=255)                                                                                                                                                                                                                                              aaaaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATAc  <  1:383483/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 GAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATAC  >  NC_000913/3637189‑3637676

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