Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	621,409	A→G	F261L (TTT→CTT)	fepD ←	ferrienterobactin ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	621,409	0	A	G	100.0%	256.1 / NA	68	F261L (TTT→CTT)	fepD	ferrienterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (43/25);  total (43/25)

ACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACT  >  NC_000913/621261‑621550                                                                                                                                                     |                                                                                                                                              aCGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAggg                                                                                                                                             >  1:81792/1‑151 (MQ=255)   gcgcAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAggg                                                                                                                                             >  1:158810/1‑149 (MQ=255)    cgcAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCa                                                                                                                                           >  2:1162891/1‑150 (MQ=255)     gcAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAAt                                                                                                                                         >  2:748516/1‑151 (MQ=255)       aGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCg                                                                                                                                       >  1:1098787/1‑151 (MQ=255)              cGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAAc                                                                                                                                >  2:1125221/1‑151 (MQ=255)               gggCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACt                                                                                                                               <  2:220102/151‑1 (MQ=255)                  cacaATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTAtt                                                                                                                            <  1:435287/151‑1 (MQ=255)                         aCCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCg                                                                                                                     >  2:1249533/1‑151 (MQ=255)                            cgcccgATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCAc                                                                                                                  >  2:414631/1‑151 (MQ=255)                               ccgATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTAc                                                                                                               >  1:414125/1‑151 (MQ=255)                                cgATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTAcc                                                                                                              >  1:241371/1‑151 (MQ=255)                                 gatgacATCGGCAAACAGCACCAGGGCAGGGGTAGCAAGTAGCGTGACGGGAAGCGACCACCGAAGATCCGCACCCACCACCCAACCCGCCATATGCGGCATCATCAGGCCCATACGGGCAATCGGGCCAACAATTGCTGTCGCACTAcc                                                                                                              >  1:707964/1‑150 (MQ=255)                                    gatATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACcaca                                                                                                           >  2:237899/1‑150 (MQ=255)                                      tatCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACaa                                                                                                          >  1:502851/1‑149 (MQ=255)                                       atCGGCCAACAGCAGCGGGGCAGGGGTAGCAAGTAGCGTGACGGGCATCGACCAGCGATGATCGGCACCTGCCAGCTTACGCGCCATATGCGGCACCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGTCGTCGCTCTACCACaaa                                                                                                         <  1:1254547/149‑1 (MQ=255)                                       atCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAg                                                                                                        <  1:701426/150‑1 (MQ=255)                                       atCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGc                                                                                                       <  1:455266/151‑1 (MQ=255)                                        tCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGc                                                                                                       >  2:206247/1‑150 (MQ=255)                                        tCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGc                                                                                                       >  2:1107001/1‑150 (MQ=255)                                         cGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCAc                                                                                                     >  1:729582/1‑151 (MQ=255)                                            cAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCATgg                                                                                                   >  2:252642/1‑150 (MQ=255)                                            cAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACgg                                                                                                   >  1:288497/1‑150 (MQ=255)                                            cAAACAGCAGCAGGGCAGCGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGt                                                                                                  >  1:980381/1‑151 (MQ=255)                                             aaaCAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGt                                                                                                  >  2:628129/1‑150 (MQ=255)                                             aaaCAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTa                                                                                                 <  2:87909/151‑1 (MQ=255)                                              aaCAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTa                                                                                                 >  2:885903/1‑150 (MQ=255)                                              aaCAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTa                                                                                                 >  2:596509/1‑150 (MQ=255)                                                cagcagcagGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAAt                                                                                               >  2:641049/1‑150 (MQ=255)                                                 agcagcagGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCg                                                                                             >  1:1217047/1‑151 (MQ=255)                                                           cagggGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGAcc                                                                                   >  2:1163435/1‑151 (MQ=255)                                                                  aGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAAc                                                                            >  2:1078982/1‑151 (MQ=255)                                                                    cAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACtg                                                                          >  1:179394/1‑151 (MQ=255)                                                                         aGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTgc                                                                     >  2:325527/1‑151 (MQ=255)                                                                         aGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTgc                                                                     <  1:400917/151‑1 (MQ=255)                                                                           cGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTgcgc                                                                   <  1:242697/151‑1 (MQ=255)                                                                           cGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTgcgc                                                                   <  1:99678/151‑1 (MQ=255)                                                                            gTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTgcgcg                                                                  <  1:572141/151‑1 (MQ=255)                                                                                cGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCAc                                                              >  2:786833/1‑151 (MQ=255)                                                                                    cAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCAcgcg                                                           >  2:909122/1‑150 (MQ=255)                                                                                     aGCGACCAGCGATGATCGGCACCCAACAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGTCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGATAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGa                                                          >  1:1250743/1‑150 (MQ=255)                                                                                       cGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTg                                                       >  2:1217900/1‑151 (MQ=255)                                                                                         aCCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGcc                                                     >  2:1139205/1‑151 (MQ=255)                                                                                          ccAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGccc                                                    >  1:330207/1‑151 (MQ=255)                                                                                            aGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGccc                                                    >  2:254546/1‑149 (MQ=255)                                                                                            aGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGccc                                                    >  2:705176/1‑149 (MQ=255)                                                                                              cGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGATAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcg                                                <  2:137424/151‑1 (MQ=255)                                                                                               gatgatCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAgcgc                                               <  1:139738/151‑1 (MQ=255)                                                                                                   atCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTc                                           >  2:920662/1‑151 (MQ=255)                                                                                                   atCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTc                                           >  1:208513/1‑151 (MQ=255)                                                                                                   atCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTc                                           <  1:414601/151‑1 (MQ=255)                                                                                                        cacccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCgg                                       >  1:1143712/1‑150 (MQ=255)                                                                                                              ccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGc                                 >  2:1207313/1‑150 (MQ=255)                                                                                                              ccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCt                                <  1:1225227/151‑1 (MQ=255)                                                                                                               caGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTg                               <  2:502879/151‑1 (MQ=255)                                                                                                               caGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCACCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTg                               >  1:1148098/1‑151 (MQ=255)                                                                                                                aGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGc                              >  1:351705/1‑151 (MQ=255)                                                                                                                    aaCGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAg                          <  1:1249495/151‑1 (MQ=255)                                                                                                                     aCGCGCCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAgg                         <  2:4842/151‑1 (MQ=255)                                                                                                                       gcgcCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCt                       <  2:1098827/151‑1 (MQ=255)                                                                                                                        cgcCATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTc                      <  1:131539/151‑1 (MQ=255)                                                                                                                           cATATGCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCaaa                   <  1:895182/151‑1 (MQ=255)                                                                                                                                gCGGCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGt               <  1:860267/150‑1 (MQ=255)                                                                                                                                  tgCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc             >  1:714190/2‑150 (MQ=255)                                                                                                                                  ggCATCATCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCa            <  2:651860/151‑1 (MQ=255)                                                                                                                                    catcatCAGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgc          <  1:350684/151‑1 (MQ=255)                                                                                                                                       catcaGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcg       >  1:450440/1‑151 (MQ=255)                                                                                                                                        atcaGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcgc      >  2:354238/1‑151 (MQ=255)                                                                                                                                        atcaGGCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcgc      <  2:891686/151‑1 (MQ=255)                                                                                                                                            ggCCAATAAGGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACt  <  1:1268377/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              ACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACT  >  NC_000913/621261‑621550

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