Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,883,678	G→A	P18S (CCA→TCA)	dnaA ←	chromosomal replication initiator protein DnaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,883,678	0	G	A	100.0%	143.0 / NA	40	P18S (CCA→TCA)	dnaA	chromosomal replication initiator protein DnaA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/21);  total (19/21)

GTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTC  >  NC_000913/3883432‑3883923                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      gtgCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAAc                                                                                                                                                                                                                                                   >  1:342713/1‑251 (MQ=255)     gggggCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCAt                                                                                                                                                                                                                                               <  2:311310/250‑1 (MQ=255)            gcgACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTTGGGCCTCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCATCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCaa                                                                                                                                                                                                                                        <  2:197623/251‑1 (MQ=255)            gcgACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCaa                                                                                                                                                                                                                                        <  1:833/251‑1 (MQ=255)                                cGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGcc                                                                                                                                                                                                                    >  1:62430/1‑251 (MQ=255)                                    tgtgGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCaaa                                                                                                                                                                                                                 <  2:38557/250‑1 (MQ=255)                                         ggTTTGCGTCACCGGTTTGGGGCCGACGTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGaaa                                                                                                                                                                                                           <  2:307565/250‑1 (MQ=255)                                                    ccGGTTTGGTGCCGACTTCAACACGCAGCTGGGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAATGTACTTGTCCCGTACCCAATCGAGGTCAAAACGGTTTGGCGCGTACAGGGCCATCGTGTTATCGCTTAGTTCCGCCTTCAAGGTGCGTATCCACATACTGAATTCTGTTGCTGATAATTCATCCTGCAATCGGGCAAGACACTGCTGTCAAAGCTAAAGGGACAcggcg                                                                                                                                                                                                <  2:57694/251‑1 (MQ=255)                                                           ggTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTcc                                                                                                                                                                                          <  1:246023/250‑1 (MQ=255)                                                            gTGCCGACTTCTAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAACACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCAc                                                                                                                                                                                        >  2:229139/1‑251 (MQ=255)                                                                         aaagCAGATGTGGGGCAACCGGTCCGCAGAAACTGGTTAGGAGTCCATTGATAATATTAAGGTAACTGGCCCGTACCCAAACCAGGACAAAACGGTTTGGCGCGGAAAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACCGAAATCTGGGGCTTATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTc                                                                                                                                                                           <  1:227499/248‑1 (MQ=255)                                                                                       gCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAgg                                                                                                                                                             <  1:226924/251‑1 (MQ=255)                                                                                         aTCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCt                                                                                                                                                           >  2:195658/1‑251 (MQ=255)                                                                                                     gAAACTGGTTAGCAGTCCATTGATATTATTACGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATg                                                                                                                                               >  1:26151/1‑251 (MQ=255)                                                                                                         cTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttg                                                                                                                                           >  2:96659/1‑251 (MQ=255)                                                                                                          tGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttgt                                                                                                                                          >  2:138109/1‑251 (MQ=255)                                                                                                          tGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttgt                                                                                                                                          >  2:315152/1‑251 (MQ=255)                                                                                                           ggTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttgtt                                                                                                                                         >  1:101590/1‑251 (MQ=255)                                                                                                                     ccATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGacgtacgt                                                                                                                                <  1:1448/250‑1 (MQ=255)                                                                                                                          gATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGaaa                                                                                                                           >  2:111598/1‑250 (MQ=255)                                                                                                                          gATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGaaa                                                                                                                           >  2:296418/1‑250 (MQ=255)                                                                                                                                     ggTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACCCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTg                                                                                                               >  2:274733/1‑251 (MQ=255)                                                                                                                                           tGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGt                                                                                                          <  1:195658/250‑1 (MQ=255)                                                                                                                                            gTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAACGACCCTGCTTGAGGGTg                                                                                                         >  2:749/1‑250 (MQ=255)                                                                                                                                            gTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGAGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTg                                                                                                         >  1:332070/1‑250 (MQ=255)                                                                                                                                                            aGGACAAACCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCAGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCACTCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGCCTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGACAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTc                                                                                        >  2:309523/1‑251 (MQ=255)                                                                                                                                                                      ggTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTACGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTtat                                                                               >  1:358677/1‑250 (MQ=255)                                                                                                                                                                           ggCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGccc                                                                         <  1:114239/251‑1 (MQ=255)                                                                                                                                                                           ggCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGccc                                                                         <  1:138109/251‑1 (MQ=255)                                                                                                                                                                                   cAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCa                                                                    >  1:334774/1‑248 (MQ=255)                                                                                                                                                                                     gggCCAGCGTGTTATCTCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGGCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCACGATCAAGATc                                                               <  2:240809/251‑1 (MQ=255)                                                                                                                                                                                                  aTCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAt                                                   <  2:114305/250‑1 (MQ=255)                                                                                                                                                                                                                 cTGCAATGGGCGTATCCACATAATGAATTCTGTGGCTGATAACTCATACTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                   >  2:232200/1‑251 (MQ=255)                                                                                                                                                                                                                   gCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                   <  1:749/249‑1 (MQ=255)                                                                                                                                                                                                                    cAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGc                                <  1:164729/251‑1 (MQ=255)                                                                                                                                                                                                                    cAATGGGCGTATCAACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGCCACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACCTTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGCCCGCGGATCATAGc                                >  2:118292/1‑251 (MQ=255)                                                                                                                                                                                                                          gCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                          <  1:76014/251‑1 (MQ=255)                                                                                                                                                                                                                                          aatgcTGTGGCTGATAACTCATCCTGGAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc          <  2:2972/247‑1 (MQ=255)                                                                                                                                                                                                                                             tCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg        <  2:116607/250‑1 (MQ=255)                                                                                                                                                                                                                                                  ggCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc  <  1:328060/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTC  >  NC_000913/3883432‑3883923

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