Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,792,541	C→T	G64E (GGA→GAA)	nlpC ←	NlpC/P60 family lipoprotein NlpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,792,541	0	C	T	84.6%	88.1 / 8.5	39	G64E (GGA→GAA)	nlpC	NlpC/P60 family lipoprotein NlpC
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (18/15);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

CATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCC  >  NC_000913/1792308‑1792786                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      cATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGGAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAAGTTGCGTGCCAATAGAGGCATGTTGTTTGGTTTCACGGGGCAGCTTCAAATCGAAACGATCGCTCATCGTCACAACCACAAATTCCGAACAGTCCTCAccg                                                                                                                                                                                                                                      <  2:407749/251‑1 (MQ=255)   ttATCAAGTGAGGACCGCATCACTCCATTGCTGGTAGAGGCGTGGATAAATTGGTTGATGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCAGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGAGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGAACGCGCATCGTCACAACCACAAACTCCGGACAGTCCACAccgcg                                                                                                                                                                                                                                    >  2:308907/1‑251 (MQ=255)                     aTCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCCCCCCGCCTCATGCCACCAt                                                                                                                                                                                                                   >  2:349741/1‑250 (MQ=255)                     aTCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGAAAATACCTCAATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCCTCCTGCCACCAt                                                                                                                                                                                                                   >  2:465481/1‑250 (MQ=255)                       cACTCCCTTGCTAGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTAAATGCAAACCATTTTTTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTTCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGGTGCAAATCGAAACGATCGCGCATCGTCCCAACCACAAATTCCTAACAGTCCACACCGCGACGCGCCATTCCACCATAAc                                                                                                                                                                                                                >  2:291714/1‑251 (MQ=255)                             cTTGCTGGTAGAGGCATCGGTTAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCAATTTGTCCGGCCCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTAACAATTTGCGTGCCAATAAAGGCTTGTTGTTTGGTTTCACGGGGCCGCTGCAAATCGAAACGATCGCGCATCGTCCCCACCACAAATCCCGAACAGTACACACCGCGCCGCGTCCTGCCACCATAACGATACg                                                                                                                                                                                                          >  2:198800/1‑251 (MQ=255)                                   ggTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAGACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGtgc                                                                                                                                                                                                    <  1:350505/251‑1 (MQ=255)                                    gTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCCACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACCAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCCCCATAACGATACGGCGtgcc                                                                                                                                                                                                   >  1:210575/1‑251 (MQ=255)                                                    aTTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGc                                                                                                                                                                                   <  2:80834/251‑1 (MQ=255)                                                             tgGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTAATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTc                                                                                                                                                                          >  1:410522/1‑251 (MQ=255)                                                                    gTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCgg                                                                                                                                                                   >  2:396563/1‑251 (MQ=255)                                                                        aTACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGACGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACAACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAAt                                                                                                                                                               >  1:482077/1‑251 (MQ=255)                                                                                      aCCATTTTGTCCGGACCCCGTTTTGAAAAAGACCATGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTTCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAAt                                                                                                                                                 <  1:294923/251‑1 (MQ=255)                                                                                         aTTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAg                                                                                                                                              <  2:231374/251‑1 (MQ=255)                                                                                                   gACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGAGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGc                                                                                                                                    >  1:445453/1‑251 (MQ=255)                                                                                                     ccccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGGTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCAt                                                                                                                                  >  1:317778/1‑251 (MQ=255)                                                                                                              gAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGTGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACTGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGgcggcg                                                                                                                         <  2:408925/251‑1 (MQ=255)                                                                                                              gAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGGCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATCACGGTAATCGAAACAGAAAGTCTGGCATTTGgcggcg                                                                                                                         >  2:470387/1‑251 (MQ=255)                                                                                                               aaaaaGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGcggcgg                                                                                                                        >  2:195631/1‑251 (MQ=255)                                                                                                                              ccTGGCAGCATCTCGTCTTTACCAATTTGCGTTCTAATTGACGCTTTTTGTCTCGTTTCACTGGGCAGCTTCAAATCGAAACGACCGCGCATCGTCCCAACCACAAATTCCGAACCGTCCACACCGCGCCTCGTCATTCCCCCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAATCCGGCAATAACGGTAATCTAATCAGAAAGTTTGTCATTTGGCGGCGGTGCTTTATGGTGGCt                                                                                                         <  1:408485/251‑1 (MQ=255)                                                                                                                                   cagcaACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACAcc                                                                                                    >  2:204167/1‑251 (MQ=255)                                                                                                                                    agcaACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGTTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACAccc                                                                                                   <  1:76630/251‑1 (MQ=255)                                                                                                                                      caACTCGTCTTTATCAATTTGCGTGCCAATCGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCgg                                                                                                 >  1:129291/1‑251 (MQ=255)                                                                                                                                                          gCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTATGGCATTTGGCGGCGGGGCTTTATGGTGGCTACACCCGGCCAGaaagaaggccgtgatc                                                                             >  2:175305/1‑237 (MQ=255)                                                                                                                                                                            gttgttTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTGTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACg                                                            <  2:317778/250‑1 (MQ=255)                                                                                                                                                                               gttTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCAt                                                         <  1:118214/250‑1 (MQ=255)                                                                                                                                                                                   ggTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCg                                                     <  1:470387/250‑1 (MQ=255)                                                                                                                                                                                                aGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCCGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCCGTGATCCAAATAAAGCGGAAACGCCATCCCTACGGTTCCTCTGt                                       >  2:40232/1‑251 (MQ=255)                                                                                                                                                                                                     cAAATCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGtttttt                                  >  1:434703/1‑251 (MQ=255)                                                                                                                                                                                                        aTCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTTTGGTGGCTACACCCGGCCCGAAGGAGTGCTGGGATCAAAATAAGGCAGAAACGCCTTCCGTACGGtagatcggaagaccggc                               >  2:395942/1‑235 (MQ=255)                                                                                                                                                                                                        aTCGAAACGATCGCGCATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCATAAACGCATTCCGTACGGt                                               <  1:395942/235‑1 (MQ=255)                                                                                                                                                                                                                   cgctcATCGGCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTaa                     <  2:210575/246‑1 (MQ=255)                                                                                                                                                                                                                     cgcATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAAttt                  >  1:109889/1‑251 (MQ=255)                                                                                                                                                                                                                     cgcATCGTCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAAttt                  <  2:12719/251‑1 (MQ=255)                                                                                                                                                                                                                          cGTCACAACCAAAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTTGCTGGTTCTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGGAACACCCGGCCAGAAACAATGCTTTGTTTCAAATAACGCAGAAACGCAATTCGTACTGTTCCTCCGTTTTTTTTTCTTGCATTAATTTCGcgt             >  2:39790/1‑251 (MQ=255)                                                                                                                                                                                                                            tCACAACCACAAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGcgtcg           >  1:239033/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCCNCAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAAtt      <  2:445453/251‑1 (MQ=255)                                                                                                                                                                                                                                    acaAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTAcc   <  1:349741/251‑1 (MQ=255)                                                                                                                                                                                                                                     caAATTCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTAccc  <  2:109889/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      CATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACGGCGTGCCATGCCAGCTTTGTAGCTGGTCGTTCAAACCGGCAATAACGGTAATCGAATCAGAAAGTCTGGCATTTGGCGGCGGTGCTTTATGGTGGCTACACCCGGCCAGAAGCAGTGCTGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCC  >  NC_000913/1792308‑1792786

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