Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,784,535	G→A	R520Q (CGG→CAG)	fadK →	short chain acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,784,535	0	G	A	97.4%	122.4 / ‑4.0	39	R520Q (CGG→CAG)	fadK	short chain acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (19/18);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GATGAGGCTGGCTATATAAAAATTACCGGACGCAAAAAAGATATTATTGTCCGCGGCGGCGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACGGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGG  >  NC_000913/1784294‑1784765                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       gATGAGGCTGGCTATATAAAAATTACCGGACGCAAAAAAGATATTATTGTCCGCGGCGGCGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCaaa                                                                                                                                                                                                                               <  1:357513/251‑1 (MQ=255)                           cGGACGCAAAAAAGATATTATTGTCCGCGGCGACGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGGTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTgg                                                                                                                                                                                                     <  1:404597/251‑1 (MQ=255)                            tgACGCAACAAAGATATTGTTGTCCGCGGCGGCGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGCTGCCTGAGTGATTGCAATATCCGATGAACGTTTAGGGGAACGATCATGCGCTTATGTCGGGCTGAAAGCGCCGCATCATTCATTATCGCTGGATGAGGCAGTGGCTTTTTTTAGCCGTAAACAGGCCGCAAAATATACATCTCCTGAACATAGCGTGGt                                                                                                                                                                                                    >  2:433683/2‑251 (MQ=255)                                       aGATATTATTGTCCGCGGCGGCGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAAc                                                                                                                                                                                         <  1:485388/251‑1 (MQ=255)                                               ttGTCCGCGGCGGCGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTCGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGAAAGAGGTAGTGGCGTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACcgcg                                                                                                                                                                                  >  1:616511/1‑250 (MQ=255)                                                  tCCGCGGCGGCGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTAACGCGAACt                                                                                                                                                                              >  1:229037/1‑251 (MQ=255)                                                     gcggcggcgAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTAcc                                                                                                                                                                           <  2:246880/251‑1 (MQ=255)                                                        gcggcgAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTc                                                                                                                                                                         >  1:340851/1‑250 (MQ=255)                                                             gAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGCCGCAAAATATAAATATCCTGAACATATCGTGGAAATCGAAAAACTACCGCGAACTACCTCAGGTa                                                                                                                                                                    >  2:577181/1‑250 (MQ=255)                                                                aaTATTAGCAGCCGTGAAATGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTAGTTGCAATGTCCTATGACCGTTTAGGTGAACGATCATGCGCTTATGTCGTGTTGAAAGCGCCGCATCACGCATTATCGCTGGAAGAGGTAGTGGCTTATTTTAGGCGTAAACAGGACGAAAAATATAAATATCCGGAACATATCGTGGTAATAGTAAAACTACCATGAACTACCGCAGTTAAAAt                                                                                                                                                                <  2:398258/251‑1 (MQ=255)                                                                     tANCAGCCTTGAAGTGGATGATAGTTTATGTCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATTTCCGATGAACGTTTATGTGTACGATCATGCACTTATGTCTTGCTGATAGCGCCGTATCATTCATTATCGCTTAAAGAGGTGGTGGCGTTTTTTAGCCGTAAACATGTCGCAAAATATAAAGATCCTGATCATATCGTGGTAATCGAAAAACTACCACTACCTACCTCAGGGAAAATTCaa                                                                                                                                                            <  2:609862/250‑1 (MQ=255)                                                                      agcagcCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACaaaa                                                                                                                                                          <  1:429897/251‑1 (MQ=255)                                                                                      aaGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATAAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGATCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGaaa                                                                                                                                           >  2:448624/1‑250 (MQ=255)                                                                                       aGATATTTTATTGCAGCATCATAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATGCAAAAGTTTTGGGTAAGAAAAg                                                                                                                                         <  2:639740/251‑1 (MQ=255)                                                                                                aTTGCAGCATCCTAAAATTCACGATGGCTGTGTGGTTACAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTACTGCAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGAAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGCACTACCTCAGGTCAAATACAAAAGTTTTTGTTAAGAAAAGATATTATgc                                                                                                                                >  2:623988/1‑251 (MQ=255)                                                                                                          ccTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATTTCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAAc                                                                                                                      >  1:96668/1‑251 (MQ=255)                                                                                                                      gATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCGTGAGATTATGTCGTGCTGAAAGCGCAGCATAATGCATTATCGCTGGAAGAGGTAGTGGGGTTTTTTAGCCGTAAACAGGTAGCAAAATATAAATATCCTGAACATATCGTGGGAATCGAAAAACGACCGCGAACTACATCAGGAAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCtg                                                                                                          >  1:406196/1‑251 (MQ=255)                                                                                                                           cTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGAAGTGCTGCAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAAATACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTATGTGa                                                                                                       >  1:91217/1‑249 (MQ=255)                                                                                                                                                 gAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTc                                                                               >  2:192524/1‑251 (MQ=255)                                                                                                                                                  aaCGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTTCCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAAATTCATCTTCg                                                                              <  2:386510/251‑1 (MQ=255)                                                                                                                                                   aCGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCgg                                                                             <  2:340851/251‑1 (MQ=255)                                                                                                                                                    cGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCggg                                                                            >  2:396721/1‑251 (MQ=255)                                                                                                                                                        tAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATc                                                                        >  1:44811/1‑251 (MQ=255)                                                                                                                                                            tGAACGATGATGCGCTTATGTCGTGCTGAAAGAGCCGCATCATTCATTATCGCTGGAAGAGGTGGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAACAACTACCGCGCACTACCTCAGGGAACATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACAt                                                                    >  2:188425/1‑251 (MQ=255)                                                                                                                                                             gAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACAt                                                                    <  2:184626/250‑1 (MQ=255)                                                                                                                                                                  aTCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGACGAGGTAGTGGCTTTTTTTAGCCGTAAACGGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAAcccc                                                              >  1:51351/1‑251 (MQ=255)                                                                                                                                                                    cATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAACCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCgg                                                            >  2:145002/1‑251 (MQ=255)                                                                                                                                                                            tATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCTGCGACTaa                                                     <  2:49001/250‑1 (MQ=255)                                                                                                                                                                                     ctaaaaGAGCCGCCTCATTCATTATCGCTGGAAGCGGTAATGGCTTTTTTTACCCGTAAACTGGTCGTACAATATAAATATCCTGAACATATCGTGGTAATCGCAAAACTCCCGCGAACTACCTCAGGTAAGATACAACAGTTTTTGTTAAGAAAAGATATTCTGCGGCGTTTAACCAAGCATCTCCGTGAAGAAATTGACTAATTTTCATCTTCGGAGCTCACATACCCCCGGCgcataaccgcccncgg                                           >  2:592420/4‑236 (MQ=255)                                                                                                                                                                                         aaGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGat                                       <  1:433928/251‑1 (MQ=255)                                                                                                                                                                                            cgccgcATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTATCCGTAAACATGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCTCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAATAGATATTATGCGGCGTTTATCGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCTCATAACCCCTGCGACTAATCgccgccatgtattta                                    <  2:28100/251‑10 (MQ=255)                                                                                                                                                                                                 catcatTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTGTTTGTTAAGAAAAGATATTAGGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACCAAACGCCGCCGGGGAtttatttta                               <  2:562827/251‑1 (MQ=255)                                                                                                                                                                                                        cagtATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAAGCGAAAAACTACCGCGAACTACCTCAGTTAATATACAAAAGTTTTTGTTAAGAAAAGATATGATGCGGCGTTTAACGCAGGATGTCTGTGAAGATATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTATACGCCGCTTGTGATTTATTTTATttctt                         <  1:637788/247‑1 (MQ=255)                                                                                                                                                                                                                  ggAAGAGGTAGTTGCTTTTTTTAGCCGGAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGTTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATATCCCCGTCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGc                <  2:639658/249‑1 (MQ=255)                                                                                                                                                                                                                      gagGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCt          >  1:327052/1‑251 (MQ=255)                                                                                                                                                                                                                      gagGTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCt          <  2:44811/251‑1 (MQ=255)                                                                                                                                                                                                                      gagGCAGTGGCTTTCTTTAGCCGTAGACAGGTCGCAAAATATAAATATCCTGAACATATCGTGAAAATCGAAAAACCACCGCGCACTACCTCAGGCAAAATAAAAAAGTTTTTGTTAAGAAACGATATTATGCAGCGTTTAACGCAGGATGACTGTGAAGAGATTGAATAAGTTACATCTTCGGGCATCACATAACCTCGGCGACTACACGCCGCCGGGGAGTTATTTTATTTCGTCAGTTTCGCCAGGCt          >  2:73141/1‑251 (MQ=255)                                                                                                                                                                                                                        ggTAGTGGCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCtt         >  1:45774/1‑250 (MQ=255)                                                                                                                                                                                                                              ggCTTTTTTTAGCCGTAAACAGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAgg  >  2:458412/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       GATGAGGCTGGCTATATAAAAATTACCGGACGCAAAAAAGATATTATTGTCCGCGGCGGCGAAAATATTAGCAGCCGTGAAGTGGAAGATATTTTATTGCAGCATCCTAAAATTCACGATGCCTGTGTGGTTGCAATGTCCGATGAACGTTTAGGTGAACGATCATGCGCTTATGTCGTGCTGAAAGCGCCGCATCATTCATTATCGCTGGAAGAGGTAGTGGCTTTTTTTAGCCGTAAACGGGTCGCAAAATATAAATATCCTGAACATATCGTGGTAATCGAAAAACTACCGCGAACTACCTCAGGTAAAATACAAAAGTTTTTGTTAAGAAAAGATATTATGCGGCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGG  >  NC_000913/1784294‑1784765

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