Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,990,715	A→G	6.2%	W18R (TGG→CGG)	hemC ←	hydroxymethylbilane synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,990,715	0	A	G	6.2%	203.1 / 3.4	64	W18R (TGG→CGG)	hemC	hydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (2/2);  total (31/33)
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.67e-01

CAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACGCCTGCGCAAAC  >  NC_000913/3990569‑3990858                                                                                                                                                   |                                                                                                                                                cAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCag                                                                                                                                                <  7:137924/148‑1 (MQ=255) cAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCag                                                                                                                                                <  7:129252/148‑1 (MQ=255)  agatctTTTACAAATAAGACTTTTCCGCCTACTTTCGCCAGCGGCGTATGAAGAATCACATCGCCGCGCGTCAGCATCGGTACCAGTTCTACGACCAGGCCCGGATTTCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCaga                                                                                                                                               <  7:130378/145‑1 (MQ=255)     tctTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCTGTAACAGTTCAACGACCAGTCCCGGATGGCTCTCAATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTg                                                                                                                                            >  2:260285/1‑148 (MQ=255)       ttttACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTAACAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTg                                                                                                                                            <  8:149222/146‑1 (MQ=255)          tACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAg                                                                                                                                        <  6:209532/147‑1 (MQ=255)            cAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGCCCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTgg                                                                                                                                     >  6:45681/1‑148 (MQ=255)                      ttttCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGCCCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCgt                                                                                                                           >  1:64566/1‑148 (MQ=255)                       tttCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGCCCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCgtg                                                                                                                          >  3:281596/1‑148 (MQ=255)                         tCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCgtgt                                                                                                                         <  6:236975/147‑1 (MQ=255)                                    tCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCtt                                                                                                               <  3:273744/146‑1 (MQ=255)                                    tCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTa                                                                                                              <  8:328041/147‑1 (MQ=255)                                     cGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTaa                                                                                                             <  4:260270/147‑1 (MQ=255)                                                tATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTaa                                                                                                 <  6:172843/148‑1 (MQ=255)                                                   cAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTg                                                                                                                                      <  4:292011/108‑1 (MQ=255)                                                   cAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTg                                                                                                                                      >  3:292011/1‑108 (MQ=255)                                                    aaGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTaa                                                                                                 <  2:255299/144‑1 (MQ=255)                                                    aaGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTaa                                                                                                 >  1:255299/1‑144 (MQ=255)                                                     aGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGc                                                                                            >  6:140252/1‑148 (MQ=255)                                                      gAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCGGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCt                                                                                           >  3:15461/1‑148 (MQ=255)                                                      gAATCACATCGCCGCGCGTCACCATCCGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCGGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCt                                                                                           >  3:253012/1‑148 (MQ=255)                                                         tCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGt                                                                                        >  7:165102/1‑148 (MQ=255)                                                           acaTCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGGCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTa                                                                                      <  8:182327/148‑1 (MQ=255)                                                             aTCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAAGATGCTTGTTAcc                                                                                    >  1:278331/1‑148 (MQ=255)                                                               cgccgcGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTg                                                                                                      >  8:246472/1‑128 (MQ=255)                                                               cgccgcGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTg                                                                                                      <  7:246472/128‑1 (MQ=255)                                                                  cgagcgTCACCATAGGTACCAGTTCAACGACCATGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTc                                                                                 <  3:65705/143‑1 (MQ=255)                                                                        tCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTatcat                                                                         <  4:261675/148‑1 (MQ=255)                                                                           ccATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCg                                                                      >  1:313160/1‑148 (MQ=255)                                                                            cATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGt                                                                     <  7:117871/148‑1 (MQ=255)                                                                               cGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGc                                                                                            <  2:296468/122‑1 (MQ=255)                                                                               cGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGc                                                                                            >  1:296468/1‑122 (MQ=255)                                                                                  tACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCa                                                               <  3:231417/148‑1 (MQ=255)                                                                                  tACCAGTTCAAAGACCAGGCCCGGATGTCTCGCCATCAACTTGTCTTTGACATTGTGTGCATTCAATAGTGCAAGTGGGCTTTTGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTTGTCCa                                                               >  4:165073/1‑148 (MQ=255)                                                                                    ccAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCAt                                                              <  4:53191/147‑1 (MQ=255)                                                                                        ttCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCAt                                                                               <  4:238315/126‑1 (MQ=255)                                                                                        ttCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCAt                                                                               >  3:238315/1‑126 (MQ=255)                                                                                        ttCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGGCTAACATGCTTGTTACCGGCATTATCATCCGTGGGCCATCCTaa                                                         >  6:22955/1‑148 (MQ=255)                                                                                         tCAAAGACCAGGCACGGATGGCTCGCCATCAAATTGTCTTTGACATAGTGTGCATGCCAGAGTGAAAGTGGGCTTTGGCGTGTTGCAATTCTTAAAACAATTTCTAACATGCTTGTTACCGTAAATATCATCAGTGGTCCATCCTac                                                         >  4:233141/1‑146 (MQ=255)                                                                                             cGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCAt                                                              <  7:208438/138‑1 (MQ=255)                                                                                             cGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCAt                                                              >  8:208438/1‑138 (MQ=255)                                                                                              gACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCt                                                   >  1:248265/1‑148 (MQ=255)                                                                                              gACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCt                                                   <  5:140252/148‑1 (MQ=255)                                                                                                    gCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGGGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCTTCCTAACATCCTTGCCag                                             <  4:281596/148‑1 (MQ=255)                                                                                                     cccGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCa                                              >  7:51130/1‑146 (MQ=255)                                                                                                       cGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGt                                                                     <  3:97002/121‑1 (MQ=255)                                                                                                       cGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGt                                                                     >  4:97002/1‑121 (MQ=255)                                                                                                          aTGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGAt                                       >  1:107982/1‑148 (MQ=255)                                                                                                           tGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTTGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATg                                      <  1:260285/148‑1 (MQ=255)                                                                                                             gCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGt                                     <  2:278331/147‑1 (MQ=255)                                                                                                                 gCCATCAACTTGTCTTTGACATAGTGTGCCTGGCGGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATAATCCGTTGTCCATCCTAACATCCTTGCCAGAGTGATGTCAg                                  <  4:253012/146‑1 (MQ=255)                                                                                                                 gCCATCAACTTGTCTTTGACATAGTGTGCCTGCCGGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAg                                  <  4:15461/146‑1 (MQ=255)                                                                                                                  ccATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTTGCAATTCTTAAAACATTGTCTAACATGCTTGTTACAGTCATTATAATCCGTGGTCCATACTAACATCCTTGCCAGAGTGAGTTCAgtgt                               >  4:157065/1‑148 (MQ=255)                                                                                                                  ccATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCAt                                                              <  6:265776/117‑1 (MQ=255)                                                                                                                  ccATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCAt                                                              >  5:265776/1‑117 (MQ=255)                                                                                                                    aTCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGtgttg                             <  5:45681/148‑1 (MQ=255)                                                                                                                     tCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGtgttgt                            >  4:220522/1‑148 (MQ=255)                                                                                                                      cAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTtgtg                           >  1:106481/1‑148 (MQ=255)                                                                                                                        aCTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCAAGAGTGATGTCAGTGTTgtggt                         >  5:528/1‑148 (MQ=255)                                                                                                                            gTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGaa                      <  2:313160/147‑1 (MQ=255)                                                                                                                                 ttaCATAGTGTGCCTTCCAGAGTTCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTa                 <  3:157065/145‑1 (MQ=255)                                                                                                                                 tGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCGGTGTTGTGGTGAAACGTa                 <  2:107982/147‑1 (MQ=255)                                                                                                                                  gACATAGTTTGCCTGCCAGAGTGCATGTGGGCTTTGTCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCTTCTTTATCATCCGTGGTCCATCCTAACATCCTTGCCTTAGTGATGTCAGTGTTGTGGTGAAACGTAg                <  3:233141/147‑1 (MQ=255)                                                                                                                                  gACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAg                <  2:64566/147‑1 (MQ=255)                                                                                                                                    cATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACg             >  8:79753/1‑148 (MQ=255)                                                                                                                                           tgcctgccAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACGCCTgcg       >  5:287521/1‑147 (MQ=255)                                                                                                                                               tgccAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACGCCTGCGCAAAc  >  5:112146/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                CAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACGCCTGCGCAAAC  >  NC_000913/3990569‑3990858

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