Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,458,121	C→A	10.3%	intergenic (+21/+19)	PP_1276 → / ← algA	hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,458,121	0	C	A	10.3%	192.5 / 8.1	66	intergenic (+21/+19)	PP_1276/algA	hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base A (3/4);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCAGC  >  NC_002947/1457980‑1458269                                                                                                                                                   |                                                                                                                                                     gACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAgggg                                                                                                                                               >  7:286542/1‑149 (MQ=255)  aCTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAg                                                                                                                                                  >  5:59349/1‑145 (MQ=255)  aCTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAg                                                                                                                                                  <  6:59349/145‑1 (MQ=255)  aCTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGGCTGACCAGGGGc                                                                                                                                              <  5:135008/149‑1 (MQ=255)      ccACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGt                                                                                                                                          >  4:15847/1‑149 (MQ=255)      ccACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGt                                                                                                                                          >  5:252485/1‑149 (MQ=255)        acaAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGt                                                                                                                                          >  8:161905/1‑147 (MQ=255)                 aGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACt                                                                                                                               >  2:121414/1‑149 (MQ=255)                 aGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACt                                                                                                                               <  3:15847/149‑1 (MQ=255)                    gCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTgcg                                                                                                                            <  4:70399/149‑1 (MQ=255)                            cTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTc                                                                                                                    >  8:154218/1‑149 (MQ=255)                             ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACCATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTc                                                                                                                    >  1:107746/1‑149 (MQ=255)                             ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCt                                                                                                                   >  6:253304/1‑149 (MQ=255)                             ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCt                                                                                                                   >  4:317176/1‑149 (MQ=255)                             ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCt                                                                                                                   >  1:132116/1‑149 (MQ=255)                                gAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGA‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCa                                                                                                                <  5:362156/149‑1 (MQ=255)                                 aaTTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCAc                                                                                                               <  1:278976/149‑1 (MQ=255)                                 aaTTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGACGGTCTTCAc                                                                                                               <  2:8495/149‑1 (MQ=255)                                    ttGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCg                                                                                                             >  8:319633/1‑148 (MQ=255)                                      gAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGaaa                                                                                                          >  7:59442/1‑149 (MQ=255)                                            gaggaTGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAAcgccgc                                                                                                    <  3:317176/149‑1 (MQ=255)                                               gatgatTTGCAGCGACCTTTCAGGGTGGTGGGTGACCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCt                                                                                                 >  4:353902/1‑149 (MQ=255)                                                    tttGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATg                                                                                            <  7:161905/149‑1 (MQ=255)                                                    tttGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATg                                                                                            >  4:55718/1‑149 (MQ=255)                                                    tttGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATg                                                                                            >  7:159189/1‑149 (MQ=255)                                                      tGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgg                                                                                           <  5:32350/148‑1 (MQ=255)                                                        cAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                         <  3:40608/148‑1 (MQ=255)                                                         aGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGt                                                                                       <  8:286542/149‑1 (MQ=255)                                                          gCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTg                                                                                      <  3:55718/149‑1 (MQ=255)                                                           cGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTgg                                                                                     <  2:132116/149‑1 (MQ=255)                                                           cGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTg                                                                                      >  2:134563/1‑148 (MQ=255)                                                           cGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTg                                                                                      >  1:246751/1‑148 (MQ=255)                                                               cTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGt                                                                                 >  7:83048/1‑149 (MQ=255)                                                               cTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGt                                                                                 >  5:164912/1‑149 (MQ=255)                                                                tttCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAg                                                                                                                                                  >  4:212884/1‑83 (MQ=255)                                                                tttCAGGGTGGAGGGTGGACCATTCGA‑CATCAACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCAGGTCAGCCCAGGGGAGAGTCAGCATTTACTGCGCGATCGTCATCACCGAAACGCCGCGCTCGATGGGGGTGGACGTg                                                                                >  2:82941/1‑149 (MQ=255)                                                                            gggtggCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATAcaca                                                                    >  6:158035/1‑149 (MQ=255)                                                                             gggggCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCCTCCCTGCCCCTCCCATGACGAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATac                                                                       <  1:121414/146‑1 (MQ=255)                                                                               tggCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                         >  4:140377/1‑125 (MQ=255)                                                                               tggCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                         <  3:140377/125‑1 (MQ=255)                                                                               tggCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATcc                                                                 <  8:59442/149‑1 (MQ=255)                                                                                    cATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                            <  2:107746/149‑1 (MQ=255)                                                                                    cATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                            <  6:164912/149‑1 (MQ=255)                                                                                    cATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                            <  7:154218/149‑1 (MQ=255)                                                                                      ttCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAgc                                                          >  5:94562/1‑149 (MQ=255)                                                                                         gA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATaca                                                                      <  8:29742/134‑1 (MQ=255)                                                                                         gA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATaca                                                                      >  7:29742/1‑134 (MQ=255)                                                                                            cATCTACTTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGACGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAgagatag                                                     >  7:244834/1‑143 (MQ=255)                                                                                            cATCTACGTCATAGCACCGAGAGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTTCGCGATGGTCTTCACCGAAACGACGCGCTCGATGGGGGTGGACGTGCGCCAATACACATCCTCTaa                                                            <  8:244834/142‑1 (MQ=255)                                                                                               cTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGGGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                                   <  6:252485/148‑1 (MQ=255)                                                                                                   gTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCAt                                              >  2:315981/1‑149 (MQ=255)                                                                                                                   ttGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCt                              >  5:170292/1‑149 (MQ=255)                                                                                                                           gCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACt                      >  7:129395/1‑149 (MQ=255)                                                                                                                               g‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCAc                  >  4:167363/1‑149 (MQ=255)                                                                                                                                  tC‑‑‑GCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGa                        <  6:254612/141‑1 (MQ=255)                                                                                                                                  tC‑‑‑GCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGa                        >  5:254612/1‑141 (MQ=255)                                                                                                                                       gCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCg              <  1:134563/149‑1 (MQ=255)                                                                                                                                       gCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTc               >  7:278586/1‑148 (MQ=255)                                                                                                                                       gCCACTGGTCAGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCt                <  7:170962/147‑1 (MQ=255)                                                                                                                                        ccccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga             <  3:353902/149‑1 (MQ=255)                                                                                                                                        ccccTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga             >  6:41343/1‑149 (MQ=255)                                                                                                                                         cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat            <  6:94562/149‑1 (MQ=255)                                                                                                                                         cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga             <  7:42380/148‑1 (MQ=255)                                                                                                                                         cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga             >  8:42380/1‑148 (MQ=255)                                                                                                                                                 aGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                                   >  5:188216/1‑102 (MQ=255)                                                                                                                                                 aGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                                   <  6:188216/102‑1 (MQ=255)                                                                                                                                                   cccAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCAGc  >  3:229723/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGA‑CATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGG‑CTC‑‑‑GCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTCCAGC  >  NC_002947/1457980‑1458269

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