Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,826,697	G→C	5.2%	W46S (TGG→TCG)	pncC →	NMN amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,826,697	0	G	C	5.2%	246.1 / 2.1	76	W46S (TGG→TCG)	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 = 8.88e-01

GGAGCACCCACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTGG  >  NC_002947/1826550‑1826838                                                                                                                                                         |                                                                                                                                              ggAGCACCCACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTgg                                                                                                                                              >  8:302520/1‑149 (MQ=255) ggAGCACCCACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTgg                                                                                                                                              <  2:158886/149‑1 (MQ=255)  gAGCACCCACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGt                                                                                                                                             >  4:39055/1‑149 (MQ=255)       cccACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAg                                                                                                                                        <  4:61601/149‑1 (MQ=255)       cccACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAg                                                                                                                                        <  4:37578/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                                                                                    <  7:278262/149‑1 (MQ=255)                                                           ccGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGc                                                                                    <  6:263270/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                                                                         <  4:285238/142‑1 (MQ=255)                                                                                 gAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCgttg                                                                         >  3:285238/1‑142 (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                      >  3:122237/1‑149 (MQ=255)                                                                                                                              tCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCcgc                      >  1:195000/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   >  7:294672/1‑149 (MQ=255)                                                                                                                                                 cTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTg   <  1:334515/149‑1 (MQ=255)                                                                                                                                                  tCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAGGCCATGGTCCGCGGGGCTCCGGCTGCCGGTgg  >  8:58421/1‑149 (MQ=255)                                                                                                                                                    ggCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATACCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTgg  <  5:385227/147‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                              GGAGCACCCACATGGACCCGATCACTACGCTCGCCACCCGCCTGGGTGAACACCTGCGCCGATTCAA‑TGCGCAGGT‑‑‑GAC‑TACCGCCGAATCCTGCACCGGTGGTGGCATTGCCGAAGCCATCACCCGCATACCTGGCAGCTCGGCCTGGTTCGAGGCCGGTTACGTCACCTATTCCAATGCCCAGAAAACCCGCCAGCTGGGGGTGCCTGAAGCGTTGTTCGGTCAGGTGGGGGCGGTCAGCCAGGAAGTGGTCGAAGCCATGGTCCGCGGAGCTCAGGCTGCCAGTGG  >  NC_002947/1826550‑1826838

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