Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,052,862:1	+A	5.3%	coding (505/1428 nt)	PP_0910 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,052,862	1	.	A	5.3%	308.6 / 9.7	75	coding (505/1428 nt)	PP_0910	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (35/36);  new base A (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CCAGGCCGATCATTGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGCGTACCAGCATAC  >  NC_002947/1052715‑1053008                                                                                                                                                     |                                                                                                                                                    ccaggccGATCATTGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCt                                                                                                                                                      >  4:158123/1‑148 (MQ=255) ccaggccGATCATTGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCt                                                                                                                                                      >  4:81465/1‑148 (MQ=255)  caggccGATCATTGGTGGTGCCACCCGCTCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGTGTGATTTGCATGTCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGTGTTTCGCGGGTGTAGCT‑Ca                                                                                                                                                   <  3:51265/149‑1 (MQ=255)       cGATCATTGGTGGTGCCAGCGGCGGAGAAAAACGCGTGTCGGGGTGCATCAGTCGATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTG‑AGCTACAC‑CCGc                                                                                                                                             <  4:371629/149‑1 (MQ=255)        gATCATTGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCa                                                                                                                                            <  1:160800/149‑1 (MQ=255)            aTTGGTGGTGCCACCGGCGCAGATAAAGGGGTGGCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CTC‑GCGCACCtt                                                                                                                                        <  2:288969/149‑1 (MQ=255)            aTTGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCtt                                                                                                                                        <  8:342550/149‑1 (MQ=255)            aTTGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCtt                                                                                                                                        <  5:240261/149‑1 (MQ=255)             ttGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTg                                                                                                                                       >  6:192977/1‑149 (MQ=255)                gtggtgCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCt                                                                                                                                    >  4:85771/1‑149 (MQ=255)                gtggtgCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCt                                                                                                                                    >  6:250445/1‑149 (MQ=255)                     gCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGCTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATCACAGGGTTTCGCGGTTGTAGCT‑CAC‑CCGCACCTTGCATTCggg                                                                                                                               >  6:382006/1‑149 (MQ=255)                      ccACCCGCGCAGAAAACCGCGGGTCGGGGTGCATCACTCCATCGCCGCGCTGGGTGATTTGCCGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGCAGGACAGGGGTTCGCGGGGGTAGCT‑CAC‑CCGCACCTTGCCTTCggg                                                                                                                               <  3:208476/148‑1 (MQ=255)                        aCCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGgcgc                                                                                                                            >  3:100655/1‑149 (MQ=255)                        aCCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGgcgc                                                                                                                            >  4:231080/1‑149 (MQ=255)                        aCCCGCGCAGAAAAACGCGTGTCGGGGGGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGGTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGgcgc                                                                                                                            <  5:315810/149‑1 (MQ=255)                         cccGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGcc                                                                                                                           <  3:158123/149‑1 (MQ=255)                                aGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCgg                                                                                                                    <  3:113060/149‑1 (MQ=255)                                    aaaCGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTAt                                                                                                                >  8:70323/1‑149 (MQ=255)                                       cgcgTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGc                                                                                                             <  5:192977/149‑1 (MQ=255)                                           tgtCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATg                                                                                                         <  8:217466/149‑1 (MQ=255)                                             tCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTc                                                                                                       <  3:133265/149‑1 (MQ=255)                                                ggggGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CACACCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCt                                                                                                     <  1:264878/149‑1 (MQ=255)                                                      aTCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTg                                                                                              >  4:371667/1‑149 (MQ=255)                                                      aTCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTg                                                                                              >  3:140370/1‑149 (MQ=255)                                                      aTCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTg                                                                                              >  7:319104/1‑149 (MQ=255)                                                      aTCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTg                                                                                              >  3:134715/1‑149 (MQ=255)                                                        cACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTggg                                                                                            >  8:33593/1‑149 (MQ=255)                                                           tCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTAc                                                                                         >  3:61266/1‑149 (MQ=255)                                                           tCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGAGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTAc                                                                                         >  3:43817/1‑149 (MQ=255)                                                             cATAGCCGAGCTGGGTGATTTGCAGTGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCTCGGCGGGTATGGCCATGTCCTGCAACTGGGTAcc                                                                                        <  5:250445/148‑1 (MQ=255)                                                             cATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCg                                                                                       >  1:13718/1‑149 (MQ=255)                                                             cATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCg                                                                                       >  8:115592/1‑149 (MQ=255)                                                                 gCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGa                                                                                   >  1:81163/1‑149 (MQ=255)                                                                  ccGAGCTCGGAGATTTGCAGGGCTGGCTTGAGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCTCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCAGGGTATGGCCATGTCCTCAAACTGCCTACCGTAGAc                                                                                  >  6:133446/1‑149 (MQ=255)                                                                        tGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGctc                                                                            <  2:81163/149‑1 (MQ=255)                                                                          ggTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTg                                                                          >  2:349573/1‑149 (MQ=255)                                                                           gTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACg                                                                                 >  4:284553/1‑141 (MQ=255)                                                                           gTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACg                                                                                 <  3:284553/141‑1 (MQ=255)                                                                               tttGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGcac                                                                      >  5:127774/1‑148 (MQ=255)                                                                                 tGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGCGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGcacca                                                                    <  4:43817/148‑1 (MQ=255)                                                                                  gCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACg                                                                  >  7:85230/1‑149 (MQ=255)                                                                                    aGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTTGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGcc                                                                >  5:141942/1‑149 (MQ=255)                                                                                       gctggctTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTgg                                                                                             <  7:60615/117‑1 (MQ=255)                                                                                       gctggctTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTgg                                                                                             >  8:60615/1‑117 (MQ=255)                                                                                           gctTGCGCAGCTCGACCTTATGGACGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTgg                                                                         <  6:152990/133‑1 (MQ=255)                                                                                           gctTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGAGGGAGTAGCT‑CAC‑CCGCACCTTGCCTTCGCGCGCCACGCCGGGAATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTgg                                                                         >  5:152990/1‑133 (MQ=255)                                                                                             ttGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGGCGGGCTCTGGCACCACGCCGAGGCGCTc                                                       >  2:187357/1‑149 (MQ=255)                                                                                             ttGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTc                                                       >  3:49022/1‑149 (MQ=255)                                                                                                 gcAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATc                                                   <  7:33593/149‑1 (MQ=255)                                                                                                 gcAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATc                                                   <  3:13132/149‑1 (MQ=255)                                                                                                   aGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCAc                                                 >  3:298607/1‑149 (MQ=255)                                                                                                     cTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACtt                                               >  3:282213/1‑149 (MQ=255)                                                                                                        gACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACttt                                              <  4:100655/147‑1 (MQ=255)                                                                                                        gACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACttt                                              <  4:298607/147‑1 (MQ=255)                                                                                                        gACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTgg                                            <  6:127774/149‑1 (MQ=255)                                                                                                        gACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTg                                             <  3:371667/148‑1 (MQ=255)                                                                                                        gACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTg                                             <  3:81465/148‑1 (MQ=255)                                                                                                         aCCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGc                                           <  7:70323/149‑1 (MQ=255)                                                                                                          ccTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGcc                                          <  4:49022/149‑1 (MQ=255)                                                                                                                   cGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCAGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACTCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGt                                 >  1:285610/1‑149 (MQ=255)                                                                                                                        ttGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTc                            >  7:295184/1‑149 (MQ=255)                                                                                                                          gAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCg                           <  4:61266/148‑1 (MQ=255)                                                                                                                           aaTGACAGGGTTTCTCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCGATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCggg                         <  5:133446/149‑1 (MQ=255)                                                                                                                            aTGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCgggg                        <  3:231080/149‑1 (MQ=255)                                                                                                                                 aGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCa                                                  <  6:8332/118‑1 (MQ=255)                                                                                                                                 aGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCa                                                  >  5:8332/1‑118 (MQ=255)                                                                                                                                    gTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAg                >  7:149383/1‑149 (MQ=255)                                                                                                                                       tCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGc              <  5:172268/148‑1 (MQ=255)                                                                                                                                       gtgcgGGTG‑AGCTACAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGAt                      >  6:273351/3‑140 (MQ=255)                                                                                                                                       gtgcgGGTG‑AGCTACAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGAt                      <  5:273351/138‑1 (MQ=255)                                                                                                                                       gtgcgGGTG‑AGCTACAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGCg             >  1:278170/3‑149 (MQ=255)                                                                                                                                              tgtAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGGCGGGGATCGACAGGCGTACCAGc      <  7:115592/149‑1 (MQ=255)                                                                                                                                                tAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGGCGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGCGTACCAGCAt    <  1:187357/149‑1 (MQ=255)                                                                                                                                                 aGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGCGTACCAGCATa   >  6:365262/1‑149 (MQ=255)                                                                                                                                                  gCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGCGTACCAGCATa   <  3:85771/148‑1 (MQ=255)                                                                                                                                                  gCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGCGTACCAGCATAc  >  6:57008/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCAGGCCGATCATTGGTGGTGCCACCCGCGCAGAAAAACGCGTGTCGGGGTGCATCACTCCATAGCCGAGCTGGGTGATTTGCAGGGCTGGCTTGCGCAGCTCGACCTGATGGCCGTCCTTGAATGACAGGGTTTCGCGGGTGTAGCT‑CAC‑CCGCACCTTGCCTTCGGGCGCCACGCCGGGTATGGCCATGTCCTGCAACTGGGTACCGTAGACGGGCTCTGGCACCACGCCGAGGCGCTCGATCACTTTGGCCAGGTAGGGTTGGTCGGGGATCGACAGGCGTACCAGCATAC  >  NC_002947/1052715‑1053008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑