Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	968,837	G→A	D158N (GAT→AAT)	lpxK →	tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,837	0	G	A	97.4%	126.7 / ‑6.0	38	D158N (GAT→AAT)	lpxK	tetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/22);  minor base C (1/0);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAAC  >  NC_000913/968603‑969056                                                                                                                                                                                                                                           |                                                                                                                                                                                                                              ggggTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGt                                                                                                                                                                                                               <  1:344887/251‑1 (MQ=255)    gTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAt                                                                                                                                                                                                            <  1:62301/251‑1 (MQ=255)       gTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTTATGCGGTAAAAGCCCTTCTTGCGCCACACCCTGGTGTGCCGATCCTCGTAACCGCAGACGGGTTACAGCCTTTCCGTCTGGCGCGTAATGTGGGAATTGTCGTTAtttt                                                                                                                                                                                                         >  1:482126/1‑249 (MQ=255)                    atGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGAGTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCg                                                                                                                                                                                             <  1:230141/250‑1 (MQ=255)                                  gAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGc                                                                                                                                                                                <  1:328669/249‑1 (MQ=255)                                     tCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCtggtg                                                                                                                                                                           >  1:111374/1‑251 (MQ=255)                                                         ggCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTTCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAACCCTTCTGGCGCAACACCCTGGTGTGCAGATCATCGTAACCGAAGACGGTTTACCGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGGGTGCGGCGCTTTGGGAATGGCTGGGGGGTGCCGGCGGGGGCaaatg                                                                                                                                                       >  1:321315/1‑248 (MQ=255)                                                              aTACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgag                                                                                                                                                  <  2:114992/251‑1 (MQ=255)                                                                           acaGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATGGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCt                                                                                                                                     <  1:337088/251‑1 (MQ=255)                                                                            caGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTTCCGGCGGGGGCGATGCGTGAGCGGGCGGGGGGCtt                                                                                                                                    >  2:15853/1‑251 (MQ=255)                                                                                  ggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTc                                                                                                                              <  2:321315/251‑1 (MQ=255)                                                                                   gTGATGAACCTGTGTTGATGTATCAACGCACTGATGCGCCTGTTGCGGTTTCGCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTGAAAGTCg                                                                                                                             <  2:66063/251‑1 (MQ=255)                                                                                     gatgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGt                                                                                                                           <  1:370094/251‑1 (MQ=255)                                                                                          aCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATg                                                                                                                      <  2:334317/251‑1 (MQ=255)                                                                                           ccTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGc                                                                                                                     >  1:166123/1‑251 (MQ=255)                                                                                              gtgtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGt                                                                                                                  <  1:15853/251‑1 (MQ=255)                                                                                                      ttATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGAGTGAGAGAGCGGGGCGCTTAAAGTCGGTTGAGGCGGTAATCGtt                                                                                                           >  1:301925/1‑249 (MQ=255)                                                                                                           aaCGCACGGATGCTCCTGGGGCGGTTTCGCCGGTTCGGTTGGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCTGTTGATGCGGTAATCGTCAAcggc                                                                                                     <  2:160631/251‑1 (MQ=255)                                                                                                                  tGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTTAAGTCGGTTGATGCGGTAATCGTCAACGGGGGGGT‑cc                                                                                             >  1:197709/1‑251 (MQ=255)                                                                                                                      gcgcCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCg                                                                                         <  1:183815/251‑1 (MQ=255)                                                                                                                            gTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGAGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAATCGTCAACGGGGGTGG‑CCCTCGCCGCgg                                                                                   >  1:294309/1‑251 (MQ=255)                                                                                                                                        cacgttcgttCTGATGCGGTAAAAGCCATGCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCa                                                                        <  2:380010/248‑1 (MQ=255)                                                                                                                                                    gATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGcc                                                          <  1:222137/251‑1 (MQ=255)                                                                                                                                                          gTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGAAATCGTCAACGGCGGTGT‑CCCTCCCACCGGTGAAATCCCCATGCATCTA‑CCGCCggggc                                                     >  1:347932/1‑248 (MQ=255)                                                                                                                                                                  cATTCTGGCGCAACACCCTGATGTTCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGa                                            <  2:368229/251‑1 (MQ=255)                                                                                                                                                                  cATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGGAATCGTCAACGGCGGGGT‑CCCTCGCAGCGGTGGAATCCCCAAGCATCTG‑CTGcccgggccgggggtgg                                            >  1:409630/1‑237 (MQ=255)                                                                                                                                                                    ttCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGT‑CACTCCCAGCGGTGGAATCCCCAAGAATCTG‑CTGCCGGGGGAGGCggggtat                                          >  1:461029/1‑246 (MQ=255)                                                                                                                                                                       tGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAACGTCGGTTGATGCGGTAATCGTCAAAGGCGGTGT‑CCCTCGCAGCGGTGGAATCCCCCTGCATCTG‑CTGCCGGGGCAGGCgggggatttt                                       >  1:485671/1‑243 (MQ=255)                                                                                                                                                                         gcgcAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGAATTTACg                                     <  1:35495/251‑1 (MQ=255)                                                                                                                                                                             aaCACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCCACGGCGGGGT‑CCCTCGCCACGGGGAAATTCCCATGCATCTG‑CTGCCGGGGCAGGCGGGGGATTTACGTaac                                 >  1:453327/1‑249 (MQ=255)                                                                                                                                                                                 cccTGAGGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTGCCGTCTGGCTCGTAATGTGGAAATTGTGGTTATGGATGGTGGGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGAATTTACGTACCGGTa                             <  2:140667/251‑1 (MQ=255)                                                                                                                                                                                        gtgCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGtt                       <  1:466856/250‑1 (MQ=255)                                                                                                                                                                                         tgCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGGAATCCCCCTGCCTCTG‑CTGCCGGGGCCGGCGGGGGATTTTCGTCACGGTACGCGTtgg                     >  1:221770/1‑250 (MQ=255)                                                                                                                                                                                          gCAGAGCATCGTAACCGACGACGGTGTACAGCATTACCGTCTGTCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTTGGTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAATTCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCTTtgtg                    <  2:242935/251‑1 (MQ=255)                                                                                                                                                                                            aGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTCATGAGGCAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGAGTGGGCGGGGGGGGGGCTTAAATGCGGTTGATTCGGTAATCGTAAAAGGCGG‑GTACCCCCCCAGCAGGGGAATCCCCATGCATCTGGCT‑CCGGGGCAGGGGGTGGAATTACGTACCCGGAcgcgcttggac                  >  2:368427/1‑244 (MQ=255)                                                                                                                                                                                                catcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGtt               <  2:166123/250‑1 (MQ=255)                                                                                                                                                                                                      aaacgacgaCGGTTTACAGCATGTCCGTCTGGCGCTTAATGTGGAAATTGTCGTTAGTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAAGCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGGCAGGCGGGGAAGTTACGTACCGGTACGCGTTGTGACGtttatcag         <  2:486246/247‑7 (MQ=255)                                                                                                                                                                                                            cgacgGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGGTGATGCGGTAATCGGCAACGGCGGGGT‑CCCCCGCAGCGGTGAAATCCCCATGCTTTTG‑CTGCCGGGGCAGGGGGTGAAATTTCCTACCGGAGCGCGGTGTGACGTTGCTCAGTTTggac  >  1:118216/1‑248 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                              GGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTG‑CTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAAC  >  NC_000913/968603‑969056

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: ‑