Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,826,702	G→C	5.3%	E48Q (GAG→CAG)	pncC →	NMN amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,826,702	0	G	C	5.3%	247.0 / 2.7	76	E48Q (GAG→CAG)	pncC	NMN amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (38/34);  new base C (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

CCCACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTGGCGCACGCTT  >  NC_002947/1826556‑1826847                                                                                                                                                        |                                                                                                                                                  cccACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAg                                                                                                                                                 <  4:37578/149‑1 (MQ=255) cccACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAg                                                                                                                                                 <  4:61601/149‑1 (MQ=255)    acaTGGACCCGATCACTACGCTCGCCACCTGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGcc                                                                                                                                              >  8:91215/1‑149 (MQ=255)        ggACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGt                                                                                                                                           >  6:242676/1‑148 (MQ=255)         gACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGtt                                                                                                                                          <  3:39055/148‑1 (MQ=255)          aCCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTAc                                                                                                                                        >  2:334515/1‑149 (MQ=255)          aCCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTCGAACCAGGCCgagc                                                                                                                                          <  7:163231/147‑4 (MQ=255)          aCCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTCGAACCAGGCCgagc                                                                                                                                          >  8:163231/1‑144 (MQ=255)           ccccATCACTACGCTCGCCACCCCCCTGGGTTAACAACTGCGCCGTTTTAATTCCTCTCGTCCATACAAACAGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTCGAACCAGGCCga                                                                                                                                            <  1:307574/149‑2 (MQ=255)              gATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCa                                                                                                                                    >  7:137280/1‑149 (MQ=255)               aTCACTACGCTCGCCACCTGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCAc                                                                                                                                   <  7:91215/149‑1 (MQ=255)                    tACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCGA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTAtt                                                                                                                              >  8:25633/1‑149 (MQ=255)                    tACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTAtt                                                                                                                              >  4:298267/1‑149 (MQ=255)                    tACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTAtt                                                                                                                              >  1:205721/1‑149 (MQ=255)                    tACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACAGAGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTAt                                                                                                                               >  1:311832/1‑148 (MQ=255)                      cGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCGA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTc                                                                                                                             <  7:25633/148‑1 (MQ=255)                      cGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTcc                                                                                                                            <  3:71543/149‑1 (MQ=255)                      cGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCATGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTc                                                                                                                             >  5:263270/1‑148 (MQ=255)                         tCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCa                                                                                                                           >  3:5645/1‑147 (MQ=255)                         tCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTCGAACCAGGCCgag                                                                                                                                           >  8:191341/1‑129 (MQ=255)                           gCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAAt                                                                                                                         <  5:98873/147‑1 (MQ=255)                            ccACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGcc                                                                                                                      <  2:16448/149‑1 (MQ=255)                              aCCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCa                                                                                                                    <  7:302520/149‑1 (MQ=255)                                        gTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGc                                                                                                          <  1:32240/149‑1 (MQ=255)                                         tGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGc                                                                                                          >  2:274633/1‑148 (MQ=255)                                             cacCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCt                                                                                                     >  7:324068/1‑149 (MQ=255)                                                   cgcCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTgggg                                                                                                 <  5:242676/147‑1 (MQ=255)                                                    gcCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTg                                                                                              >  4:317505/1‑149 (MQ=255)                                                     ccGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGc                                                                                             <  6:263270/149‑1 (MQ=255)                                                     ccGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGc                                                                                             <  7:278262/149‑1 (MQ=255)                                                           aatcggcgcAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGaa                                                                                        <  8:84711/144‑1 (MQ=255)                                                            aa‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAg                                                                                       <  3:298267/148‑1 (MQ=255)                                                            aa‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGc                                                                                      >  3:295505/1‑149 (MQ=255)                                                                       t‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGt                                                                            >  2:158792/1‑149 (MQ=255)                                                                           gAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCgttg                                                                                  >  3:285238/1‑142 (MQ=255)                                                                           gAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCgttg                                                                                  <  4:285238/142‑1 (MQ=255)                                                                            aC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTc                                                                           <  4:295505/148‑1 (MQ=255)                                                                               tACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAgg                                                                        >  1:282654/1‑149 (MQ=255)                                                                                 ccgccgAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTg                                                                      >  2:197115/1‑149 (MQ=255)                                                                                 ccgccgAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTg                                                                      >  1:92246/1‑149 (MQ=255)                                                                                    ccgAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTgggg                                                                   <  8:137280/149‑1 (MQ=255)                                                                                       aaTCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCg                                                                <  1:197115/149‑1 (MQ=255)                                                                                        aTCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAAt                                                                                                                         >  8:114016/1‑91 (MQ=255)                                                                                        aTCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAAt                                                                                                                         <  7:114016/91‑1 (MQ=255)                                                                                            tGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGACTGAAGCGTTGTTCGGTCAGGTGGTGGCGGTCAg                                                           >  8:200453/1‑149 (MQ=255)                                                                                             gCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGc                                                          >  8:39109/1‑149 (MQ=255)                                                                                               aCCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGc                                                          >  2:196841/1‑147 (MQ=255)                                                                                               aCCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCa                                                        <  2:299746/149‑1 (MQ=255)                                                                                                ccGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAg                                                       <  2:205721/149‑1 (MQ=255)                                                                                                       tggCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTgg                                                <  1:158792/149‑1 (MQ=255)                                                                                                         gCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTc                                              >  5:230874/1‑149 (MQ=255)                                                                                                           aTTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCg                                             >  6:116513/1‑148 (MQ=255)                                                                                                           aTTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCg                                             >  3:44643/1‑148 (MQ=255)                                                                                                            ttGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGaa                                           <  2:282654/149‑1 (MQ=255)                                                                                                             tGCCGAATGCATCATCCGCATCCCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCTCCTATTCCAATGCCCAGATAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCTGTCTGGTGGGGGGGGTCAGCCAGGTAGTGGTCGaa                                           <  8:324068/148‑1 (MQ=255)                                                                                                               ccGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGc                                         >  2:254134/1‑148 (MQ=255)                                                                                                               ccGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGGGGGGGCgg                                                               >  7:235789/1‑126 (MQ=255)                                                                                                                       aTCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCcg                                >  5:102308/1‑149 (MQ=255)                                                                                                                        tCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCcgc                               >  1:195000/1‑149 (MQ=255)                                                                                                                        tCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCcgc                               >  3:122237/1‑149 (MQ=255)                                                                                                                           cccGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTc                                  >  7:318484/1‑143 (MQ=255)                                                                                                                           cccGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTc                                  <  8:318484/143‑1 (MQ=255)                                                                                                                           cccGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATACCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCcgcg                              >  6:385227/1‑147 (MQ=255)                                                                                                                                aTACGTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCa                       <  4:5645/149‑1 (MQ=255)                                                                                                                                 tACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAg                      <  3:317505/149‑1 (MQ=255)                                                                                                                                  aCCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAgg                     >  1:80103/1‑149 (MQ=255)                                                                                                                                   ccTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAgg                     <  6:230874/148‑1 (MQ=255)                                                                                                                                       gCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGc                 <  7:39109/148‑1 (MQ=255)                                                                                                                                       gCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGc                 <  1:254134/148‑1 (MQ=255)                                                                                                                                        cAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCa               <  4:122237/149‑1 (MQ=255)                                                                                                                                           cTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTg            >  1:20971/1‑149 (MQ=255)                                                                                                                                           cTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTg            <  1:334515/149‑1 (MQ=255)                                                                                                                                           cTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTg            >  7:294672/1‑149 (MQ=255)                                                                                                                                            tCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAGGCCATGGTCCGCGGGGCTCCGGCTGCCGGTgg           >  8:58421/1‑149 (MQ=255)                                                                                                                                              ggCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATACCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTgg           <  5:385227/147‑1 (MQ=255)                                                                                                                                                     gtCGAGGCCGGGTACGTCCCCTCTTCCAATGCCCAGAAAACCCTCCAGCTGGGGGTGCCTGAAGCGGTGGTCGGTCAGGTGGGGTCGCTCAGCCAGCAAGTGGTCGCAGCCATGGTCCGCGCCGCTCCGGCTGCCAGTGGCGCCCGcgt  >  8:241034/2‑147 (MQ=255)                                                                                                                                                        |                                                                                                                                                  CCCACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTGGCGCACGCTT  >  NC_002947/1826556‑1826847

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