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	87.8%	91.0 / 4.9	41	A74T (GCT→ACT)	pssA	phosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (14/22);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

TGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTG  >  NC_000913/2722703‑2723191                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      tGAACAGAGAAGAAATGAACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCAGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGCAACTCGCACAGGATGACGGTGACCACGACATTCTGAACGAGTTCTATGAGACTaaa                                                                                                                                                                                                                                                  >  2:105217/1‑249 (MQ=255)   aaCAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGc                                                                                                                                                                                                                                              >  2:28855/1‑251 (MQ=255)    aCAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGAAATAAACATCAACAACGCCTTGCCCAACTACCCAAGATTTCCCAATTAGTTTAGGATGTCGAGGTCATTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCATCTCATTTGCATTGGCGTCCTTTATCGCGAACAGGATGACGTTGGCAAAGGCATTCTTAACGTGTGGTATGAGACTAAAAGGCa                                                                                                                                                                                                                                             <  2:213896/251‑1 (MQ=255)     cAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGCGGCTAAAAGGccg                                                                                                                                                                                                                                            >  2:215263/1‑249 (MQ=255)                     tgtgATGTTGTAAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGATGATGATGTCGATTTCTTTTACGCTACCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAAAGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGa                                                                                                                                                                                                                              >  1:390021/1‑249 (MQ=255)                       tgATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGTtgtg                                                                                                                                                                                                                          >  1:335867/1‑251 (MQ=255)                              gtCAAAAGGTAATCGTAATAAACATCAACAACACCTTGCCCAACGACCCAAGATTTCTCAAACAGTTGATGATGTCGATTTCTTTTACGCACCCGCCGACGTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGGATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTg                                                                                                                                                                                                                    <  1:321054/250‑1 (MQ=255)                                                  aaCATCAACAACACCTTGCCAAACTACCCAATATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGca                                                                                                                                                                                               <  2:430917/251‑1 (MQ=255)                                                   aCATCAACAACCCCTTGCCCAACTACCCAAGATTTCTCAATCATTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATGCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGcac                                                                                                                                                                                              <  2:52189/251‑1 (MQ=255)                                                    cATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGcaca                                                                                                                                                                                             <  1:450776/251‑1 (MQ=255)                                                            cccccTTTCCCAACTCCCCAAGATTTCTCAATCAGTTGAGTATGTCGATTCCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGAc                                                                                                                                                                                     <  2:55169/248‑1 (MQ=255)                                                                   gCCCAACGACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTg                                                                                                                                                                               <  1:215263/250‑1 (MQ=255)                                                                                                     atgTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTAc                                                                                                                                            <  1:53850/251‑1 (MQ=255)                                                                                                     atgTCGATTTCTTTTACGCTCCCGCCGACTTCCAGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTAc                                                                                                                                            >  1:496767/1‑251 (MQ=255)                                                                                                       gTCGATTGCTTTTACTCCCCCGACGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTg                                                                                                                                          <  2:289554/251‑1 (MQ=255)                                                                                                           aTTTCGTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGc                                                                                                                                      <  1:266733/251‑1 (MQ=255)                                                                                                                        tCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGaa                                                                                                                          >  1:88444/1‑250 (MQ=255)                                                                                                                           cgccgACGTCCGGGAGACGCGGCTGGAAAAAATAGTCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGACCAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGGCCGGAACTGGATGTGCGGGTGCTGGTCGACTTGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGGACTGCTGCATGGCGCAGGAAAATc                                                                                                                      <  2:74651/251‑1 (MQ=255)                                                                                                                             ccgACCTCCGGGAGACTCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATGTGCCTTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTAGGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCg                                                                                                                    <  2:163650/251‑1 (MQ=255)                                                                                                                                         agaCGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGtt                                                                                                        <  1:332998/251‑1 (MQ=255)                                                                                                                                            cGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATGTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCg                                                                                                     <  1:369746/251‑1 (MQ=255)                                                                                                                                              ctgctgGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGt                                                                                                   >  1:179485/1‑251 (MQ=255)                                                                                                                                                gctgGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCTCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACTCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACTCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCAGGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGtt                                                                                                  <  1:161497/250‑1 (MQ=255)                                                                                                                                                         aaTAGCCAGCGCGAAGCAGCGCATTTGAATTGTCGCCCGGGATCTCGAACAGGATGACGGTGGAAAAGGCATTCTGAACGCGGTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTAGGGCGTTc                                                                                        <  2:284322/251‑1 (MQ=255)                                                                                                                                                              ccAGCGCGAATCAGCGCATTTGCATTGGCGCCCTGTATCTCGAACAGGAGGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCAGCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatc                                                                                   <  2:127127/251‑1 (MQ=255)                                                                                                                                                                  cgcgAAGCAGCGCCTTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGGGCACAAAGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTAGGGCGTTCcaatcact                                                                                >  1:145898/1‑248 (MQ=255)                                                                                                                                                                    cgAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGCTGTTCCGGTTTTTGGCGTTCCAATCAATAcc                                                                             >  2:426184/1‑250 (MQ=255)                                                                                                                                                                      aaGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTACCACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCg                                                                           >  2:139303/1‑251 (MQ=255)                                                                                                                                                                          agcGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGaa                                                                       <  1:491677/251‑1 (MQ=255)                                                                                                                                                                                tttGCATTGTCGCCCTGAATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAACGTCCGGCACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGGACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGCTCCAATCACTACTCGTGAAGCCCtt                                                                 >  1:58506/1‑251 (MQ=255)                                                                                                                                                                                                      cGAACAGGATGACGGTGGCAAAGGCATTCTGAACGTGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGc                                            <  2:264663/250‑1 (MQ=255)                                                                                                                                                                                                      cGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTaaa                                               <  1:206339/247‑1 (MQ=255)                                                                                                                                                                                                        aaCAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCttt                                         <  1:457360/251‑1 (MQ=255)                                                                                                                                                                                                           aGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTatc                                      <  2:179485/251‑1 (MQ=255)                                                                                                                                                                                                                          aaGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGt                        >  1:255703/1‑250 (MQ=255)                                                                                                                                                                                                                          aaGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTa                       <  1:478676/251‑1 (MQ=255)                                                                                                                                                                                                                          aaGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTa                       <  1:435555/251‑1 (MQ=255)                                                                                                                                                                                                                           agtcATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGCCCGGAACTGGATGTGCGGGTGCTGGGCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACGAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATTTTCCGGTTTATGGCGTTCCAATTAATACGCGTGAATCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCTTCGACGATAGCGTAc                      <  2:289351/248‑1 (MQ=255)                                                                                                                                                                                                                                     aaCGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCgg            <  1:28855/251‑1 (MQ=255)                                                                                                                                                                                                                                           ttGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGcc        >  1:86055/1‑249 (MQ=255)                                                                                                                                                                                                                                                tGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTg  >  2:4013/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTG  >  NC_000913/2722703‑2723191

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