Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	397,808	A→G	E390E (GAA→GAG)	sbmA →	peptide antibiotic/peptide nucleic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	397,808	0	A	G	100.0%	251.3 / NA	66	E390E (GAA→GAG)	sbmA	peptide antibiotic/peptide nucleic acid transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/32);  total (34/32)

TACGATTACGCTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAG  >  NC_000913/397673‑397936                                                                                                                                        |                                                                                                                                 tACGATTACGCTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATg                                                                                                                              >  1:869930/1‑140 (MQ=255)         cGCTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGAt                                                                                                                      <  1:847207/140‑1 (MQ=255)           cTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATgg                                                                                                                    >  1:570917/1‑140 (MQ=255)              ggCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGa                                                                                                                 >  2:768611/1‑140 (MQ=255)              ggCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGa                                                                                                                 >  1:845316/1‑140 (MQ=255)                 cTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACaa                                                                                                              <  1:36679/140‑1 (MQ=255)                       aCGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCa                                                                                                        <  1:768611/140‑1 (MQ=255)                       aCGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCa                                                                                                        <  1:819258/140‑1 (MQ=255)                        cGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAg                                                                                                       >  1:1072881/1‑140 (MQ=255)                           aGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGaa                                                                                                    <  2:570917/140‑1 (MQ=255)                            gATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAg                                                                                                   <  2:63140/140‑1 (MQ=255)                             aTTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGt                                                                                                  <  2:869930/140‑1 (MQ=255)                              ttACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTa                                                                                                 >  1:612907/1‑140 (MQ=255)                                       tttttGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATAcc                                                                                        >  2:745459/1‑140 (MQ=255)                                       tttttGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATAcc                                                                                        <  1:165740/140‑1 (MQ=255)                                       tttttGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATAcc                                                                                        >  2:609363/1‑140 (MQ=255)                                           tGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGa                                                                                    >  1:918139/1‑140 (MQ=255)                                             gTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGc                                                                                  <  1:745459/140‑1 (MQ=255)                                               cAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTa                                                                                >  2:587746/1‑140 (MQ=255)                                                 ggTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTaaa                                                                              <  2:612907/140‑1 (MQ=255)                                                  gTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTaaaa                                                                             <  1:587746/140‑1 (MQ=255)                                                    tCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAg                                                                           <  2:1072881/140‑1 (MQ=255)                                                     cgcgGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAgg                                                                          <  2:747474/140‑1 (MQ=255)                                                     cgcgGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAgg                                                                          <  2:526347/140‑1 (MQ=255)                                                       cgGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGgag                                                                        <  2:417601/140‑1 (MQ=255)                                                            cTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTa                                                                   >  2:3635/1‑140 (MQ=255)                                                                 cAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTc                                                              <  1:349755/140‑1 (MQ=255)                                                                  aGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTcg                                                             <  2:400892/140‑1 (MQ=255)                                                                      ccTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTcgcgcg                                                         >  1:862615/1‑140 (MQ=255)                                                                       cTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAgg                                                                          <  2:153259/122‑1 (MQ=255)                                                                       cTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAgg                                                                          >  1:153259/1‑122 (MQ=255)                                                                         gATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCa                                                      >  1:272387/1‑140 (MQ=255)                                                                         gATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCa                                                      >  1:187472/1‑140 (MQ=255)                                                                            tAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATc                                                   <  2:187472/140‑1 (MQ=255)                                                                                          cacTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTg                                     >  1:1040968/1‑140 (MQ=255)                                                                                              ggTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCAt                                 >  1:8076/1‑140 (MQ=255)                                                                                               gttgagttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCtt                                          <  2:837062/130‑1 (MQ=255)                                                                                               gttgagttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCtt                                          >  1:837062/1‑130 (MQ=255)                                                                                                ttgagttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCAttg                               >  1:57284/1‑140 (MQ=255)                                                                                                    gttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTgg                           <  2:918139/140‑1 (MQ=255)                                                                                                    gttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTgg                           >  2:457640/1‑140 (MQ=255)                                                                                                    gttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTgg                           <  1:457640/140‑1 (MQ=255)                                                                                                    gttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTgg                           >  2:701486/1‑140 (MQ=255)                                                                                                         tGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGtt                      >  1:231004/1‑140 (MQ=255)                                                                                                         tGTCTATCTACAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGtt                      >  1:469934/1‑140 (MQ=255)                                                                                                                   cAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAg            <  2:845316/140‑1 (MQ=255)                                                                                                                   cAAACGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAg            <  2:231004/140‑1 (MQ=255)                                                                                                                     aaCGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg          >  2:407873/1‑140 (MQ=255)                                                                                                                     aaCGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg          <  1:3635/140‑1 (MQ=255)                                                                                                                       cGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        >  2:981041/1‑140 (MQ=255)                                                                                                                       cGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        <  1:981041/140‑1 (MQ=255)                                                                                                                       cGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        <  2:272387/140‑1 (MQ=255)                                                                                                                       cGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        >  2:233326/1‑140 (MQ=255)                                                                                                                       cGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        >  2:216155/1‑140 (MQ=255)                                                                                                                       cGTCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        <  1:407873/140‑1 (MQ=255)                                                                                                                         tCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggggg      >  2:298966/1‑139 (MQ=255)                                                                                                                         tCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggggg      >  2:98737/1‑139 (MQ=255)                                                                                                                         tCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg          <  2:578495/136‑1 (MQ=255)                                                                                                                         tCTGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg          >  1:578495/1‑136 (MQ=255)                                                                                                                           tGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        <  2:588436/136‑1 (MQ=255)                                                                                                                           tGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg        >  1:588436/1‑136 (MQ=255)                                                                                                                           tGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGgtgt    <  1:233326/140‑1 (MQ=255)                                                                                                                           tGCGCAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGGTTTAGCGGGgtgt    >  1:282553/1‑140 (MQ=255)                                                                                                                             cgcAGCTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGGgtag  >  1:1132030/1‑140 (MQ=255)                                                                                                                                  cTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAg            <  2:554732/125‑1 (MQ=255)                                                                                                                                  cTTTGAGCATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAg            >  1:554732/1‑125 (MQ=255)                                                                                                                                        |                                                                                                                                 TACGATTACGCTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAG  >  NC_000913/397673‑397936

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