Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	822,527	C→T	E405K (GAA→AAA)	clsB ←	cardiolipin synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	822,527	0	C	T	97.6%	127.6 / NA	42	E405K (GAA→AAA)	clsB	cardiolipin synthase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/22);  minor base G (0/1);  total (19/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

CATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTATCGCGCAGCGTCTGGTTAAAATGACGATCGT  >  NC_000913/822294‑822769                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    cATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCTGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCATCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGtt                                                                                                                                                                                                                                     <  1:510792/250‑1 (MQ=255)    agcagcCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTcc                                                                                                                                                                                                                                  <  2:470752/250‑1 (MQ=255)     gcagcCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTTTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCACCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCAt                                                                                                                                                                                                                                <  1:461315/251‑1 (MQ=255)           gTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATGGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACGCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCgg                                                                                                                                                                                                                          <  1:502672/251‑1 (MQ=255)                 cataTAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTTCCCCCGTTTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTtgtgt                                                                                                                                                                                                                    <  2:18268/248‑2 (MQ=255)                 cAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTtgtg                                                                                                                                                                                                                     <  2:204834/250‑1 (MQ=255)                             gacCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCaa                                                                                                                                                                                                         >  1:167829/1‑250 (MQ=255)                                    aGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGCAGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACACTCACCGCGATAAACACCCGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGCGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCCACCTGCGcc                                                                                                                                                                                                 >  2:456229/1‑251 (MQ=255)                                         gACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACg                                                                                                                                                                                            >  1:522339/1‑251 (MQ=255)                                            cGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTgg                                                                                                                                                                                         >  1:444564/1‑251 (MQ=255)                                                             cgcgATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCTCTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTTCGCCAGACGTGGCGTGTGTGCCGGAa‑gcc                                                                                                                                                                       <  2:374274/251‑1 (MQ=255)                                                               cgATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAa‑gccag                                                                                                                                                                     <  1:487486/251‑1 (MQ=255)                                                                          gcgGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGAGATTTACTAATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCAATTGTCGGTTGTGCGGGCGGCTCAAACTGCGCCAGACCTTGCGTGTTTGCCAGAA‑GCCAGCCAACCCgcgc                                                                                                                                                          >  2:379979/1‑251 (MQ=255)                                                                                      ccAGACCTCTTCCCAGTCCACTTTTTTGGCTTAGACCACCAGCAACACAATCACCGGGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCCAAGCGCCTGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCAGTGTCGGTTGTGCAGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGc                                                                                                                                               <  1:611172/250‑1 (MQ=255)                                                                                      ccAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGACAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGCAGGGGCg                                                                                                                                              >  1:357957/1‑251 (MQ=255)                                                                                                 cccAGTCAACTTTCTTGACGAAGACCACCAGCAACAAAATCACCGCGATAAAAAACAGCCCGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCCGAC‑GACAGCCAACCAACGCCGGGAAGTGGCGGAAAAAgggg                                                                                                                                    >  2:325896/1‑247 (MQ=255)                                                                                                     gTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGCTAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGc                                                                                                                               >  2:419542/1‑251 (MQ=255)                                                                                                       ccACTTTTTTGGCGTAGAACACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCGTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCa                                                                                                                             <  2:435589/251‑1 (MQ=255)                                                                                                        cgcggtttgggCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTGTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCAGTGTCGGTTGTGCGGGCGGAGCAACCTGCGCCAGACGGGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGGACGCCAg                                                                                                                            <  2:260189/242‑1 (MQ=255)                                                                                                                               gCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCATACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTccacc                                                                                                     >  1:432927/1‑251 (MQ=255)                                                                                                                                          cccGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTg                                                                                          <  1:244165/250‑1 (MQ=255)                                                                                                                                          aCCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTg                                                                                          <  1:398108/251‑1 (MQ=255)                                                                                                                                                         aGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCATGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGCACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTc                                                                            >  1:218399/1‑250 (MQ=255)                                                                                                                                                                 gAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGGGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGTCGTAAAAAGTGGAACGCCAGCACTCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACct                                                                    <  1:45998/250‑1 (MQ=255)                                                                                                                                                                 gAGGATCTGCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTGCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACctg                                                                   <  1:606235/251‑1 (MQ=255)                                                                                                                                                                        ttcttTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGGGCCGGAA‑GCCAGCCCACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGGTCCACCAGGTGCGTTTGGGCAGCATGGTTTCCTCCACCTGCTGACa                                                             >  2:511598/1‑250 (MQ=255)                                                                                                                                                                               cTAAGCGCCCGAGCGGGTGTGATTTCCTCATCCGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCACCCTGCGCCAGACGTGGCGTGTGTGGCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGGGCGTTTGGGCCGCATGGTTTCATCCACCTGCTGTCAATCTGCgg                                                     >  2:51280/1‑251 (MQ=255)                                                                                                                                                                                tAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGGGTGTGCCGGGA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGc                                                    >  1:231458/1‑251 (MQ=255)                                                                                                                                                                                 aaGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCgg                                                     <  1:367150/249‑1 (MQ=255)                                                                                                                                                                                  aGCGCCGGCGAGGGTGTGCTTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCAGATCCTGCGTTTCCATGGTCGGGTGTGCGGGCGGGTCAACCTGCGCCAGACGGGGCGTGTGTGCCGGAAAGCC‑TCCAAACAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGAACGCTTTTGGGCAGGTTCCACCAGGGGCATTTGGGTAGCAATGCTTCAGCCACCTGCTGACAATCTGCGGCa                                                   >  2:173536/1‑250 (MQ=255)                                                                                                                                                                                  aGCGCCAGCGCGGGTGTGCTTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCCGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGGGCGTTTGGGCAGAATGGTTTTATCCACCTGCTGACAAACTGCGGccc                                                  >  2:394053/1‑249 (MQ=255)                                                                                                                                                                                        aGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACACTCTGCGGCAATAATGc                                            >  2:193434/1‑251 (MQ=255)                                                                                                                                                                                            cgGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGtt                                        <  1:620330/251‑1 (MQ=255)                                                                                                                                                                                             gggTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGATGTGCGCGCGGCTCAACCTGCGCCAGACGTGTCGTGTGTGCCGGAA‑GCCAGCCCACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCCCGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCAGGGTTTCCTCCACCTGCTGGCAATCTGCGGCAAAAATGCCGTTc                                       >  1:607502/1‑251 (MQ=255)                                                                                                                                                                                              ggTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTTCGGCAATAATGCCGTTCa                                      <  1:150484/251‑1 (MQ=255)                                                                                                                                                                                                   gATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTTGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTa                                 <  2:12034/251‑1 (MQ=255)                                                                                                                                                                                                   gATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTGTCCCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTa                                 <  2:498476/251‑1 (MQ=255)                                                                                                                                                                                                                     ttACCCCCGTGTTTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTATCGCGCAGCGTCTGGtt                >  1:39646/1‑250 (MQ=255)                                                                                                                                                                                                                          cccGTGTCTTCAGTTTTTACCCGATCCTGCGTTTCCATTGTCTGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGTAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTATCGCGCAGCGTCTGGTTAAAATg          <  1:166010/251‑1 (MQ=255)                                                                                                                                                                                                                                ttttCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTATCGCGCAGCGTCTGGTTAAAATGACGATc    >  1:620245/1‑251 (MQ=255)                                                                                                                                                                                                                                  ttCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTATCGCGCATCGTCTGGTTAAAATGACGATCg   <  1:193434/250‑1 (MQ=255)                                                                                                                                                                                                                                   tCAGTTTTTACCCGATCCTGCGTTTCCATTGTCGGGTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTATCGCGCAGCGTCTGGTTAAAATGACGATCGt  <  1:394053/250‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    CATAGCAGCCGTAAATCAGATAGCTGACGACCACCAGCCCGACCGCACTAAGCAGCGCAACGCGATTGTAGTCGCGGATGACCTTCCAGACCTCTTCCCAGTCCACTTTTTTGGCGTAGACCACCAGCAACACAATCACCGCGATAAAAAACAGCCAGGTGAGGATCTTCTTTGCTAAGCGCCAGCGCGGGTGTGATTTACTCATCAGGGTTTTACCCCCGTGTTTTCAGTTTCTACCCGATCCTGCGTTTCCATTGTCGGTTGTGCGGGCGGATCAACCTGCGCCAGACGTGGCGTGTGTGCCGGAA‑GCCAGCCAACCAGCGCCGGGAAGTGGCGTAAAAAGTGGAACGCCAGCACGCTTTTGGTCAGGTTCCACCAGGTGCGTTTGGGCAGCATGGTTTCATCCACCTGCTGACAATCTGCGGCAATAATGCCGTTCAGATTATCGCGCAGCGTCTGGTTAAAATGACGATCGT  >  NC_000913/822294‑822769

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