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.8%	155.7 / NA	46	L277L (CTG→CTT)	appB	cytochrome bd‑II subunit 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (19/26);  minor base A (1/0);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGG  >  NC_000913/1039880‑1040363                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              gATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACgg                                                                                                                                                                                                                                           >  1:126512/1‑251 (MQ=255)              tCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGAttt                                                                                                                                                                                                                               <  1:106798/250‑1 (MQ=255)                          aCCAACCGCGCTGCACTGTGGGGGATGCTGTGCTTTTTTCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGGTAATGAAACTGGGTGCAGTGCTACCGGGTGCCGTGATGAATAATTTTGTCGAGTCTCCCGTTTTG‑TGTATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGccc                                                                                                                                                                                                                  <  1:447924/251‑1 (MQ=255)                                       cACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCCGGGTGCTAACGGTCCTTCCAATCCGTTTATGAAACTGGTGGCAGTGCTACCTGGGGCCTGGGTGAATAATTTTGTCGGGTCGCCCGGTTTG‑TGGATCTTCCCGCTGCTGGGGTTCTTCTGCCCCTTGCTTACGGTGATGGGGATTTATCTTGGGCGCCCGGGTTGGGGGttt                                                                                                                                                                                                     >  2:214586/1‑251 (MQ=255)                                          tGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTTATGGCTTTGTACTGCTCGCCCCGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGGTTTATCGTGGTCGCCCGGGTTGGGGGTTTTTg                                                                                                                                                                                                  >  2:202705/1‑251 (MQ=255)                                              ggTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAAGCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATgg                                                                                                                                                                                              <  1:87939/251‑1 (MQ=255)                                                 gATGCTGTGCTTTTTGCTGGAGGGATAATGGCTGTGGGTCGGGATTGATGGCTTTTTACTGGTCGCCCAGGGTGCTAAAGGGTCTTCCAATACGTTAATGAAACTGGTGGCAGTGCTACCTGGGGCCTGGATGGATAATTTTGTCGAGTCGCCCGTTTTGGTGG‑TCTTCCCGCTGCTGGGGTTCTTCTGGCCATTGCTTACGGGGATGGCGGTTTATCTTGGTCGCCCGGGGTGGGGGTTTTTGATGGcat                                                                                                                                                                                           >  2:100577/1‑251 (MQ=255)                                                              ttGCTGGCGGGTTACTGGCTGTGGGTCGGTATTTATTGCTTTGTCCTGCTCGCCCAGGATGCTAACGTTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTGTGATGGCATCATTGATGCAAt                                                                                                                                                                               <  2:84419/250‑1 (MQ=255)                                                               tGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAAtt                                                                                                                                                                              <  1:333563/250‑1 (MQ=255)                                                                   ggCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATGTCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTTGGGATTTTTGATGGCATCATTGATGCAATTTgg                                                                                                                                                                           <  2:270377/249‑1 (MQ=255)                                                                    gCGGGTTCCTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCACGCTCAGGATGCGAACGGTCCTTCCAATCCGTTAATGAAACGTGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATGTGGCGt                                                                                                                                                                        <  2:53676/251‑1 (MQ=255)                                                                    gCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGt                                                                                                                                                                        <  1:33927/251‑1 (MQ=255)                                                                    gCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGGGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGt                                                                                                                                                                        <  1:202705/251‑1 (MQ=255)                                                                       ggTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCACCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGAt                                                                                                                                                                     <  1:9847/251‑1 (MQ=255)                                                                                               gATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAAttt                                                                                                                                                                             >  1:268306/1‑219 (MQ=255)                                                                                               gATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAAttt                                                                                                                                                                             <  2:268306/219‑1 (MQ=255)                                                                                                    cTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCt                                                                                                                                        >  1:69082/1‑251 (MQ=255)                                                                                                        gTACTGGTCGCCCAGGATGCTACCGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCGACCTGGTGACTGGATGAATAAGTTTGTCAAGTAGCCCGTTTTG‑TGGATCTTACCGCTGCTGGGATTCTTCTGCCCCTTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGATTGGGGCTTTTTGATGGCATCATTGATGCAATTTGGAGTGATTTTCAGGGCAGACAttaccgatttcgcctttgt                                                                                                                                    >  2:177493/1‑233 (MQ=255)                                                                                                           cTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCAt                                                                                                                                 >  1:163554/1‑251 (MQ=255)                                                                                                              atCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGACGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGTTTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGcc                                                                                                                              >  2:247652/2‑251 (MQ=255)                                                                                                                       gATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGt                                                                                                                     >  2:438380/1‑251 (MQ=255)                                                                                                                                ggTCCTTCAAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGAGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGCTCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGGTTTATCGGGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCCCGCTGTTCCCCTTTGTCATGCCGTCACGCGTGAGTCCGAt                                                                                                            >  2:258847/1‑251 (MQ=255)                                                                                                                                  tCCTTCCAATCCGTTAATGCAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTCCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGTATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGAtct                                                                                                          >  2:435545/1‑251 (MQ=255)                                                                                                                                   ccTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGAtct                                                                                                          <  2:69082/250‑1 (MQ=255)                                                                                                                                       ccAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAg                                                                                                      <  1:80484/250‑1 (MQ=255)                                                                                                                                       ccAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGc                                                                                                     >  1:79851/1‑251 (MQ=255)                                                                                                                                                tAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGtt                                                                                             <  1:332957/250‑1 (MQ=255)                                                                                                                                                       aCTGGTGGAAGTGCTACATGGTGACTGGAGGAATAATTTTGACGAGTCGCACGATTTG‑TGGATCTTCGCGCTGCTGGGATTATTCTGCCCATTGCTAACGGTGAGGGCGATTTATCGTGGTCGCCCGGGTTGGGGCTTTTTGATGGCATCATTGATGCAATTTGGCGTGATGTGCCAGGCAGGCATAACGCCGTTCCCCTATGTAAAGACGTCAAGCGTGAGTCCGATCTCCACGCTGACGCTGTGGGAcg                                                                                     >  2:232144/1‑250 (MQ=255)                                                                                                                                                        cTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGCTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGAc                                                                                      >  1:253811/1‑249 (MQ=255)                                                                                                                                                          ggtggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGt                                                                                   <  2:89564/250‑1 (MQ=255)                                                                                                                                                          ggaggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGtc                                                                                  >  2:144048/4‑250 (MQ=255)                                                                                                                                                           gtggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTAc                                                                                 >  2:164451/1‑251 (MQ=255)                                                                                                                                                             ggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑AGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGGGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACtt                                                                               >  1:227456/1‑251 (MQ=255)                                                                                                                                                                      aCCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAg                                                                       <  1:348645/250‑1 (MQ=255)                                                                                                                                                                          ggTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGAc                                                                  <  2:156991/251‑1 (MQ=255)                                                                                                                                                                             gCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCt                                                               <  1:221881/251‑1 (MQ=255)                                                                                                                                                                                       ataatTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTGGATGGCATCATTGATGCAATTTGGCGGGTTTTTCACGGGAGGCATCACGCTGTTCCCCTTTGTCATGCCGGCAAGCGTGAGTGCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATg                                                     <  1:59972/251‑1 (MQ=255)                                                                                                                                                                                                           gTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATAtt                                 <  2:79851/251‑1 (MQ=255)                                                                                                                                                                                                           gTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATAtt                                 <  1:214586/251‑1 (MQ=255)                                                                                                                                                                                                                g‑tgGATCCTCCCGCTGCTGCTCTACTTCTGCCCATTGCTTCCGGTGCTGGCGATTTATCGTGGTCGCCCGGGTTGGGGCTTTTTGATGGCATCATTGCTGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCCGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCAGACGCTGAGCATTATGTTGGTAATCGTGATGATATTTTTGc                            >  2:50886/1‑251 (MQ=255)                                                                                                                                                                                                                   ggATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGccc                          <  2:335521/251‑1 (MQ=255)                                                                                                                                                                                                                                gctgGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATTGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGctc              <  1:38835/250‑1 (MQ=255)                                                                                                                                                                                                                                gctgGGATTCTGCTGCCCAGTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCCGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGctct             <  2:60307/251‑1 (MQ=255)                                                                                                                                                                                                                                   gggATTCTTCTGCACATTGCTTACGGTGATTGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCAGCATTGATGCAATTTGGCGTGATTGTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGGCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGGCAGCTGACGCTGAGCATGATGTTGGTAATCGTGCTGATATTATTGCCCATTGTTTTGCTCTac           <  2:431032/250‑1 (MQ=255)                                                                                                                                                                                                                                    ggATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGCTCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTACTCGTGCTGATATTTTTTCCCATTGTGTTGCTCTacac         >  2:21933/1‑251 (MQ=255)                                                                                                                                                                                                                                              gCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTGTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTgg  <  1:258847/248‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              GATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTG‑TGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGG  >  NC_000913/1039880‑1040363

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