Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,720,668	Δ1 bp	coding (716/2661 nt)	pka →	protein lysine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,720,666	0	C	.	98.2%	244.6 / ‑0.3	57	coding (714/2661 nt)	pka	protein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (29/26);  total (31/26)
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

ATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAGCCCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATGCGCC  >  NC_000913/2720518‑2720805                                                                                                                                                     |                                                                                                                                             atgaATTGCTTGACTATCTGGCACGCGACAGTAAACCCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAGc                                                                                                                                              <  7:131398/149‑1 (MQ=255) atgaATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAGc                                                                                                                                              <  4:260465/149‑1 (MQ=255) atgaATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGGGATTAAAAGCGGACGTAGc                                                                                                                                              >  5:293563/1‑149 (MQ=255)     aTTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAGccggc                                                                                                                                          <  7:175386/149‑4 (MQ=255)        gCTTGACTATCTGGCACGCGAGAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGCCGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAACCCGATTCTGGTGATTAAAAGCGGACGTCG‑Ccggcgg                                                                                                                                       >  8:340862/1‑148 (MQ=255)        gCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑Ccggcgg                                                                                                                                       <  6:300919/148‑1 (MQ=255)        gCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑Ccggcgg                                                                                                                                       <  7:103360/148‑1 (MQ=255)         cTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGca                                                                                                                                     >  4:295706/1‑149 (MQ=255)                     gCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCt                                                                                                                         >  7:122144/1‑148 (MQ=255)                     gCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTc                                                                                                                        <  4:177455/149‑1 (MQ=255)                       aCGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCaa                                                                                                                      >  5:120220/1‑149 (MQ=255)                        cgcgACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGTCGGCACAGCGAC‑TGCTCAac                                                                                                                     <  2:335182/149‑1 (MQ=255)                           gACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAacac                                                                                                                   >  5:21550/1‑148 (MQ=255)                           gACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACacg                                                                                                                  >  5:261782/1‑149 (MQ=255)                                tAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACgg                                                                                                              >  7:354743/1‑148 (MQ=255)                                  aaaCCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAg                                                                                                           <  3:349209/149‑1 (MQ=255)                                    aCCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGa                                                                                                         <  3:103402/149‑1 (MQ=255)                                               ccTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑Tg                                                                                                                           <  4:58666/120‑1 (MQ=255)                                               ccTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑Tg                                                                                                                           >  3:58666/1‑120 (MQ=255)                                                   ctctATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCAt                                                                                            <  1:60818/147‑1 (MQ=255)                                                   ctctATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATgg                                                                                          >  5:226366/1‑149 (MQ=255)                                                     ctATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGa                                                                                        >  7:205069/1‑149 (MQ=255)                                                      tATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAGCCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAAt                                                                                                       <  6:325160/133‑1 (MQ=255)                                                      tATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAGCCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAAt                                                                                                       >  5:325160/1‑133 (MQ=255)                                                        tctcGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGc                                                                                     >  3:202008/1‑149 (MQ=255)                                                          tcGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGTTTAAAA‑CGGGCCTTGGCGGGC‑GCCCACCGGCATCAAAAACCGGACGGCAGGAATGGACCCGCCAGGGGATGCgg                                                                                   >  3:100722/1‑149 (MQ=255)                                                          tcGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCgg                                                                                   >  2:288019/1‑149 (MQ=255)                                                            gAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCgg                                                                                   <  8:122144/147‑1 (MQ=255)                                                            gAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCgg                                                                                   <  6:120220/147‑1 (MQ=255)                                                            gAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGc                                                                                  >  1:257443/1‑148 (MQ=255)                                                              aCAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTAt                                                                               <  1:12278/149‑1 (MQ=255)                                                                 gTTAAGCTACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCa                                                                            <  6:21550/149‑1 (MQ=255)                                                                           gcgcgacgcTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGt                                                                  <  4:358163/149‑1 (MQ=255)                                                                             gcgacgcTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGttt                                                                >  1:344765/1‑149 (MQ=255)                                                                             gcgacgcTTTGTTTCGGCGGCCCGTAGGGCCGCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGttt                                                                <  7:218953/149‑1 (MQ=255)                                                                               gacgcTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGTACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGttt                                                                <  5:187052/147‑1 (MQ=255)                                                                               gacgcTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCAGGTGATTAAAAGCGGACGTAG‑CAGGCGGCACAGCGAC‑TGCTCAAAACGACGGCAGGAATGGACACGGCATGGGTTGCGGCTATTCAGCGTGCCGGTtagt                                                              >  3:71801/1‑149 (MQ=255)                                                                                    tttGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCgg                                                         >  2:182225/1‑149 (MQ=255)                                                                                      tGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCgg                                                         <  8:354743/147‑1 (MQ=255)                                                                                      tGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCgg                                                         <  2:344765/147‑1 (MQ=255)                                                                                       gTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTa                                                      >  8:114654/1‑149 (MQ=255)                                                                                                    gTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACg                                         <  8:205069/149‑1 (MQ=255)                                                                                                     tAGTGCCTCGCGTAATAAACAGATTCTGGTGATAAAAAGAGGACGTA‑CCCGCCTGCACAGCGAC‑AGCTCAACACGACGGCAGGAATGGACCCGGCATCGGATGCGCCTATTCAGCGTGCCGCTttgtt                                                             >  8:349466/1‑128 (MQ=255)                                                                                                       gTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGc                                      >  6:158751/1‑149 (MQ=255)                                                                                                         gCCGCGCGTAATAACCCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTg                                    <  4:141708/149‑1 (MQ=255)                                                                                                          ccTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCATGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGt                                   <  6:113017/149‑1 (MQ=255)                                                                                                            tCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACATGACACCCACGAGCTGttt                                 >  6:160621/1‑149 (MQ=255)                                                                                                              gcgTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTc                               >  7:127246/1‑149 (MQ=255)                                                                                                                  aataaACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTcggcg                           >  5:253501/1‑149 (MQ=255)                                                                                                                   ataaACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTcggcg                           <  1:288019/148‑1 (MQ=255)                                                                                                                   ataaACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTcggcg                           <  2:257443/148‑1 (MQ=255)                                                                                                                    taaACCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGt                         <  4:202008/149‑1 (MQ=255)                                                                                                                      aaacgtttCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTg                        <  4:100722/142‑1 (MQ=255)                                                                                                                       aCCGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTgg                       >  3:331063/1‑148 (MQ=255)                                                                                                                        ccGATTCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGa                      >  4:54402/1‑148 (MQ=255)                                                                                                                            ttCTGGTGATTAAAAGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAAccc                 >  3:205403/1‑149 (MQ=255)                                                                                                                                       aaaaGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGc                                                           <  6:357959/96‑1 (MQ=255)                                                                                                                                       aaaaGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGc                                                           >  5:357959/1‑96 (MQ=255)                                                                                                                                        aaaGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATgc     >  2:25550/1‑149 (MQ=255)                                                                                                                                          aGCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATgcgc   <  6:261782/149‑1 (MQ=255)                                                                                                                                           gCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCAc                                          <  6:373839/109‑1 (MQ=255)                                                                                                                                           gCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCAc                                          >  5:373839/1‑109 (MQ=255)                                                                                                                                           gCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATgcgac  >  2:308146/1‑147 (MQ=255)                                                                                                                                           gCGGACGTAG‑CCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATGCGcc  >  7:41220/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                             ATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATTCTGGTGATTAAAAGCGGACGTAGCCCGGCGGCACAGCGAC‑TGCTCAACACGACGGCAGGAATGGACCCGGCATGGGATGCGGCTATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATGCGCC  >  NC_000913/2720518‑2720805

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