Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	289,256	G→C	6.5%	intergenic (‑94/+45)	yagI ← / ← argF	CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	289,256	0	G	C	6.5%	180.4 / 3.5	59	intergenic (‑94/+45)	yagI/argF	CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTC  >  NC_000913/289109‑289399                                                                                                                                                    |                                                                                                                                                gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtg                                                                                                                                                 <  2:185879/148‑1 (MQ=255)  ttGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgtg                                                                                                                                               >  6:154023/1‑149 (MQ=255)  ttGAACAGGTCGAGGATCTGCAACGCACGTGCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGAACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtggg                                                                                                                                               <  4:248460/149‑3 (MQ=255)    gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTgg                                                                                                                                              >  6:23880/1‑148 (MQ=255)      aCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                                                                                                                                           <  4:143797/149‑1 (MQ=255)      aCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                                                                                                                                           <  8:330323/149‑1 (MQ=255)        aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGagag                                                                                                                                         <  8:46297/149‑1 (MQ=255)         ggTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGACCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                        <  6:113477/149‑1 (MQ=255)          gTCGTGGATCTGGAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                       <  4:128956/149‑1 (MQ=255)          gTCGAGGATCTTCAACGCACGTTCAACAGACTGAATAATAGGCATAATGCTGGCATGTACACGATGGAATTAAGGCGATAACGTACCTTTTTCGGCATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGAGGTAGAggg                                                                                                                                       >  6:89819/1‑149 (MQ=255)            cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAgggag                                                                                                                                     <  6:48964/149‑3 (MQ=255)             gAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                    >  8:55798/1‑149 (MQ=255)             gAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                    <  8:205237/149‑1 (MQ=255)               ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                  <  6:26890/149‑1 (MQ=255)                 aTCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTAGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                >  5:261003/1‑149 (MQ=255)                    tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                   <  8:270186/143‑1 (MQ=255)                    tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTCTGGGAGAGGGCCgg                                                                                                                                   >  7:270186/1‑143 (MQ=255)                      cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                               <  5:23880/145‑1 (MQ=255)                      cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                            <  6:120271/148‑1 (MQ=255)                          gCACGTTCAACAGACTTAATAATCGCCATACTGCGGGAATGTCCACGCTGTAATTACGGCGAAAACTTACTTTTTTCGGTATTAAAAATCAATGAATTGGAGCCGGGGTTCTCCCTCTCCGTGGGGGAGAGGGCCGGGGTGAGGGCAcc                                                                                                                       <  8:319118/149‑1 (MQ=255)                            aCGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAc                                                                                                                        >  6:295590/1‑146 (MQ=255)                            aCGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAc                                                                                                                        <  5:295590/146‑1 (MQ=255)                            aCGTTCAACAGACTGAATAATCGCCATAATGCTGGCCTGCCCACGCTGTAATTACGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTATCCCTGTGGGAGAGGGCCGGGGTGTGGGCACCaa                                                                                                                     <  5:2559/149‑1 (MQ=255)                             cGTTAAACAGACTTAATAAGCGGCATACTGCGGGCTTGGCCACGCGGTAATTACGGCCAAAACGTCCCTTTTTCGGGAGGAAAAATCAATTAAATGGAGCCGGGGTTCTCCCTCTCCCGGGGGGAGTGGGCCGGGGTGAGGGCACCAAc                                                                                                                    <  2:84823/149‑1 (MQ=255)                              gTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcg                                                                                                                   >  8:282456/1‑149 (MQ=255)                               ttCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgc                                                                                                                  >  6:333908/1‑149 (MQ=255)                                 cAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGAGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgcg                                                                                                                 >  1:29897/1‑148 (MQ=255)                                                tCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCaa                                                                                                                     >  1:48181/1‑129 (MQ=255)                                                tCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCaa                                                                                                                     <  2:48181/129‑1 (MQ=255)                                                  ggCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGCTGAGGGCAACAACGCGCAGCGGACCCAATCACTc                                                                                               >  4:329410/1‑149 (MQ=255)                                                            tGGCTTGTCCACGCTGTAATTAAGGCCAAACCTTACTTTTTTCGGGATTAAAAATCAATTAATTGGCGCCGGTTTTCCCCCTCTCCGTGGGGTAGTGGGCGGGGGTGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGtt                                                                                     <  4:152688/149‑1 (MQ=255)                                                                   gccACGCTGGAATTACGGCGAAAACGTACCTTTTTCGCGAGGAAAAATCAATGAAATGGCGCCGGTGTTCTCCCTCTCCCTGGGGGAGTGGGCGGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatc                                                                               <  2:29897/147‑1 (MQ=255)                                                                   cccACGCTGGAATTACGGCTAAAACTTACCTTTTGCGGTATTAAAAATCAATCAAAGGGACCCGGTGTGCGCCCTCTCCCTGTGGGAGAGGGCCGGGGGGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatc                                                                               <  8:181328/147‑1 (MQ=255)                                                                     cACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                        <  6:81071/89‑1 (MQ=255)                                                                     cACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                        >  5:81071/1‑89 (MQ=255)                                                                      aCGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCAtcatca                                                                           >  8:12109/1‑149 (MQ=255)                                                                        gCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGATCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc                                                                         >  3:127849/1‑149 (MQ=255)                                                                          tGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGc                                                                       >  7:54490/1‑149 (MQ=255)                                                                          tGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGc                                                                       >  4:76780/1‑149 (MQ=255)                                                                           ggAATTACGGCCAAAACGTACCTTTTTCGGGATCAAAAATCAATGAAGTGGAGCGGGTGTTCCCCCTCTCCCGGTGGGTGTGGGCCGGGTGGAGGGCACCATCGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc                                                                      <  6:137471/149‑1 (MQ=255)                                                                              aTTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCACTCACTCCCCAAGCGTTGCCATCATCACCGCCtt                                                                    >  5:334663/1‑148 (MQ=255)                                                                               taaaGGCGAAACCGTACCTTTTTCGTGATTAAAATTAAATCAACGGCGGCCGGTGTTCGCCCTCACCCTGGGGGAGAGGGCCGGGGTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTaa                                                                  <  6:82549/147‑1 (MQ=255)                                                                                     ctaacaCGTACATTTTGCGGTATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCGCTCCCTGGGGGATTGGGCGGGGGTTAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                            <  1:176977/144‑1 (MQ=255)                                                                                     ctaaaaCGTAGTTTTTTCGGGTTGAAAAATAATTGAATTGGAGCCGGGGTTGGCCCGCTCCGTGTGGGTGAGGGCCGGTGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                            <  6:261003/147‑1 (MQ=255)                                                                                     caaaacCGTACATTTTGGGGTTTTAAAAATCAATCAATGGGCGCGGTTGTTCTCCCTAGCCGTGTGGCATGGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                            <  7:282456/147‑1 (MQ=255)                                                                                     cGAAAACTTACCTTTTTGGTTATGAAAAATCATTGAAAGGGACCCGGGTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGTGGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                            <  3:329410/149‑1 (MQ=255)                                                                                     cGAAAACGTACCTTTTTGGGGTTTAAAAATCAATGAATTGGAGCCGGTGTTCCCCTTCTCCCGTTGGTTTTGGGCGGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                            <  6:190063/149‑1 (MQ=255)                                                                                            gTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATctg                                                     >  7:145022/1‑147 (MQ=255)                                                                                                cTTTTTCGGGATGAAAAATCAATGACATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGttt                                                 <  5:154023/149‑1 (MQ=255)                                                                                                 tttttCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCGGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGtttt                                                >  5:23174/1‑149 (MQ=255)                                                                                                              aaaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCCTATGCATCCGGTTTTCCGCCTGGTCGaa                                   >  1:264574/1‑149 (MQ=255)                                                                                                              aaaaTCAATGAAATGGAGCCGGTGTTCTCCATCTCCCTGTGGGAGAGGGCCCGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACACCCCAAGCGTTCCCATCTTCACCGCCTTAATCGTAAGCATCCGGTTTTCCCCCTGGTCCaa                                   >  8:254228/1‑149 (MQ=38)                                                                                                                    aaTGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                             >  1:126490/1‑149 (MQ=255)                                                                                                                                      ttCTCCCTC‑CCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCata           <  3:108270/148‑3 (MQ=255)                                                                                                                                       cggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCaaa           >  1:117386/4‑148 (MQ=255)                                                                                                                                          ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACAcc       >  4:195160/1‑149 (MQ=255)                                                                                                                                           cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt      >  8:157657/1‑149 (MQ=255)                                                                                                                                            ctctcgcgGTGGGTGCGGGCCGGGGTCAGGCCCCCAACGCGCAGCGGACCCAATCCCTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc     <  2:264574/149‑1 (MQ=255)                                                                                                                                            ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc     >  2:228234/1‑149 (MQ=255)                                                                                                                                             tctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGCCTCAAACACcgcg    >  5:226689/1‑149 (MQ=255)                                                                                                                                               tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtc  >  7:330699/1‑149 (MQ=255)                                                                                                                                               tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtc  >  5:87183/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                GTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTC  >  NC_000913/289109‑289399

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