Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,470,410	(T)7→6	coding (1280/1332 nt)	gtrS →	serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,470,404	0	T	.	100.0%	269.8 / NA	61	coding (1274/1332 nt)	gtrS	serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (39/22);  total (39/22)

TAATAAAATAAAGGCATCTTTATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCC  >  NC_000913/2470255‑2470544                                                                                                                                                      |                                                                                                                                             taataaAATAAAGGCATCTTTATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtt                                                                                                                                             >  1:750421/1‑151 (MQ=255)       aataaaGGCATCTTTATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                                                                       >  2:516030/1‑149 (MQ=255)          aaaGGCATCTTTATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAAttt                                                                                                                                   >  1:976944/1‑151 (MQ=255)                aTCTTTATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaaa                                                                                                                             >  2:293565/1‑151 (MQ=255)                  cTTTATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAaata                                                                                                                           >  1:482517/1‑151 (MQ=255)                     tATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAAtt                                                                                                                        >  1:790538/1‑151 (MQ=255)                     tATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAAtt                                                                                                                        >  2:77356/1‑151 (MQ=255)                        ttGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCaa                                                                                                                     >  1:943125/1‑151 (MQ=255)                          gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAg                                                                                                                   >  2:645251/1‑151 (MQ=255)                                 ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                            >  2:507968/1‑151 (MQ=255)                                       ttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                                                                      <  2:477250/151‑1 (MQ=255)                                                  cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                           >  2:1164229/1‑151 (MQ=255)                                                  cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                           >  1:1184704/1‑151 (MQ=255)                                                    aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                                                                         <  1:9388/151‑1 (MQ=255)                                                         aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                    >  1:953795/1‑151 (MQ=255)                                                            tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                                                                  <  1:267261/150‑1 (MQ=255)                                                                 tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                                                                            <  2:941475/151‑1 (MQ=255)                                                                 tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                                                                            >  2:1053220/1‑151 (MQ=255)                                                                 tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                                                                            <  2:1098956/151‑1 (MQ=255)                                                                 tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                                                                            <  2:482344/151‑1 (MQ=255)                                                                   tttCGGGGAAGGATATAGAGATTTAAGAATGCATTTATTTGGAATGTATTTTTCGTTCGACCTGTGCTTATACATAACAGTCG‑TTTTTTAATTGATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTGTAAATTg                                                                          <  2:881932/151‑1 (MQ=255)                                                                            aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                 >  2:466725/1‑151 (MQ=255)                                                                             ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                 >  1:653698/1‑150 (MQ=255)                                                                             ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                 >  2:498075/1‑150 (MQ=255)                                                                             ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                 >  1:1230809/1‑150 (MQ=255)                                                                               atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                                                              >  1:348959/1‑151 (MQ=255)                                                                               atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                                                              >  1:377540/1‑151 (MQ=255)                                                                                      gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg                                                        >  2:952471/1‑150 (MQ=255)                                                                                                 gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                            >  1:1237061/1‑151 (MQ=255)                                                                                                  caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                           <  1:676478/151‑1 (MQ=255)                                                                                                  caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                           <  1:298605/151‑1 (MQ=255)                                                                                                    tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                         <  1:779184/151‑1 (MQ=255)                                                                                                    tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                         <  2:750105/151‑1 (MQ=255)                                                                                                    tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                         <  2:656103/151‑1 (MQ=255)                                                                                                    tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                         <  1:133578/151‑1 (MQ=255)                                                                                                    tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                         <  1:294403/151‑1 (MQ=255)                                                                                                       atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCaa                                      >  1:203179/1‑151 (MQ=255)                                                                                                        tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                     <  2:207297/151‑1 (MQ=255)                                                                                                           ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTa                                  >  2:535131/1‑151 (MQ=255)                                                                                                                gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                              >  1:187106/1‑150 (MQ=255)                                                                                                                   tttttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                          <  2:932889/151‑1 (MQ=255)                                                                                                                              aCCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCGAATCATGGCGCGCACAAGCTATAATACCAACCTAATTACTCCTCCTCCTagag               <  2:942776/151‑1 (MQ=255)                                                                                                                              aCCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTagag               <  2:851923/151‑1 (MQ=255)                                                                                                                              aCCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATA‑CAACCTAATTTCTCCTCCTCTTagag               <  2:1265046/150‑1 (MQ=255)                                                                                                                                cTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTg             >  1:525997/1‑151 (MQ=255)                                                                                                                                   tGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa            >  1:1023037/1‑149 (MQ=255)                                                                                                                                   tGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa            >  2:289666/1‑149 (MQ=255)                                                                                                                                   tGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAca          >  1:1195027/1‑150 (MQ=255)                                                                                                                                        ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTagag               <  1:592005/141‑1 (MQ=255)                                                                                                                                        ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTagag               >  2:591753/1‑141 (MQ=255)                                                                                                                                        ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat      >  1:1096158/1‑150 (MQ=255)                                                                                                                                        ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat      >  1:337318/1‑150 (MQ=255)                                                                                                                                        ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat      >  1:991890/1‑150 (MQ=255)                                                                                                                                        ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat      >  1:1237517/1‑150 (MQ=255)                                                                                                                                        ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat      >  1:185060/1‑150 (MQ=255)                                                                                                                                          aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct    >  1:340206/1‑150 (MQ=255)                                                                                                                                          aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct    >  2:1159533/1‑150 (MQ=255)                                                                                                                                          aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct    >  2:359957/1‑150 (MQ=255)                                                                                                                                          aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct    >  1:178313/1‑150 (MQ=255)                                                                                                                                           cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAGCTAATTTCTCCTCCTCTTAGAGTGACTATAtctc   >  2:1192767/1‑150 (MQ=255)                                                                                                                                           cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc  <  1:712810/151‑1 (MQ=255)                                                                                                                                           cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc  <  1:1185993/151‑1 (MQ=255)                                                                                                                                             tAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtctc   >  2:904032/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                             TAATAAAATAAAGGCATCTTTATTTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCC  >  NC_000913/2470255‑2470544

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