Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,421,374	C→T	S17L (TCA→TTA)	yhdY →	putative ABC transporter membrane subunit YhdY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,421,374	0	C	T	100.0%	239.1 / NA	68	S17L (TCA→TTA)	yhdY	putative ABC transporter membrane subunit YhdY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/38);  total (30/38)

TCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCG  >  NC_000913/3421130‑3421620                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       gtgaatgacctctatgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGcc                                                                                                                                                                                                                                                   <  1:569382‑M2/13‑1 (MQ=255)             taggagctggtactggctgaagtagaacagnccctgttggacatggtgaggcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaatgcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgcctttgttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGgtgt                                                                                                                                                                                                                                       <  2:362225‑M2/25‑1 (MQ=255)              atgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGtgtg                                                                                                                                                                                                                                      <  1:194302‑M2/26‑1 (MQ=255)              atgagcgggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTgg                                                                                                                                                                                                                                     <  2:701341‑M2/27‑1 (MQ=255)                gcgctggtcctggccgaagtagaaccgcccctgttggacatggtgatgcaatacacccgtgataaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgctggattaaaaaggcgctcctcggcatggggaagcgccttttttataggtggcacaaaggaCATAACTTAAGCCGCGTGATGGTGTggg                                                                                                                                                                                                                                    <  2:441971‑M2/28‑1 (MQ=255)                  gctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGc                                                                                                                                                                                                                                 >  1:424121‑M2/221‑251 (MQ=255)                      gtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaaa                                                                                                                                                                                                                             <  2:522129‑M2/35‑1 (MQ=255)                      gtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaaa                                                                                                                                                                                                                             <  1:332475‑M2/35‑1 (MQ=255)                       tactggctgaagtagaacagcccctgttggacatggttatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcctcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcctgaactaattcaggttagctaaatgcttgcttaaaaaggcgctcctcggcatggggaagcgccttttttataggtgtcacaaaggaCCTAACTTCAGCCACGCGATGGGGTGGGTGCGaaa                                                                                                                                                                                                                             >  2:552671‑M2/216‑250 (MQ=255)                          tggctgaagtagaacagcccctgttggacatggtgatgcaatacacccggggtaaccagacccgtgctgcgctgatgatgggcataaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcctgaaataattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAAt                                                                                                                                                                                                                         >  1:408977‑M2/213‑251 (MQ=255)                            gctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCt                                                                                                                                                                                                                       <  1:515077‑M2/41‑1 (MQ=255)                                        cagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTg                                                                                                                                                                                                           >  1:155582‑M2/199‑251 (MQ=255)                                         agcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTgg                                                                                                                                                                                                          <  2:67429‑M2/54‑1 (MQ=255)                                         agcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGGGTGGGTGCGAAAAAATCTGTTCTCCAGCTgg                                                                                                                                                                                                          >  1:440594‑M2/198‑251 (MQ=255)                                            ccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGc                                                                                                                                                                                                       >  1:43688‑M2/195‑251 (MQ=255)                                               tgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAAt                                                                                                                                                                                                    >  1:250457‑M2/192‑251 (MQ=255)                                                 ttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgacgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaaaa                                                                                                                                                                                                                            >  2:532297‑M2/190‑225 (MQ=255)                                                 ttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgacgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaaaa                                                                                                                                                                                                                            <  1:532297‑M2/36‑1 (MQ=255)                                                         ggtgatgcaatacacccgtggtaaccagaccngtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGa                                                                                                                                                                                          <  2:334247‑M2/70‑1 (MQ=255)                                                         ggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCctgctg                                                                                                                                                                                           >  1:331952‑M2/182‑250 (MQ=255)                                                               gcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTg                                                                                                                                                                                    <  1:700618‑M2/76‑1 (MQ=255)                                                                      gcccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaatgcagggtagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtggcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTGTTGGCTGCat                                                                                                                                                                             <  2:626088‑M2/83‑1 (MQ=255)                                                                            ggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaacaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttatttaaaggtgtcacaaaagaCATAACTTAAGCCGCGCGATGGTGTGGGGGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGtt                                                                                                                                                                       >  1:698126‑M2/163‑251 (MQ=255)                                                                                accagaaccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgcaaaaattcggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactaggcatggggaaacgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGAGATGGTGTGGGTGCGAAAAAATCTGTTCTCAAGCTGGAGCAAGAGCCTGCTGCCTATTGGCTGCATATGGTTGAtg                                                                                                                                                                   >  2:551084‑M2/159‑251 (MQ=255)                                                                                   agacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTgg                                                                                                                                                                >  1:598668‑M2/156‑251 (MQ=255)                                                                                   agacccgtgctgcgctgatgatgggcatcaacagtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGGGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGCGCAATAGCCTGCTGACTATTGGCTGCATCTGGTTGATGTgg                                                                                                                                                                >  2:556018‑M2/156‑251 (MQ=255)                                                                                     acccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCCTATGGTTGATGtgggg                                                                                                                                                              >  2:257998‑M2/154‑250 (MQ=255)                                                                                     acccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggagtcacaaaggaCATAACTTAAACCGCGCGATGGTGTGGGTGCGAAAAAATCTGATCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGGTGTggg                                                                                                                                                               >  1:358878‑M2/154‑250 (MQ=255)                                                                                       ccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAAc                                                                                                                                                            >  1:149072‑M2/152‑251 (MQ=255)                                                                                        cgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACGTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACt                                                                                                                                                           <  1:559943‑M2/101‑1 (MQ=255)                                                                                          tgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTg                                                                                                                                                          >  2:437380‑M2/149‑250 (MQ=255)                                                                                                   gatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCAcc                                                                                                                                                  <  1:523386‑M2/110‑1 (MQ=255)                                                                                                    atgatgggcatcaaccgtggtacgctgcgtaaaaaatggaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGGTGATGTGGGAACTGATCCCACCGtt                                                                                                                                               <  1:414587‑M2/113‑1 (MQ=255)                                                                                                       atgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCt                                                                                                                                            <  1:683594‑M2/116‑1 (MQ=255)                                                                                                          ggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGaa                                                                                                                                         <  1:227191‑M2/119‑1 (MQ=255)                                                                                                           gcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGGGGGAACTGATCCCACCGTTGCTGAAc                                                                                                                                        <  1:131434‑M2/120‑1 (MQ=255)                                                                                                             atcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggagtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTg                                                                                                                                      >  2:579110‑M2/130‑251 (MQ=255)                                                                                                               caaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgatgaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTTACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTggg                                                                                                                                    <  1:112647‑M2/124‑1 (MQ=255)                                                                                                                         acgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGTAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGCTCCCACCGTTGCTGAACTGGGCATTTTTGc                                                                                                                           >  1:28697‑M2/118‑250 (MQ=255)                                                                                                                          cgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAg                                                                                                                         <  1:683939‑M2/135‑1 (MQ=255)                                                                                                                           gctgagtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGTCTATTGTCTGCATATGGGTGATGTGGGAACTGATCCCACCTTTGCTGAACTGGGCATTTTTGCAgg                                                                                                                        <  1:607958‑M2/136‑1 (MQ=255)                                                                                                                                      aattgaaaaaatacggcatgaactaattcagnttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGtt                                                                                                             <  2:333721‑M2/147‑1 (MQ=255)                                                                                                                                         tgaaaaaatacggcattaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGATCTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTgg                                                                                                           <  1:572812‑M2/149‑1 (MQ=255)                                                                                                                                                       atgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCg                                                                                             >  1:303081‑M2/88‑250 (MQ=255)                                                                                                                                                         gaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACt                                                                                          >  1:349036‑M2/86‑251 (MQ=255)                                                                                                                                                          aactaattcaggttagctaaatgcttgattaaacaggcgctactcggcatggggacgcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTg                                                                                         >  2:526‑M2/85‑251 (MQ=255)                                                                                                                                                                   aggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGc                                                                                >  1:59877‑M2/76‑251 (MQ=255)                                                                                                                                                                            aaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCt                                                                        <  1:287525‑M2/184‑1 (MQ=255)                                                                                                                                                                                    gattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTc                                                                <  1:33221‑M2/192‑1 (MQ=255)                                                                                                                                                                                     attaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCtt                                                              <  2:53938‑M2/194‑1 (MQ=255)                                                                                                                                                                                      ttaaaaaggcgctactcggcatggggaagcgcattttttataggtgtcacaaagaaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGAAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGATGAAATGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGAGCAGACAGCTCAAAAGCCAGTGCGTGTTGGGActtt                                                             >  2:93693‑M2/57‑250 (MQ=255)                                                                                                                                                                                             ggcgctacccggcatggggaagcgtcttttttataggtgtcccaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCGGTTCTCCAGCGGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTTACTGGGTTGGCTCAACGCATGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATTCacg                                                      <  2:398082‑M2/202‑1 (MQ=255)                                                                                                                                                                                                    cgcggcatgggtaagcgccttttgtataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGAGCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGAtt                                               <  2:235696‑M2/209‑1 (MQ=255)                                                                                                                                                                                                      cggcatggggaagcgccttttttagaggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAAGCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGAt                                                <  1:313983‑M2/208‑1 (MQ=255)                                                                                                                                                                                                        gcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACATAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGCTTtggg                                           >  2:75133‑M2/39‑250 (MQ=255)                                                                                                                                                                                                         catggggaagcgccgtgtttataggtgtcacaaaggaCATAACTTAATCCGCGCGATGGTGTGGGTGCGAAAAAAGCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTc                                          <  2:149072‑M2/214‑1 (MQ=255)                                                                                                                                                                                                          atggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGCACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCAAAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGtcc                                         >  2:288335‑M2/37‑250 (MQ=255)                                                                                                                                                                                                                  gcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATAGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCAGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACGGCACAAAAGCTTGTGCCTGTTGGGTCTTCATCCACGAACGATGTGGTCAGTTtatgt                                 <  2:541349‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                        tttgtataggtgtcacaaagtaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTTTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATAGGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATggg                            <  2:303081‑M2/228‑1 (MQ=255)                                                                                                                                                                                                                            tataggggtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGTTTTGGTCAGTTTATGTATGGGCttt                        >  2:99267‑M2/19‑250 (MQ=255)                                                                                                                                                                                                                                 gtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCcac                  <  2:299251‑M2/238‑1 (MQ=255)                                                                                                                                                                                                                                   gtcacaaaggaCATAACTTAAGCCGCGCGTTGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCACTAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGACCGCTNTGGTCAGTTTATGTATGGGCTTTACCCacac                >  2:337201‑M2/12‑251 (MQ=255)                                                                                                                                                                                                                                      acaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGAc             <  1:257998‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                       caaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCG‑AAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCa           <  1:422374‑M2/244‑1 (MQ=255)                                                                                                                                                                                                                                          aggaCATAACTTAAGCCGCGCGATGGTGGGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAAc         <  2:155582‑M2/247‑1 (MQ=255)                                                                                                                                                                                                                                             aCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTTCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCt      <  1:288335‑M2/250‑1 (MQ=255)                                                                                                                                                                                                                                               aTAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTgg    >  1:91450/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aaCTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAACAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCg  >  2:58284/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       TCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCG  >  NC_000913/3421130‑3421620

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