Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,499,393	T→G	intergenic (+192/‑55)	cysG → / → yhfL	siroheme synthase/DUF4223 domain‑containing lipoprotein YhfL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,499,393	0	T	G	100.0%	168.1 / NA	47	intergenic (+192/‑55)	cysG/yhfL	siroheme synthase/DUF4223 domain‑containing lipoprotein YhfL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/22);  total (25/22)

TGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTGCTTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGG  >  NC_000913/3499162‑3499633                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 tgGCCTGCGCGATACACTGAACTGGTTCTCCAACCATTAATTTACACCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCCGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTAACCGTTACAGCGCCTTACAAAATTAAGTCATTCGCTCACTGCCTACAATTAAATCTTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCatcgaccagat                                                                                                                                                                                                                                >  2:342029‑M1/1‑239 (MQ=255)        cgcgATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgcc                                                                                                                                                                                                                         >  1:35303‑M1/1‑232 (MQ=255)              aaaCTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcg                                                                                                                                                                                                                   <  2:452870‑M1/250‑25 (MQ=255)              aaaCTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcg                                                                                                                                                                                                                   <  1:468556‑M1/250‑25 (MQ=255)                     actgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggc                                                                                                                                                                                                           <  1:134024‑M1/251‑33 (MQ=255)                                aaCCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttaaatccgt                                                                                                                                                                                                >  2:192202‑M1/1‑208 (MQ=255)                                                          aaTTCTGGCCTTCTTAATGACTCTTTTTATCTATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgtttt                                                                                                                                                                      >  2:233216‑M1/1‑182 (MQ=255)                                                                  ccTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATCTTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaccgcgg                                                                                                                                                              >  2:331378‑M1/1‑174 (MQ=255)                                                                            gACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtca                                                                                                                                                     <  1:551009‑M1/250‑87 (MQ=255)                                                                                          atTCAGCAAATAAAACAAGGCCCCTTAATAAATATTATCTGCTAATTGGTTAGATTCTCCGTTCCAGCGACTTGAAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTAGCGTTCTTAACGGGCAGTCGTCGTATATTCctcgaccagataactgccattgcgcccacggcggtcaatccgtatcgctccgcgtgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaac                                                                                                                                      <  2:243967‑M1/251‑102 (MQ=255)                                                                                                            aGCCCCTTAATAAATATTATCTGCTAATGTGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaccgcggcctgcgtcaggccaacggatgaacggcattcatttgccagct                                                                                                                    <  2:425259‑M1/251‑120 (MQ=255)                                                                                                                     ataaataTTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatgg                                                                                                           >  2:169155‑M1/1‑123 (MQ=255)                                                                                                                           attatCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggagttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttc                                                                                                      >  2:139658‑M1/1‑117 (MQ=255)                                                                                                                                  gCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgatt                                                                                              >  1:477143‑M1/1‑110 (MQ=255)                                                                                                                                   cTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtctcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaccggatgaacggcattcatttgccagctcatctatggttttcaatctattc                                                                                             >  2:552364‑M1/1‑109 (MQ=255)                                                                                                                                     aaTGGTTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattac                                                                                            <  1:225013‑M1/250‑144 (MQ=255)                                                                                                                                        gggTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgca                                                                                        <  1:449261‑M1/251‑148 (MQ=255)                                                                                                                                         ggTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttaaatcgattactgca                                                                                        >  2:2293‑M1/1‑103 (MQ=255)                                                                                                                                                   aCCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcaccc                                                                              >  2:199765‑M1/1‑93 (MQ=255)                                                                                                                                                         aCAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgcctccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaccaaat                                                                       >  2:25652‑M1/1‑87 (MQ=255)                                                                                                                                                             cGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgctcttttcacccaacaaataacg                                                                   <  1:102098‑M1/251‑169 (MQ=255)                                                                                                                                                             cGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacg                                                                   >  2:238987‑M1/1‑83 (MQ=255)                                                                                                                                                                           aaGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcccccaacaaataacgggcca                                                              >  2:371708‑M1/1‑69 (MQ=255)                                                                                                                                                                           aaGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggcca                                                              <  1:371708‑M1/242‑174 (MQ=255)                                                                                                                                                                             gTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataaagggccagcgtcgtgatg                                                   >  2:555058‑M1/1‑67 (MQ=255)                                                                                                                                                                                   cGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacg                                             >  1:138110‑M1/1‑61 (MQ=255)                                                                                                                                                                                       cACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgcct                                          <  2:138110‑M1/250‑194 (MQ=255)                                                                                                                                                                                          tGACGCCAATTAAATATTACTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataaatgccattgcgcacacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgacagctcatctatggttttcaatcgattaccgcacttttcacacaacaaataacggtccagcgtcatgaggacgacgacttctt                                      >  2:492099‑M1/1‑54 (MQ=255)                                                                                                                                                                                           gACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttt                                     <  1:435442‑M1/251‑199 (MQ=255)                                                                                                                                                                                           gACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggacaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttt                                     >  1:346387‑M1/1‑53 (MQ=255)                                                                                                                                                                                               ccAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgacaagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccaacgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttg                                 >  2:402310‑M1/1‑49 (MQ=255)                                                                                                                                                                                                     aaaTATTCCTGTATCGTTCTTAACGAGCAGTATATGTATATGCctctaccagataactgccattgcgtccacggcgttcaatccggatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggtgttcaagcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgtcttcttttttgataaac                           <  2:367163‑M1/251‑209 (MQ=255)                                                                                                                                                                                                         atTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcgtttttgataaacgtt                        <  1:139658‑M1/250‑212 (MQ=255)                                                                                                                                                                                                          ttCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctctaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgagtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattacctcacttttcacccaacaaataacggaccagcgtcgtgatgacgacgccttcttttttgataaacgttcg                      >  1:6518‑M1/1‑38 (MQ=255)                                                                                                                                                                                                           tCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcg                      <  2:304407‑M1/250‑214 (MQ=255)                                                                                                                                                                                                            ccTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgac                    <  1:438171‑M1/251‑216 (MQ=255)                                                                                                                                                                                                            ccTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataacagacattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcatcgtgatgacgacgccttcttttttgataaacgttcgac                    >  1:576452‑M1/1‑36 (MQ=255)                                                                                                                                                                                                            ccTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCcccgaccagataactgccattgcgcccacggcggtcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgac                    <  1:102707‑M1/251‑216 (MQ=255)                                                                                                                                                                                                                 aTCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataactggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcata               <  1:494835‑M1/251‑221 (MQ=255)                                                                                                                                                                                                                 aTCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgac                                                 <  1:439503‑M1/217‑187 (MQ=255)                                                                                                                                                                                                                 aTCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgac                                                 >  2:439503‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                   cGTTCTTACCGAGCCGTCTTCGTATATTCctcgaccagataaccgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttctcccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataa              <  2:124706‑M1/250‑222 (MQ=255)                                                                                                                                                                                                                   cGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccccggcgttcaatccgtagcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataat             <  2:419986‑M1/251‑223 (MQ=255)                                                                                                                                                                                                                        ttAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttcataaacgttcgcagcataataaatt        >  2:540152‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                             gAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataataaattttccg   >  1:529341‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                              aGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataataaattttcagt  <  1:559017‑M1/251‑234 (MQ=255)                                                                                                                                                                                                                               gCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcatactaaattttcagt  >  2:424195‑M1/1‑17 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 TGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTGCTTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGG  >  NC_000913/3499162‑3499633

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