Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,624,205	G→T	P195P (CCC→CCA)	ytjB ←	protein Smp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,624,205	0	G	T	100.0%	197.7 / NA	59	P195P (CCC→CCA)	ytjB	protein Smp
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/37);  total (22/37)

AACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCG  >  NC_000913/4624071‑4624341                                                                                                                                       |                                                                                                                                         aaCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGa                                                                                                                                     >  2:1096318/1‑140 (MQ=255)       gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATtgctg                                                                                                                               >  1:530367/1‑140 (MQ=255)       gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATtgctg                                                                                                                               >  2:377273/1‑140 (MQ=255)              aaGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCgg                                                                                                                        <  1:710012/140‑1 (MQ=255)               agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCggg                                                                                                                       <  1:377273/140‑1 (MQ=255)               agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCggg                                                                                                                       <  1:778060/140‑1 (MQ=255)               agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCggg                                                                                                                       <  2:1170578/140‑1 (MQ=255)               agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCggg                                                                                                                       <  2:567482/140‑1 (MQ=255)               agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCggg                                                                                                                       <  2:589026/140‑1 (MQ=255)                  gTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGc                                                                                                                    <  1:626220/140‑1 (MQ=255)                   tCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCg                                                                                                                   >  1:589128/1‑140 (MQ=255)                     agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGt                                                                                                                      >  1:268582/1‑135 (MQ=255)                     agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGt                                                                                                                      <  2:268582/135‑1 (MQ=255)                     agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGtt                                                                                                                 <  1:209462/140‑1 (MQ=255)                     agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGtt                                                                                                                 <  2:638088/140‑1 (MQ=255)                     agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGtt                                                                                                                 >  2:828507/1‑140 (MQ=255)                     agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGtt                                                                                                                 >  1:286421/1‑140 (MQ=255)                     agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGtt                                                                                                                 <  2:938093/140‑1 (MQ=255)                       agaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCTCTCTTTTGTCTCGCTTTCTTGCTCTGCCGGCACCGGTTTGCTGGCGGCTAACAGGAATGGCGATTGCTGCCAGCGGGTGcttgta                                                                                                               <  2:1140518/140‑6 (MQ=255)                       agaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTa                                                                                                               <  2:907502/140‑1 (MQ=255)                       agaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTa                                                                                                               <  2:589128/140‑1 (MQ=255)                           gAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCt                                                                                                           >  1:721827/1‑140 (MQ=255)                            aGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTg                                                                                                          >  1:38931/1‑140 (MQ=255)                                 gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTAcc                                                                                                             >  2:1198033/1‑132 (MQ=255)                                 gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTAcc                                                                                                             <  1:1198033/132‑1 (MQ=255)                                 gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTagca                                                                                                     <  1:23409/140‑1 (MQ=255)                                 gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTagca                                                                                                     >  1:1056538/1‑140 (MQ=255)                                      aTGTGGACATAACTATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTTCCCTGTATCAGCGTg                                                                                                <  1:1122478/140‑1 (MQ=255)                                      aTGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTg                                                                                                <  1:109216/140‑1 (MQ=255)                                         tgGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCgg                                                                                             <  2:401731/140‑1 (MQ=255)                                           gACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGt                                                                                           <  2:721827/140‑1 (MQ=255)                                           gACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGt                                                                                           <  1:962049/140‑1 (MQ=255)                                           gACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGt                                                                                           >  1:63084/1‑140 (MQ=255)                                           gACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGt                                                                                           <  2:1171849/140‑1 (MQ=255)                                           gACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGCGCGGGt                                                                                           >  1:586630/1‑140 (MQ=255)                                           gACATAACGATTTCCTTTCTTGTAGTAATGGGCCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGCTTACAGGAATGGCGCTTGCCGACAGCCGGTGCGGTTACACCCTAACAGCGTGCGGGt                                                                                           >  2:90210/1‑140 (MQ=255)                                                aaCGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCa                                                                                      >  2:327539/1‑140 (MQ=255)                                                aaCGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCa                                                                                      >  2:164468/1‑140 (MQ=255)                                                aaCGATTCCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCa                                                                                      >  1:925082/1‑140 (MQ=255)                                                 aCGATTTCCTTTCTTGTAGTAATGGTTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCgggg                                                                                           >  2:917298/1‑133 (MQ=255)                                                 aCGATTTCCTTTCTTGTAGTAATGGTTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGt                                                                                           <  1:917298/134‑1 (MQ=255)                                                     tttCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTaca                                                                                 >  2:838999/1‑140 (MQ=255)                                                       tCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACAcc                                                                               <  1:861369/140‑1 (MQ=255)                                                       tCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACAcc                                                                               <  1:838999/140‑1 (MQ=255)                                                                    taATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTgag                                                                  <  1:327539/140‑1 (MQ=255)                                                                    taATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTgag                                                                  <  2:63084/140‑1 (MQ=255)                                                                           tCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGGCAGCGGGGGCGTTTACCCTGTTGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAgcagca                                                           >  1:332057/1‑140 (MQ=255)                                                                           tCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGCTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCCATTGGCCATGAGAgcagca                                                           >  2:496069/1‑140 (MQ=255)                                                                                 ttttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCa                                                     <  1:828507/140‑1 (MQ=255)                                                                                 ttttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCa                                                     <  2:651334/140‑1 (MQ=255)                                                                                  tttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTTCCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAg                                                    <  1:496069/140‑1 (MQ=255)                                                                                   ttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAgg                                                   >  1:564347/1‑140 (MQ=255)                                                                                   ttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAgg                                                   >  2:495538/1‑140 (MQ=255)                                                                                         cGCTTTCTTCCTCTTCCCGCACCGGTTTGCTGGTGGTTAACATGAATGCCGATTGCCGTCCGGCTGTTCGCCTTACCTGTAGCAGCGTGCGGGGTGGCACGACTCCTATTTCCAGTGAGGACCACGACATCAGGCGGaaa                                             >  2:1119878/1‑140 (MQ=255)                                                                                          gCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCGAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTaaaa                                            <  1:90210/140‑1 (MQ=255)                                                                                                   ctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAgtgg                                   <  2:564347/140‑1 (MQ=255)                                                                                                      ttccGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGt                                         >  1:14129/1‑131 (MQ=255)                                                                                                      ttccGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGt                                         <  2:14129/131‑1 (MQ=255)                                                                                                        ccGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTa                              >  1:193407/1‑140 (MQ=255)                                                                                                        ccGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTa                              <  2:530367/140‑1 (MQ=255)                                                                                                          ggCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATc                            >  2:25799/1‑140 (MQ=255)                                                                                                              ccGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCAcc                        <  2:925082/140‑1 (MQ=255)                                                                                                                 gTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGt                     <  1:1096318/140‑1 (MQ=255)                                                                                                                     gCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTggg                 <  1:164468/140‑1 (MQ=255)                                                                                                                            gTTAACAGGAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACACCGATTGCCCGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCggg          >  2:81133/1‑139 (MQ=255)                                                                                                                                    gAATGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCg  <  1:503108/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         AACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCG  >  NC_000913/4624071‑4624341

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