Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	823,045	G→A	P232L (CCC→CTC)	clsB ←	cardiolipin synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	823,045	0	G	A	97.6%	114.8 / NA	41	P232L (CCC→CTC)	clsB	cardiolipin synthase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/18);  minor base C (1/0);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGAGCGGATCGAGATTACTGGACCCTACTGTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGGGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCAGCTCAAACTGGAGAATATCTTCGACAATCGGC  >  NC_000913/822806‑823281                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             tGAGCGGATCGAGATTACTGGACCCTACTGTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCACACCCCGCGCCGAGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGAAAGAAGt                                                                                                                                                                                                                                     >  2:197664/1‑249 (MQ=255)  gAGCGGATCGAGATTACTGGACCCTACTGTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGTGCGGCGGCTGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTa                                                                                                                                                                                                                                    <  2:177922/249‑1 (MQ=255)  gAGCGGATCGAGATTACTGGACCCTACTGTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAgg                                                                                                                                                                                                                                  >  1:348724/1‑251 (MQ=255)            agaTTACTGGACCCGACTGTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGAt                                                                                                                                                                                                                        <  2:168893/251‑1 (MQ=255)                         cTACTGTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCccgc                                                                                                                                                                                                           >  1:284393/1‑251 (MQ=255)                              gTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGAAGGCGTTGGCGATAATCACTTCCCGCCgcgc                                                                                                                                                                                                      >  1:335108/1‑251 (MQ=255)                                cGCCCAGTGATCGTCCATCAATGCCCCTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGGGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGcc                                                                                                                                                                                                     <  1:216825/250‑1 (MQ=255)                                     aGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCCAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCtgagtg                                                                                                                                                                                               >  2:182671/1‑251 (MQ=255)                                            gTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCAttt                                                                                                                                                                                        >  1:254700/1‑251 (MQ=255)                                                   aaTGCCACTTTGCCGTGGAGCGGGGGGCGGCGGTACTCAAAACCCTGAACGCCGCTTTTAACCAGATAGTTATACAGCAAGCGCGCCCCAACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCGGCTTTACGCAAGGCGTGTAAAAATCAATACCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTCCCCGCCGCGCCTGAGTGAGCATTTTCAaata                                                                                                                                                                                 <  1:77551/251‑1 (MQ=255)                                                        cACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGc                                                                                                                                                                            >  1:121292/1‑251 (MQ=255)                                                            ttGCCGTGGAGCGGGCGGCGGCGTTACTCAAAAACCTGAACGCCGCCTTTAACCAGATATTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCCCCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGt                                                                                                                                                                          <  1:351562/249‑1 (MQ=255)                                                            ttGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTc                                                                                                                                                                        >  1:283028/1‑251 (MQ=255)                                                               ccGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTAAACGCCGCCTTTAACCAGATAGTTATCCAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAat                                                                                                                                                                     <  1:288072/251‑1 (MQ=255)                                                                                   gTACTCAAAAACCTGAACGCCGCCTTTACCCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCCCCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCGAATAATGGCGTTCAATATCATCGCGATGTTCTTCGt                                                                                                                                                 <  1:14020/251‑1 (MQ=255)                                                                                     aCTCAAAACCCAGAACGCCGCCTTAAACCAGAGAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTa                                                                                                                                               <  2:103145/251‑1 (MQ=255)                                                                                       ccAAAACCCTGAACGCCGCCTTTACCCAGATAGTTATCCAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCATTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATc                                                                                                                                             <  2:120781/250‑1 (MQ=255)                                                                                          aaaaCCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATcgcg                                                                                                                                          <  1:372092/251‑1 (MQ=255)                                                                                              ccTGAACGCCGCCTTTAACCAGATACTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGCATGATCAGTTTGATCCGCACCCCGTGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGTTGTTCTTCGTTATCGCGcc                                                                                                                                        >  1:318209/1‑249 (MQ=255)                                                                                                            ttAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTg                                                                                                                        >  2:323944/1‑251 (MQ=255)                                                                                                                        ttATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCgg                                                                                                            >  2:212253/1‑251 (MQ=255)                                                                                                                          ataCAGCAAGCGCGCCCCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGTTCCGCCCCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGc                                                                                                           <  1:161986/250‑1 (MQ=255)                                                                                                                              aGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCg                                                                                                      >  2:325241/1‑251 (MQ=255)                                                                                                                                 aaGCGCGCCCCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAGGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGtt                                                                                                   <  1:40993/251‑1 (MQ=255)                                                                                                                                                         ggCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCg                                                                           >  1:189621/1‑251 (MQ=255)                                                                                                                                                                 cGGTTCGCCCTGAATGATAAGTTTGATCCGCCCCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGc                                                                   <  2:79071/251‑1 (MQ=255)                                                                                                                                                                       gCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGCCGTCGCCCCCCGCGTCGTg                                                             >  2:289041/1‑251 (MQ=255)                                                                                                                                                                        cccTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCAATTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGc                                                            >  1:336100/1‑251 (MQ=255)                                                                                                                                                                            gAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCCGCGTCGTGCcgcc                                                        >  1:353440/1‑250 (MQ=255)                                                                                                                                                                                          aTCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGGTGACGGCCCCCCCAACGCCGTGCCCCGCCCTGTCcaggccg                                          >  2:176271/1‑249 (MQ=255)                                                                                                                                                                                           tCCGCACACCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGATATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCGCGTGGTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCcaggcagg                                         >  1:275103/1‑251 (MQ=255)                                                                                                                                                                                                          gTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCAGCTCAAAc                          <  1:332211/251‑1 (MQ=255)                                                                                                                                                                                                           tGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGCGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATatcatc                                                                                                                                                               >  2:346418/1‑117 (MQ=255)                                                                                                                                                                                                           tGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAAGCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATatcatc                                                                                                                                                               <  1:346418/117‑1 (MQ=255)                                                                                                                                                                                                                         aGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCCGCTCAAACTGGAGAATATCTTCg           >  1:373003/1‑251 (MQ=255)                                                                                                                                                                                                                         aGGCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCCTTGTGCTGACGTCGCCCCCCGCGTCGTGCCGCGCTTTTTCCAGGCAGGTTTTCCAGCCCCAACTGGAGAATATCATCg           >  2:311113/1‑251 (MQ=255)                                                                                                                                                                                                                          ggCGTGTAAAAATCGAAGGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCAGCTCAAACTGGAGAATATCTTCGa          <  2:16575/251‑1 (MQ=255)                                                                                                                                                                                                                           gCGTGTAAAAATCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCAGCTCAAACTGGAGAATATCTTCGAc         <  1:371833/251‑1 (MQ=255)                                                                                                                                                                                                                                tAAAAATCGCTAGCCGAGGAAGAAGTAGGCGTTGGCGATACTCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGCTGACGTCGCCCCCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCACCTCAAACCGGAGAATATCTTCGACCATCg    >  2:26227/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaaaaTCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCTCTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCATGTTTTCCAGCTCAAACTGGAGAATATCTTCGACAATCg    <  2:121627/250‑1 (MQ=255)                                                                                                                                                                                                                                  aaaaTCGATAGCCGAGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCAGCTCAAACTGGAGAATATCCTCGACAATCGGc  >  1:144172/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             TGAGCGGATCGAGATTACTGGACCCTACTGTCGCCCAGTGATCGTCCATCAATGCCACTTTGCCGTGGAGCGGGCGGCGGCGGTACTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGTGTAAAAATCGATAGCCGGGGAAGAAGTAGGCGTTGGCGATAATCACTTCCCGCCGCGCCTGAGTGAGCATTTTCAAATAATGGCGTTCAATATCATCGCGATGTTCTTCGTTATCGCGCCAGACCAGCAATACCTGCGCTTCTCCCGGCTGGCGGTTCTCTTCCGCTTTGTGATGACGTCGCCACCAGCGTCGTGCCGCGCTCTGTCCAGGCAGGTTTTCCAGCTCAAACTGGAGAATATCTTCGACAATCGGC  >  NC_000913/822806‑823281

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