Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	863,376	A→T	11.2%	A170A (GCA→GCT)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	863,376	0	A	T	11.2%	127.8 / 6.1	60	A170A (GCA→GCT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (32/19);  new base T (3/4);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCC  >  NC_002947/863230‑863523                                                                                                                                                   |                                                                                                                                                    gCACCGCTTTATCTACAGTTCGTCGACAGCAGGGTGTCCACTGTGTTCGAGCAGGCAAAAGATACGGTCAGCGAGCACAAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTAGATCGAGGCGCTGGcagc                                                                                                                                                   <  2:295520/149‑1 (MQ=255)       cTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAgcgc                                                                                                                                             >  5:66317/1‑149 (MQ=255)            tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGc                                                                                                                                        >  5:150034/1‑149 (MQ=255)             ccACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                       >  7:5532/1‑149 (MQ=255)             ccACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                       >  6:101454/1‑149 (MQ=255)             ccACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                       >  5:223/1‑149 (MQ=255)             ccACAGTGCGCGCACCGCAGGGGGACCACGGTGAGGGAGGAGGCAAAACAGACGGGCTGCGAGCACCCGTTACGGTTGCAGATGCGCGTGCTGCAGCCGTTCTCTGTGCACCCGGCGTAGTTCGCGGCGCTGGCCGCCTGCGCCCGGc                                                                                                                                        <  4:359910/148‑1 (MQ=255)              cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcctcgatg                                                                                                                                      <  6:28648/149‑7 (MQ=255)                  gTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCCGCCAGCGCCCGGCCGTAcc                                                                                                                                  >  1:305913/1‑149 (MQ=255)                   ttCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCt                                                                                                                                 >  5:386868/1‑149 (MQ=255)                     cGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATCAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTg                                                                                                                                >  5:130698/1‑148 (MQ=255)                     cGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCACATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTg                                                                                                                                >  7:113103/1‑148 (MQ=255)                               aCGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTa                                                                                                                     <  2:62077/149‑1 (MQ=255)                                 ggTGACCACTGTGATCGAGCAGGCAAAACATCCGGTCAGCTAGCACCAGTTCCCGTTCCAGATGCGCGTGCTGCAGCCGTTCTATTACCACCCGGCGTACATCGAGCCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACg                                                                                                                   <  6:130698/149‑1 (MQ=255)                                 ggTGACCACTGTGAGCGACCAGTCAAAACAGACGCTCTCCTAGCACCAGTTACCGTTGCTGATGCTCGTGGTGCAGCGTTTCTTTGAGCACCGGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACg                                                                                                                   <  4:75693/149‑1 (MQ=255)                                             gATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcctcg                                                                                                                                         >  7:185736/1‑112 (MQ=255)                                             gATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcctcg                                                                                                                                         <  8:185736/115‑4 (MQ=255)                                                   gCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCt                                                                                                 <  1:347207/149‑1 (MQ=255)                                                       gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCa                                                                                             >  3:24868/1‑149 (MQ=255)                                                        cAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCAt                                                                                            >  4:328918/1‑149 (MQ=255)                                                        aaaaaCAGACGGTCAGAGAGCACTAGTTACGGTTGAAGATGCGCGTGTTGATGCCGTTATATGAGCACCCGGCTTACATGGAGGCGGTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCa                                                                                             <  8:113103/147‑1 (MQ=255)                                                                   gTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCgg                                                                                  >  4:81005/1‑148 (MQ=255)                                                                     cAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAc                                                                               >  7:73260/1‑149 (MQ=255)                                                                        cGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGCGGCGCTGGCAGCCAGCGCCAGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACgg                                                                             >  2:162669/1‑148 (MQ=255)                                                                          atcACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGCGCCCCCGGCGTACAGCGTGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                          <  2:437828/147‑1 (MQ=255)                                                                          accACCAGTTCGCGTTGCTGAGGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCTTACTTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                          <  4:24868/147‑1 (MQ=255)                                                                                 gTTACCGTTGCAGATGCGCGTGCTGCAGACGTTCTATGAGAACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGACCGGCCGTACCTGGAGCAGGGTTACGAACACATGCTACAGAGCTTCCATCGTTTGCAGGAACGGAAACTgaag                                                                    >  5:91283/1‑148 (MQ=255)                                                                                     ccGTTGCTGAGGCGCGTGCTGCAGCCGTTCTATGACCACCCGGCGTACATCGCGGCGCGGCCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAAc                                                                <  3:175514/148‑1 (MQ=255)                                                                                      cGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCCGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTg                                                              >  8:323636/1‑149 (MQ=255)                                                                                       gTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTgg                                                             >  8:140579/1‑148 (MQ=255)                                                                                           cacctGCGCGTGCGGCCGCCGTTCTATCCCCACCCGGCGTACATGGACGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTc                                                          <  2:305913/144‑1 (MQ=255)                                                                                            aGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGtttac                                                         >  3:263471/1‑146 (MQ=255)                                                                                            aGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTcc                                                         >  8:176037/1‑148 (MQ=255)                                                                                              aTGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCg                                                        >  2:201768/1‑147 (MQ=255)                                                                                                  ggggggtgcCGCCGTTCTATGAGCACCGGGCGTACATCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAgggg                                                  <  4:263471/143‑1 (MQ=255)                                                                                                     tgctgcAGCCGTTCTATGACCCCCCGGCGTACATCGCGGCTCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTc                                               <  8:5532/149‑1 (MQ=255)                                                                                                     tgctgaAGCGTTTCTTTGCCCACCCGGCGTACATCCAGCCGCTGTCAGCCACCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTc                                               <  5:101454/149‑1 (MQ=255)                                                                                                      gctgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTc                                               >  3:1672/1‑148 (MQ=255)                                                                                                          cAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCa                                          >  2:300511/1‑149 (MQ=255)                                                                                                            gCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACg                                        >  2:364628/1‑149 (MQ=255)                                                                                                             ccGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGa                                       >  5:216062/1‑149 (MQ=255)                                                                                                              cGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAc                                      <  7:176037/149‑1 (MQ=255)                                                                                                              cGTTCTATCAGCACCCGGCGTACAGCGCGGGGCGGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAc                                      <  4:1672/149‑1 (MQ=255)                                                                                                               gTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAcc                                     >  7:388112/1‑149 (MQ=255)                                                                                                                ttCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                    >  6:362594/1‑149 (MQ=255)                                                                                                                  cTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTg                                   <  3:269787/148‑1 (MQ=255)                                                                                                                  cTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTg                                   >  4:269787/1‑148 (MQ=255)                                                                                                                   tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcg                                 >  2:241110/1‑149 (MQ=255)                                                                                                                    aagAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACACTTTCCGAAGGGCGTCAAGCACGACCTgcgc                                >  2:318239/3‑149 (MQ=255)                                                                                                                     tccgCACCCGGCGTACATCGAGGCGCGGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                <  6:386868/145‑1 (MQ=255)                                                                                                                     tGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                <  6:216062/148‑1 (MQ=255)                                                                                                                     tGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgagcg                               >  3:28605/1‑149 (MQ=255)                                                                                                                       accACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCt                             <  3:81005/147‑1 (MQ=255)                                                                                                                                       ccgggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg             >  7:107457/4‑149 (MQ=255)                                                                                                                                       cGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg             >  5:372707/1‑149 (MQ=255)                                                                                                                                       cGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACAACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg             >  1:313450/1‑149 (MQ=255)                                                                                                                                        cgggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACgg            <  8:107457/147‑1 (MQ=255)                                                                                                                                         aGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt           >  7:244335/1‑149 (MQ=255)                                                                                                                                          ggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCAGCAACGGTg          >  5:242706/1‑149 (MQ=255)                                                                                                                                              cTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAc      <  3:328918/149‑1 (MQ=255)                                                                                                                                                 gCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc   <  8:73260/149‑1 (MQ=255)                                                                                                                                                   agccagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgcc  >  3:23933/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCC  >  NC_002947/863230‑863523

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