Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,722,946	G→A	A74T (GCT→ACT)	pssA →	phosphatidylserine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,722,946	0	G	A	94.6%	102.3 / NA	40	A74T (GCT→ACT)	pssA	phosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (16/19);  minor base C (2/0);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTG  >  NC_000913/2722711‑2723191                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      gaagaaATGCACTGTGATGTTGTCAAAATTTAAGCGTAATACACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGATGATGATGTCGATTTCTTTTACGCTCCCTCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGCACAGGATGACGGTGGCAAAGGCCTTCTGACCGCGTTGTATGAGACGAAAAGGCAGCGTc                                                                                                                                                                                                                                        >  2:36897/1‑251 (MQ=255)           aCTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGa                                                                                                                                                                                                                              <  1:376923/251‑1 (MQ=255)                    gtgataaaaTTTAAGCGTAATAAACATCAACAACACCTTGACCAACTACCCAAGATTTCACAATCAGATGATGATGGCATTTTCTGTTACGCGGCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGt                                                                                                                                                                                                                     <  2:257722/246‑1 (MQ=255)                           aaTTTAAGCGTAATAAACATCAACAACACCTTGCCCACCTACCCAATATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGTATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTc                                                                                                                                                                                                              <  1:309772/251‑1 (MQ=255)                                aaGCGTAATAAACATCAACAACACCTTGCCCAACTACACAAGATTTCACAACAAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCTGAACTGGATGTGCGGGTGCTGGTCGACt                                                                                                                                                                                                          <  1:549055/250‑1 (MQ=255)                                       ataaACATCAAAAACACCTTGCCCCACTACCCAAGATTTCTCAATAAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCCGCGCATTTGCATTGTCGCCCTGTATCTCGAACCGGCTGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCGGCGTCCGAAACTGGATGAGCGGGTGCTGGTCGACTGGAATc                                                                                                                                                                                                    >  2:425317/1‑249 (MQ=255)                                        aaaacacaaacacccccTTGCCCAACTCCCCAAGATTTCGCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTTGCATCGTg                                                                                                                                                                                                 <  2:350675/236‑1 (MQ=255)                                         aacaataaacaacCCCTGGCCCAACTCCCCAAGATGTCTAAATCGGTTGATGATTTCGATTTCTGTTACGCTCCCGCCGTCTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCTCATTTGCATTTTCGCCCTGTATCTCGAACAGGATGACGGTGGTAAAGGCATTCTGAACGCGTTGGATGAGTCTAAAAGGCAGCGTCCGGAACTGGATGGGCGGGTGCTGGTCGACTGGCATCGTGc                                                                                                                                                                                                <  2:291146/244‑1 (MQ=255)                                         aaaaagaaaaaacacagtGCCCAACGAACCAAGATGTATAAATCAGTTGATGATGTAGATGTCTGTTACGATCCCGCCGACTTCCGGGAGACGATGCGGGAAAAAAGAGCCAGAGAGAAGCAGCGCATTTGCATTGTCGTCCTGTATCTCGAACAGGATGACGGTGGCAAAGGAAGTATGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGc                                                                                                                                                                                                <  1:607498/234‑1 (MQ=255)                                                    acacCTGGCCCAACTACCCAAGATTGCTCAATCAGTTGATGATGTCGATTTCTTTGTCGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGAc                                                                                                                                                                                     <  1:98973/251‑1 (MQ=255)                                                             ccAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGCTTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGgcg                                                                                                                                                                            >  1:183745/1‑251 (MQ=255)                                                                    cccAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGa                                                                                                                                                                                                                                    <  1:348701/188‑1 (MQ=255)                                                                    cccAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGa                                                                                                                                                                                                                                    >  2:348701/1‑188 (MQ=255)                                                                        agctttCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGCAGCAGCGCATTTGCATTGTCGCCCTGAATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACGGGCATCGTGCACAACGTGGACGCATGGGCGCTGCGGCTTCt                                                                                                                                                                 >  2:557044/4‑251 (MQ=255)                                                                           tttCTCAATCAGTTGATGATGACGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCGTCTGGAAAAAATAGCCCACGCGCATCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATAACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGCCTAAAAAACCGCGTCCGACACTGGATGAGCGGTAGCTGGTCAGCTGGCATCTTGCACACCGTGGACGCATTGGCTCTGCGGCATCtccc                                                                                                                                                              >  2:377346/1‑248 (MQ=255)                                                                                      gTTGATGATGTCGCTTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGACGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGAAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTACCGc                                                                                                                                                      >  2:347394/1‑248 (MQ=255)                                                                                       ttGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGAc                                                                                                                                                  >  1:420175/1‑251 (MQ=255)                                                                                                 cGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATATCCAGCGCGAAGCAGCGCATTTGCATTGGCGCACTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACGGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTATCGCTGACTGGTACTGcc                                                                                                                                        <  1:602838/251‑1 (MQ=255)                                                                                                         tttACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATg                                                                                                                                   <  1:408535/248‑1 (MQ=255)                                                                                                            aCGCGCCCGCAGACTACCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTGGCATTGTCGCCCTGTAGCTCGAACAGGATGACGGTGGCAAAGGCATTCGGAACGCGTTGTAGGAGACTAAAAGGCAGCGTCCGGAACGGGATGGGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTTGACGCAGTGGCGCTGCGGGATCGAACACTAACGCTGACGGGTACTGCCGCAGGGTGCAg                                                                                                                             <  2:71505/251‑1 (MQ=255)                                                                                                                                 agaCGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGCCGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGAAGATGtt                                                                                                        >  2:613649/1‑251 (MQ=255)                                                                                                                                          tgGAAAAAAAAGCCAGCGCGAAGCAGCGCGTTTGCATTGTCGCCGTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGAGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACCACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAATATCCGGGCGGAGATGTTCCGGTTTAt                                                                                               >  2:353186/1‑251 (MQ=255)                                                                                                                                                 aaTAGCCAGCGCGAAGCAGCGCATTTGCATTGGCGCCCTGTATCTTGAACAGGATGCCGGTGGCAAAGGCATGCTGAACGCGTTGTATGAGTCTAAAAGGCATCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTTCACAACGGTGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTTCTGCCGCATGGCGCAGTAAAATCCGGGCGTAGATTTTCCGGTTTATGGCggtt                                                                                        <  2:430907/251‑4 (MQ=255)                                                                                                                                                             gaatcagcttaTTTGCATTGCCGCCCTGTATCTCGAACTGGATTAAGGTGGTTAACGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCAGGAACTGGATGTGCGGGTGCTGGTCGACGGGCTTCGTGCACAACGTGGTCGCATTGGCGCTTCTGCATCTAACACTAACGCTGACTGTTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATTGCGTTCcaatcattcctc                                                                            <  2:80161/241‑7 (MQ=255)                                                                                                                                                             gAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGCACAGGATGAAGGTGGCACAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTACCACTAACGATGACTGGTACTGCCGCATGGCGCAGGCAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCccatcacaactc                                                                            >  2:278928/1‑245 (MQ=255)                                                                                                                                                              aaGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACATGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCTTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCg                                                                           <  1:484907/251‑1 (MQ=255)                                                                                                                                                                 cagcGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTg                                                                         >  1:422124/1‑250 (MQ=255)                                                                                                                                                                         ttGCATTGTCGCCCTGTATCTCGAACCGGATGACGATGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGAACAACGTGGCCGCATAGGAGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCCGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTGCCAATCCATACTCGTGAGGCCCTTg                                                                >  2:96910/1‑251 (MQ=255)                                                                                                                                                                           gCATGGTCGCCCGGTATCACGAACAGGATGACGGTGGCAAAGGCGTTCTGAACGCGTTGTATGAGACGAAAAGGCAGCGGCCGGATCTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGGCTGGTACTGCCGCATGTCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCTAATCAATACTCGTGAAGCCCTTga                                                               <  2:480533/250‑2 (MQ=255)                                                                                                                                                                                  cGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGAc                                                                                                                                                                                     <  1:667471/125‑1 (MQ=255)                                                                                                                                                                                  cGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGAc                                                                                                                                                                                     >  2:667471/1‑125 (MQ=255)                                                                                                                                                                                            ctcGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGATGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATCCTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAgg                                             >  2:510932/1‑251 (MQ=255)                                                                                                                                                                                                 aCAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCttt                                         >  1:544420/1‑250 (MQ=255)                                                                                                                                                                                                        gACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGCCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                        >  1:201880/1‑228 (MQ=255)                                                                                                                                                                                                        gACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcga                                 <  1:511249/251‑1 (MQ=255)                                                                                                                                                                                                        gACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGGGCGGAAATGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACGCGTGAAGCCCTTGGTGTTCTg                                                        <  2:201880/228‑1 (MQ=255)                                                                                                                                                                                                              ggCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGATGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAg                           <  1:510932/251‑1 (MQ=255)                                                                                                                                                                                                                        ttCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTCTAAAGGCTTTATCATAGACGATAGCGTACTTtat                 >  2:373723/1‑251 (MQ=255)                                                                                                                                                                                                                          ttGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTTGATGTGCGGGTGCTGGTCGACGTGCAGCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATTGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAAGACTCGTGAAGCCCTTTGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTtata                <  2:545798/249‑1 (MQ=255)                                                                                                                                                                                                                                       aTGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTg  <  1:399208/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      GAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTG  >  NC_000913/2722711‑2723191

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