Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	498,683	G→A	5.2%	R210H (CGT→CAT)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,683	0	G	A	5.2%	241.2 / 2.2	77	R210H (CGT→CAT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

TTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGT  >  NC_000913/498537‑498830                                                                                                                                                   |                                                                                                                                                    ttACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGtt                                                                                                                                                   >  2:374443/1‑149 (MQ=255) ttACGCCGATAACAACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGtt                                                                                                                                                   <  6:265866/149‑1 (MQ=255)    cGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGcc                                                                                                                                                >  7:46515/1‑149 (MQ=255)    cGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACACGCAACTTTGcc                                                                                                                                                >  7:96540/1‑149 (MQ=255)     gCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCg                                                                                                                                               >  2:316870/1‑149 (MQ=255)              cACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCACAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGa                                                                                                                                       >  6:227496/1‑148 (MQ=255)                cGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTc                                                                                                                                    <  1:5568/149‑1 (MQ=255)                      cATTAATGCACTGGTGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAAc                                                                                                                              >  1:355593/1‑149 (MQ=255)                       aTTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCTCAACGTTGCCGCACAACGACTCGTGAACt                                                                                                                             <  1:236237/149‑1 (MQ=255)                       aTTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACt                                                                                                                             >  6:17863/1‑149 (MQ=255)                       aTTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACATTGCCGCACAACGACTCGTGAACt                                                                                                                             >  4:670/1‑149 (MQ=255)                       aTTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACATTGCCGCACAACGACTCGTGAACt                                                                                                                             <  3:670/149‑1 (MQ=255)                              cACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGt                                                                                                                                  <  5:309739/137‑1 (MQ=255)                              cACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGt                                                                                                                                  >  6:309739/1‑137 (MQ=255)                               aCTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCtt                                                                                                                        >  6:102876/1‑146 (MQ=255)                               aCTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCtt                                                                                                                        <  5:102876/146‑1 (MQ=255)                                        cAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGa                                                                                                            <  2:314313/149‑1 (MQ=255)                                         aGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCTCAACGACTCGTGAACTGGCTTCCGCATTGGGGAt                                                                                                           >  7:188123/1‑149 (MQ=255)                                         aGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGAt                                                                                                           <  2:355593/149‑1 (MQ=255)                                          gCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACACGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATg                                                                                                          >  6:59409/1‑149 (MQ=255)                                                gcgcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGAt                                                                                                           >  3:200760/1‑142 (MQ=255)                                                gcgcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGAt                                                                                                           <  4:200760/142‑1 (MQ=255)                                                 cgcTTCTTTTCCCAAACATGGCGAACCGGATCTGCTACTGCTCTATTTTCATGGCATTCGCCAGAGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCTCAACGTCTCGTTAACTGGCTTCCGCATTGGGGATTGCACCgg                                                                                                   <  8:188123/149‑1 (MQ=255)                                                  gcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGca                                                                                                                                             <  6:39717/106‑1 (MQ=255)                                                  gcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGca                                                                                                                                             >  5:39717/1‑106 (MQ=255)                                                  gcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCgg                                                                                                   >  3:195075/1‑148 (MQ=255)                                                  gcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCgg                                                                                                   >  6:332994/1‑148 (MQ=255)                                                      cttTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGaaa                                                                                                <  5:17863/147‑1 (MQ=255)                                                         ttGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGaaaaa                                                                                              <  8:46515/146‑1 (MQ=255)                                                           gCCAAACATGGTGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGtga                                                                                          <  3:295567/148‑1 (MQ=255)                                                           gCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGtga                                                                                          <  1:374443/148‑1 (MQ=255)                                                            ccAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTgg                                                                                                               <  1:328066/126‑1 (MQ=255)                                                            ccAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTgg                                                                                                               >  2:328066/1‑126 (MQ=255)                                                            ccAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGTGATGGCACCGGAAAAAGtgat                                                                                         >  2:267571/1‑148 (MQ=255)                                                            ccAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGtgatg                                                                                        <  7:331230/149‑1 (MQ=255)                                                             cAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGtgatga                                                                                       >  4:309887/1‑149 (MQ=255)                                                                     gCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCttt                                                                               <  5:124792/149‑1 (MQ=255)                                                                       gAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCa                                                                             <  1:316870/149‑1 (MQ=255)                                                                          ccGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCTCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGTGTTTGGCACCGGAAAAAGTGATGATGACCTTTCAGTc                                                                          <  5:353979/149‑1 (MQ=255)                                                                            ggATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTcgc                                                                        >  5:106557/1‑149 (MQ=255)                                                                             gATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTcgcg                                                                       >  4:266765/1‑149 (MQ=255)                                                                                  gCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGACGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGTGGATGGCACCGGAAAAATAGATGATGAGCTTTCATTCGCGCTTTg                                                                  >  1:284631/1‑149 (MQ=255)                                                                                   cTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCtt                                                                    <  2:93895/146‑1 (MQ=255)                                                                                      cTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCg                                                              <  2:94891/149‑1 (MQ=255)                                                                                       tGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCgg                                                             <  5:227496/149‑1 (MQ=255)                                                                                       tGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCgg                                                             >  7:164489/1‑149 (MQ=255)                                                                                       tGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGACGCACAACGACTCGTGAACTGGCTTCCGAATTCGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCgg                                                             >  7:300477/1‑149 (MQ=255)                                                                                         ctctctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                                                            <  5:332994/148‑1 (MQ=255)                                                                                          tctctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCTATGATTTCCCGCATCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTc                                                                          <  3:100614/133‑1 (MQ=255)                                                                                          tctctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACTTTGCCGCACAAAGACTCGTCAACAGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATTACATTTCAGTATCGCTTTGGTCGGGaa                                                          >  8:78745/1‑149 (MQ=255)                                                                                          tctctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGTGGATTGCACCCGAAAAAGTGATGATGACCTTTCAGTc                                                                          >  4:100614/1‑133 (MQ=255)                                                                                            tctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCTCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGaa                                                          <  8:96540/147‑1 (MQ=255)                                                                                            tctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAAcc                                                        >  7:343825/1‑149 (MQ=255)                                                                                            tctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAAcc                                                        >  7:322122/1‑149 (MQ=255)                                                                                              ttATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCGTTCAGTCGCGCTTTGGTCGGGAACCCt                                                      >  4:7674/1‑149 (MQ=255)                                                                                                   aTGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCt                                                  >  5:210978/1‑148 (MQ=255)                                                                                                           ccccAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTgg                                                                 <  5:26593/125‑1 (MQ=255)                                                                                                           ccccAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCGTTgg                                                                 >  6:26593/1‑125 (MQ=255)                                                                                                                         aGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGc                           >  8:39245/1‑149 (MQ=255)                                                                                                                          gATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCt                          <  3:266765/149‑1 (MQ=255)                                                                                                                          gATGAAGGCGATGATTACCCGCAACATTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGc                           >  4:342427/1‑148 (MQ=255)                                                                                                                          gATGAAGGCGATGATTACCCGCAACATTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGc                           <  3:342427/148‑1 (MQ=255)                                                                                                                                    atgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTc                <  7:78745/149‑1 (MQ=255)                                                                                                                                    aggatTACCCGCAACGTTGCCTCACAACGACTCGGGAACTGGCGTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTc                <  5:303006/147‑1 (MQ=255)                                                                                                                                        ttACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCtt                                          >  6:298684/1‑119 (MQ=255)                                                                                                                                        ttACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCtt                                          <  5:298684/119‑1 (MQ=255)                                                                                                                                        ttACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCg               >  5:136986/1‑146 (MQ=255)                                                                                                                                            ccGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGaaa         >  5:244323/1‑148 (MQ=255)                                                                                                                                              gCACCGTTGCCGACCAACGACTCGGGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGCTGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAg      <  5:257506/149‑1 (MQ=255)                                                                                                                                              gCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCTCATTGGGGATGTCACCGTAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAg      <  1:267571/149‑1 (MQ=255)                                                                                                                                              gCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCCACGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGCTGCCTTCTACCGCCGAAACGCTGAAAAGGCTCGGAGAAAAAg      >  6:257506/1‑149 (MQ=255)                                                                                                                                              gCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAg      <  7:39245/149‑1 (MQ=255)                                                                                                                                              gCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAg      >  8:132269/1‑149 (MQ=255)                                                                                                                                                aaCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTg                                                                                                                >  6:17542/1‑41 (MQ=255)                                                                                                                                                aaCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTg                                                                                                                <  5:17542/41‑1 (MQ=255)                                                                                                                                                 aCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCg   <  6:210978/149‑1 (MQ=255)                                                                                                                                                  cGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGt  <  8:164489/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGT  >  NC_000913/498537‑498830

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