Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,040,126	G→T	L277L (CTG→CTT)	appB →	cytochrome bd‑II subunit 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,040,126	0	G	T	97.7%	144.4 / NA	43	L277L (CTG→CTT)	appB	cytochrome bd‑II subunit 2
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (20/22);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGGAGCT  >  NC_000913/1039880‑1040367                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  gATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCCGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGGCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACgg                                                                                                                                                                                                                                               >  1:228455/1‑251 (MQ=255) gATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACgg                                                                                                                                                                                                                                               <  1:51229/251‑1 (MQ=255)       tCTGCGTTCACAGCTGGCGCCCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCCCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGGGGCCGGGCTACCCGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTTTGGGTCTTCCCGCTGCTGGGATTCTTCTGCCCATTTCTTACGGTGATgg                                                                                                                                                                                                                                         >  2:287696/1‑251 (MQ=255)             ttCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTAATGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTCCCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGAttt                                                                                                                                                                                                                                   >  1:447300/1‑251 (MQ=255)                              aaCGCGCTGCACTGCTGGTGATGCTGTGCTTTTTGCTGGCGGGATACTGGCTGGGGGTAGGTATTGATGGCTTTGTACTGCTCGCCCAGGGTGCTAATGGTCCTTCCAATCCGTTAATGAAACTGGGGGCAGTGCTCCCTGGTGCCTGGGTGAATAATTTTTTCGAGTCGCCCGTTTTGTGGATCTTTCCGCTGCTGGGGTTCTTTTGCCCATTGCTTACGGTGATGGCGATTTGTCGTGGTCGCCCgggg                                                                                                                                                                                                                  >  2:285546/1‑250 (MQ=255)                                           gTTGGTGATGCTGTGCTTTTTCCTGGCGGGGTACTGGCAGTGTGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACAGTTCTTCCAATCCGTTAATGAAACTGGTGGCATTTCTACCTGGGGCCTGGATGAATAATTTTGGCGAGTCGCCCGTGTTGTGGATCCTGCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTg                                                                                                                                                                                                      <  1:453680/250‑1 (MQ=255)                                              ggTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGCTGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGGGATGGCGATTTATCGTGGTCGCCCGGGGTGGGGATTTTTGATgg                                                                                                                                                                                                  >  2:175490/1‑251 (MQ=255)                                              ggTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATg                                                                                                                                                                                                   <  1:45719/250‑1 (MQ=255)                                                      tgtgCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGGTGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCAttt                                                                                                                                                                                          >  1:234095/1‑250 (MQ=255)                                                       gtgCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTg                                                                                                                                                                                          <  2:228455/250‑1 (MQ=255)                                                                        gTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGAtt                                                                                                                                                                        <  2:333976/251‑1 (MQ=255)                                                                           aCTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTc                                                                                                                                                                     <  1:90417/251‑1 (MQ=255)                                                                             tGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCAc                                                                                                                                                                   >  1:345474/1‑251 (MQ=255)                                                                             tGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCAc                                                                                                                                                                   >  2:356021/1‑251 (MQ=255)                                                                               gCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAAGCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTGCTGCCCAGTGCTTACTGTGATTGCGATTTATCGGGGTCGCCCGGTTTGGGGATTTTTGATTTCATCACTGATGCAATTTGTCGTGTTGTTCACgg                                                                                                                                                                 <  2:298963/251‑1 (MQ=255)                                                                                         ggTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCAc                                                                                                                                                       <  1:287696/251‑1 (MQ=255)                                                                                              tGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTg                                                                                                                                                   <  1:111115/250‑1 (MQ=255)                                                                                                aTGGCTTTGTACTGCGCGCCCAGGATGCTAACGGTCCTTCCAAGCCGTTAATGAAACTGGTGGCAGTTCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTGATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGGTTCACGCTGTTc                                                                                                                                                <  2:241534/251‑1 (MQ=255)                                                                                                  ggCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTccc                                                                                                                                              <  1:356021/251‑1 (MQ=255)                                                                                                             gCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCAt                                                                                                                                     <  1:302069/249‑1 (MQ=255)                                                                                                                cGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGt                                                                                                                                >  1:403185/1‑251 (MQ=255)                                                                                                                      ggATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCg                                                                                                                          >  2:192998/1‑251 (MQ=255)                                                                                                                        aTGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTg                                                                                                                        <  1:334016/251‑1 (MQ=255)                                                                                                                               cGGTCCTACCAATCCGATACTGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGCTCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCg                                                                                                                  >  2:280704/1‑250 (MQ=255)                                                                                                                                        cAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGcc                                                                                                        >  2:152255/1‑251 (MQ=255)                                                                                                                                              gTTAATGAAACTGGTGGCAGTGCTACCTGGAGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCACTTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGt                                                                                                  >  1:42031/1‑251 (MQ=255)                                                                                                                                                tAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGACGCCCGTTTTGTGGATCTTCCCGCTGCTGGGCTTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGAAGTtg                                                                                                >  1:65003/1‑251 (MQ=255)                                                                                                                                                  aTGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTtgtg                                                                                              <  1:449844/251‑1 (MQ=255)                                                                                                                                                    gAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTggg                                                                                            >  1:273278/1‑251 (MQ=255)                                                                                                                                                               cAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCACGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTcc                                                                                 >  1:196594/1‑251 (MQ=255)                                                                                                                                                                   gCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTc                                                                             <  1:311731/251‑1 (MQ=255)                                                                                                                                                                           gTGCCTGGATGAATAATTTTGTCGCGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCGGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACg                                                                     <  2:403185/251‑1 (MQ=255)                                                                                                                                                                             gCCTGGATGAAGAATTTTGTCGAGTCGCCCGTTTTGTGGATCTGCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCTAGTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCAGGCCGGCAAGCGTGAGGCCGATCTCCAGCCTGACGTTTTGGGACAGGACTTCCAGTCAGCTGACGCt                                                                   <  2:218966/251‑1 (MQ=255)                                                                                                                                                                              ccTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTg                                                                  <  2:345474/251‑1 (MQ=255)                                                                                                                                                                                        taatTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGGGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGCTTTTTGATGGCATCATTGATGCACTTTGGCGTGCTTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGt                                                        >  1:47791/1‑251 (MQ=255)                                                                                                                                                                                          atTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTg                                                      <  1:318621/251‑1 (MQ=255)                                                                                                                                                                                          atTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCGTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCTCGCTGTTCCCCTTTGGCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTg                                                      <  2:9968/251‑1 (MQ=255)                                                                                                                                                                                                    gTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGc                                            >  1:54714/1‑251 (MQ=255)                                                                                                                                                                                                                  ggATCTTCCCGTTGCTGGGATTATTCTGCCCATTGCTTACGGTGATGGCGAGTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGGCCGATCTCCAGCTTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGcc                               <  2:316426/250‑1 (MQ=255)                                                                                                                                                                                                                           cgttgctgggttttctctGCCCAGTGCTGAGGGTGTTGGCTATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTTCTTCCAGTCAGCTGACTCTTAGTATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTg                     <  2:315630/236‑1 (MQ=255)                                                                                                                                                                                                                            gctgctGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCAACACGCTGTTCCACTTTGTCATGCCGCCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTAATTCCAGTCAGCAGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGc                    >  1:463220/1‑251 (MQ=255)                                                                                                                                                                                                                                  gggATTCGTCTGCCCATTTCTTGTGGTGATGGCGATTTATTGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGAGCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATGATGTTGGTAAGCGTGCTGATATTTTTGCCTATTGTGTTGCTCTac               <  2:392832/250‑1 (MQ=255)                                                                                                                                                                                                                                              cccATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGCTCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGCGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGGAGct  >  2:189659/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  GATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGGAGCT  >  NC_000913/1039880‑1040367

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