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%	161.7 / NA	43	S17L (TCA→TTA)	yhdY	putative ABC transporter membrane subunit YhdY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/20);  total (23/20)

ACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAA  >  NC_000913/3421128‑3421612                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               aggttaatgacatctatgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaagccacccgtggtaaccagacccgtgctgcgctgatgatgggcatgaaccgtggtacgttgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaaggcttgattaaaaaggcgctactcggcattgggaagcgtcttttttataggtgtcacaaaggaCATTACTTAAGc                                                                                                                                                                                                                                            <  2:42607‑M2/12‑1 (MQ=255)     gagtgacctatatgagctggtactggctgaaggagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattaaggttagctaaatgcttgtttaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCcgcg                                                                                                                                                                                                                                        >  2:351771‑M2/236‑251 (MQ=255)              ctatgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagaaccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccgtttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGgt                                                                                                                                                                                                                                 <  2:146404‑M2/23‑1 (MQ=255)                                tgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGt                                                                                                                                                                                                             >  1:229619‑M2/209‑251 (MQ=255)                                      agaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTcc                                                                                                                                                                                                        <  2:411659‑M2/48‑1 (MQ=255)                                          cagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTg                                                                                                                                                                                                   >  1:340167‑M2/199‑251 (MQ=255)                                                         atggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCctgct                                                                                                                                                                                    >  1:255746‑M2/184‑251 (MQ=255)                                                                      acacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATATCCTGCTGACTATTGGCTg                                                                                                                                                                        <  2:171227‑M2/80‑1 (MQ=255)                                                                          ccgtggtaaccagacacgtgatgcgctgatgatgggcatcaaccgtggtacgatgtgtaaaaaattgaaaaaatacggcatgaactaattcaggttagcaaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTTTGGGTGCGAAAAAATCTGGTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTgcatct                                                                                                                                                                   >  2:158217‑M2/167‑249 (MQ=255)                                                                           cgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATg                                                                                                                                                                  >  2:76904‑M2/166‑251 (MQ=255)                                                                                 aaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGa                                                                                                                                                             <  1:338054‑M2/91‑1 (MQ=255)                                                                                         ccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAAc                                                                                                                                                    >  1:410954‑M2/152‑251 (MQ=255)                                                                                             gctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGa                                                                                                                                                 <  2:421105‑M2/103‑1 (MQ=255)                                                                                                gcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATcc                                                                                                                                              <  1:375736‑M2/106‑1 (MQ=255)                                                                                                         atgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCt                                                                                                                                    >  1:221914‑M2/136‑251 (MQ=255)                                                                                                             gcatcaaccgtgatacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttaactaaatgcttgattaaaaaggcgctactcggcatggggaaacgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTGGATGTGGGAACTGATCCCACCGTTGCTGAAc                                                                                                                                <  1:223588‑M2/120‑1 (MQ=255)                                                                                                               atcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaatgccggctagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCTCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAGCTg                                                                                                                              <  1:197847‑M2/122‑1 (MQ=255)                                                                                                                         gtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTg                                                                                                                    >  1:255736‑M2/120‑251 (MQ=255)                                                                                                                          tacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgattgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGACTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGc                                                                                                                   >  2:248403‑M2/119‑251 (MQ=255)                                                                                                                              ctgcgtacaaaattgaaaaaatacggcatggactatttcaggttagctaaatgattgattaaaaaggcgctactcggcatggggaagcgccttttttatcggtgtcacaaaggaAATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTATCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCCTATGGTTGATGTGGGAACTGGTCCCACCGTTGCTGAACTGGGCATTTTTGCaggg                                                                                                               >  2:143582‑M2/115‑250 (MQ=255)                                                                                                                                  gtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAAc                                                                                                           >  2:383203‑M2/111‑251 (MQ=255)                                                                                                                                   taaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACt                                                                                                          >  1:186466‑M2/110‑251 (MQ=255)                                                                                                                                      aaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTggg                                                                                                       >  1:93108‑M2/107‑251 (MQ=255)                                                                                                                                            gaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCt                                                                                                 >  1:440576‑M2/101‑251 (MQ=255)                                                                                                                                            gaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCt                                                                                                 >  1:325795‑M2/101‑251 (MQ=255)                                                                                                                                                    acggcatgaaataattcaggttagctaaatgcttgattaaaaaggcgctactcggaatggggaagagccttttttataggtgtcacacaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGCTCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGt                                                                                         >  2:140351‑M2/93‑251 (MQ=255)                                                                                                                                                         atgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCg                                                                                     >  2:191737‑M2/88‑250 (MQ=255)                                                                                                                                                           gaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACt                                                                                  <  2:282796‑M2/166‑1 (MQ=255)                                                                                                                                                            aactaagtcaggttagctaaatgcttgattaaaaaggcgctactcggcatggtgaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAAGCTAACTGGGTTGGCTCAACGCGTGCCGACTg                                                                                 <  2:45717‑M2/167‑1 (MQ=255)                                                                                                                                                                       gttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttagaggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCg                                                                      <  2:245625‑M2/178‑1 (MQ=255)                                                                                                                                                                         tagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGt                                                                    >  1:401070‑M2/72‑251 (MQ=255)                                                                                                                                                                         tacctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGAGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCATGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCGGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCACCGCGTGCCGAATGCACAAAAGCCGGt                                                                    >  2:285330‑M2/72‑251 (MQ=255)                                                                                                                                                                          agctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGt                                                                    <  2:287433‑M2/180‑1 (MQ=255)                                                                                                                                                                               aatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGt                                                              >  1:381732‑M2/66‑251 (MQ=255)                                                                                                                                                                                                                  aagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGAGCGATGGTGTGGGTGCGAAAAAAGCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtat                           <  2:272582‑M2/221‑1 (MQ=255)                                                                                                                                                                                                                    gcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgt                         <  2:325795‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                     cgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgta                        <  1:158063‑M2/224‑1 (MQ=255)                                                                                                                                                                                                                       ccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgtatg                      <  1:159019‑M2/226‑1 (MQ=255)                                                                                                                                                                                                                       ccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgtatg                      <  2:27702‑M2/226‑1 (MQ=255)                                                                                                                                                                                                                           ttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCt                  >  2:408942‑M2/22‑251 (MQ=255)                                                                                                                                                                                                                                taggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAAAGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTAcc             >  1:213441‑M2/17‑251 (MQ=255)                                                                                                                                                                                                                                        acaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGAc     <  1:76904‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                           aaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCaa  <  2:93108‑M2/246‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               ACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAA  >  NC_000913/3421128‑3421612

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