Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,486,216	A→C	5.1%	S30R (AGC→CGC)	PP_3978 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,486,216	0	A	C	5.1%	313.6 / 5.6	98	S30R (AGC→CGC)	PP_3978	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/48);  new base C (2/3);  total (47/51)
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.06e-01

CAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAA  >  NC_002947/4486071‑4486362                                                                                                                                                   |                                                                                                                                                   cAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCt                                                                                                                                                 >  8:359159/1‑149 (MQ=255)   aGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGc                                                                                                                                               <  5:131213/149‑1 (MQ=255)     aTCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCTCGCTGccg                                                                                                                                             >  8:111015/1‑148 (MQ=255)      tCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCTCGCTGccgg                                                                                                                                            <  1:105204/149‑3 (MQ=255)       cAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCAt                                                                                                                                            <  2:397056/148‑1 (MQ=255)         aTGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATg                                                                                                                                           >  3:85109/1‑147 (MQ=255)          tgctgcGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCTCGCTGccgggtca                                                                                                                                        <  7:111015/149‑7 (MQ=255)            ctgcGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAgc                                                                                                                                       <  3:281021/148‑1 (MQ=255)             tgcGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCAt                                                                                                                                            <  2:236086/142‑1 (MQ=255)             tgcGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCAt                                                                                                                                            >  1:236086/1‑142 (MQ=255)               cgcCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAgccgcc                                                                                                                                   >  2:41843/1‑149 (MQ=255)                gcCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGccc                                                                                                                                  >  3:275671/1‑149 (MQ=255)                 ccACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCt                                                                                                                                 <  7:359159/149‑1 (MQ=255)                     tACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCTCGCTGccgggtca                                                                                                                                        <  1:365804/138‑7 (MQ=255)                     tACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCTCGCTGccgggtca                                                                                                                                        >  2:365804/1‑132 (MQ=255)                     tACAATTACATCGTTGACGACCCGATCATTACATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGATATGAGCTACACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATg                                                                                                                             >  2:317965/1‑149 (MQ=255)                                  ttGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCTGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCAcc                                                                                                                <  8:353780/149‑1 (MQ=255)                                       gacCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGc                                                                                                           >  1:258266/1‑149 (MQ=255)                                       gacCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGc                                                                                                           >  6:18450/1‑149 (MQ=255)                                        acCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCt                                                                                                          <  8:354380/149‑1 (MQ=255)                                         cccGATCATACCATAGGTATGTGCCACCCGTGGTCTAGCCACAAACCCAACATCGATCGCATCAGCTATGACCTCCACAGCCACACAAAGCATGACCCGGC‑AGCCAGCTGCCAAGAAGCCGCCCTGAACCAGCAtt                                                                                                                      >  3:95543/1‑136 (MQ=255)                                          ccGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCttt                                                                                                        <  7:175316/149‑1 (MQ=255)                                            gATCAATCCATAGGTATGTTCCTCCCGTTGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATTAGCTCCACAGTCACTCAAAGGATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCGTGATGGAGCATTTGCACCGTGGCTTTaa                                                                                                      <  7:329909/149‑1 (MQ=255)                                             aTCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAg                                                                                                     <  3:82463/149‑1 (MQ=255)                                              tCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTg                                                                                                                    <  8:353005/133‑1 (MQ=255)                                                aTTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACg                                                                                                  <  8:62572/149‑1 (MQ=255)                                                  tCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTgg                                                                                                            >  6:91849/1‑137 (MQ=255)                                                  tCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTgg                                                                                                            <  5:91849/137‑1 (MQ=255)                                                    cATAGGTATGTGGCTCCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATTAGCTCCACAGGCACACAAAGCATGACCCGGCAAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACggg                                                                                                <  1:317965/148‑1 (MQ=255)                                                    cATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGAt                                                                                              >  2:381700/1‑149 (MQ=255)                                                    cATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGAt                                                                                              >  3:157572/1‑149 (MQ=255)                                                     aTAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATc                                                                                             >  1:303313/1‑149 (MQ=255)                                                         gTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACCGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATa                                                                                         >  6:227998/1‑149 (MQ=255)                                                         gTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATa                                                                                         >  6:274016/1‑149 (MQ=255)                                                           aTGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATAcc                                                                                       >  5:19108/1‑149 (MQ=255)                                                              tgCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGAc                                                                                                   >  7:353261/1‑134 (MQ=255)                                                              tgCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGAc                                                                                                   <  8:353261/134‑1 (MQ=255)                                                                ccACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACa                                                                                   <  4:85109/148‑1 (MQ=255)                                                                 cACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGt                                                                                 <  2:258266/149‑1 (MQ=255)                                                                 cACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGt                                                                                 <  5:182362/149‑1 (MQ=255)                                                                  aCCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTg                                                                                <  2:285996/149‑1 (MQ=255)                                                                     cGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTa                                                                                                       >  1:402864/1‑123 (MQ=255)                                                                     cGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTa                                                                                                       <  2:402864/123‑1 (MQ=255)                                                                      gtggtgTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCt                                                                            <  1:41843/149‑1 (MQ=255)                                                                            tAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCttt                                                                                                        <  7:157869/115‑1 (MQ=255)                                                                            tAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCttt                                                                                                        >  8:157869/1‑115 (MQ=255)                                                                              gCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCg                                                                    <  7:118770/149‑1 (MQ=255)                                                                               ccACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGt                                                                   <  6:19108/149‑1 (MQ=255)                                                                                cacaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGt                                                                   >  4:22893/1‑148 (MQ=255)                                                                                cacaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGt                                                                   >  6:147955/1‑148 (MQ=255)                                                                                 acaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGaa                                                                                                                                         >  5:1541/1‑77 (MQ=255)                                                                                 acaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGaa                                                                                                                                         <  6:1541/77‑1 (MQ=255)                                                                                 acaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGt                                                                   >  4:145467/1‑147 (MQ=255)                                                                                 acaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGt                                                                   >  7:390471/1‑147 (MQ=255)                                                                                 acaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTaa                                                                 <  8:390471/149‑1 (MQ=255)                                                                                 acaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTaa                                                                 <  3:22893/149‑1 (MQ=255)                                                                                 acaAACACAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGt                                                                   >  8:184499/1‑147 (MQ=255)                                                                                          aCATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAg                                                                                                                          <  2:156937/83‑1 (MQ=255)                                                                                          aCATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAg                                                                                                                          >  1:156937/1‑83 (MQ=255)                                                                                          aCATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGc                                                        <  2:303313/149‑1 (MQ=255)                                                                                             tCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                                     <  5:18450/149‑1 (MQ=255)                                                                                             tCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                                     >  8:158353/1‑149 (MQ=255)                                                                                              cGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGcc                                                                                                                                              <  1:47818/59‑1 (MQ=255)                                                                                              cGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGcc                                                                                                                                              >  2:47818/1‑59 (MQ=255)                                                                                              cGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAAc                                                    >  3:225099/1‑149 (MQ=255)                                                                                               gATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAAc                                                    <  1:381700/148‑1 (MQ=255)                                                                                                 tGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTggg                                                                              <  8:136875/120‑1 (MQ=255)                                                                                                   gcagcaGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCAcc                                               >  2:326801/1‑149 (MQ=255)                                                                                                     agcagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCtt                                             >  6:135038/1‑149 (MQ=255)                                                                                                     agcagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCtt                                             >  2:145870/1‑149 (MQ=255)                                                                                                        aacTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCg                                          <  1:145870/147‑1 (MQ=255)                                                                                                          cTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGaa                                        <  7:184499/149‑1 (MQ=255)                                                                                                          cTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGaa                                        >  5:100263/1‑149 (MQ=255)                                                                                                           tATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAAc                                       >  2:63064/1‑149 (MQ=255)                                                                                                            aTGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACt                                      >  8:22625/1‑149 (MQ=255)                                                                                                             tGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTTGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTc                                     <  4:275671/149‑1 (MQ=255)                                                                                                                 cTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAg                                  >  6:340005/1‑148 (MQ=255)                                                                                                                    cacaGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCAt                               >  8:6849/1‑148 (MQ=255)                                                                                                                         gCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAggg                         <  5:274016/149‑1 (MQ=255)                                                                                                                             acaAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTa                             <  8:83391/141‑1 (MQ=255)                                                                                                                             acaAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTa                             >  7:83391/1‑141 (MQ=255)                                                                                                                              caAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATg                    <  6:7148/149‑1 (MQ=255)                                                                                                                               aaaGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTTCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATgg                   >  1:307488/1‑149 (MQ=255)                                                                                                                               aaaGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATgg                   >  4:128161/1‑149 (MQ=255)                                                                                                                                aaGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGgt                  <  7:22625/149‑1 (MQ=255)                                                                                                                                    aTGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCTCCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTcc              <  8:53965/149‑1 (MQ=255)                                                                                                                                    aTGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCGCCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTcc              >  7:53965/1‑149 (MQ=255)                                                                                                                                     tGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTcc                                    >  7:197850/1‑126 (MQ=255)                                                                                                                                     tGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTcc                                    <  8:197850/126‑1 (MQ=255)                                                                                                                                      gACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCAtt                              <  4:159206/131‑1 (MQ=255)                                                                                                                                      gACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCAtt                              >  3:159206/1‑131 (MQ=255)                                                                                                                                       aCCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGc                        >  8:234073/1‑136 (MQ=255)                                                                                                                                       aCCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGc                        <  7:234073/136‑1 (MQ=255)                                                                                                                                        cccGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGt          <  5:147955/149‑1 (MQ=255)                                                                                                                                           ggC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACGGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTAGATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTAtt       <  5:227998/149‑1 (MQ=255)                                                                                                                                               aGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTAtttt     >  6:165509/1‑148 (MQ=255)                                                                                                                                                 cGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTaa  <  5:135038/149‑1 (MQ=255)                                                                                                                                                 cGAGCTGCAATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTaa  <  3:128161/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGC‑AGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAA  >  NC_002947/4486071‑4486362

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