Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,884,981	G→C	6.6%	A53G (GCA→GGA) ‡	PP_5159 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,884,981	0	G	C	6.6%	225.7 / 4.3	72	A53G (GCA→GGA) ‡	PP_5159	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (3/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

AGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGG  >  NC_002947/5884837‑5885127                                                                                                                                                 |                                                                                                                                                   aGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCggg                                                                                                                                                <  5:369952/149‑1 (MQ=255)    gcTTGGTTGATTGCCAGGCGCAGCTGATCCGTCGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGGCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                             <  2:229883/148‑1 (MQ=255)    ccTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                             <  1:319639/149‑1 (MQ=255)     cTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTg                                                                                                                                            <  1:133530/149‑1 (MQ=255)     cTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCGCAGGGGTTg                                                                                                                                            <  3:131461/149‑1 (MQ=255)       tGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAg                                                                                                                                          >  7:260039/1‑149 (MQ=255)         gTTGATTGCCAGGCGCAGCTGATCCATGGCGGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCTCGAATGGTTTCGCTCACCTCGTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTTAGt                                                                                                                                         <  7:107355/148‑1 (MQ=255)         gTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGtt                                                                                                                                        <  2:13443/149‑1 (MQ=255)          ttGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCTCAACCCCGCAGGGGTTgggcgg                                                                                                                                       <  5:93291/149‑6 (MQ=255)          ttGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTc                                                                                                                                       <  1:285282/149‑1 (MQ=255)          ttGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTc                                                                                                                                       >  6:205400/1‑149 (MQ=255)             aTTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCT‑GTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTg                                                                                                                                                    >  8:245024/1‑132 (MQ=255)             aTTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCT‑GTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTg                                                                                                                                                    <  7:245024/132‑1 (MQ=255)                      cgcAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATc                                                                                                                           <  4:185040/149‑1 (MQ=255)                       gcAGCTGACCCATGGCCGCAAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCa                                                                                                                          <  6:308078/149‑1 (MQ=255)                        cAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCAc                                                                                                                         >  5:62722/1‑149 (MQ=255)                         aGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACg                                                                                                                        >  6:143621/1‑149 (MQ=255)                               tCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCg                                                                                                                  <  8:260039/149‑1 (MQ=255)                                 cATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCgg                                                                                                                 >  2:149369/1‑148 (MQ=255)                                      ccgccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCg                                                                                                           >  7:254128/1‑149 (MQ=255)                                       cgccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCgg                                                                                                          <  5:332963/149‑1 (MQ=255)                                           aGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGc                                                                                                      <  1:149369/149‑1 (MQ=255)                                             gTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGc                                                                                                      >  1:269308/1‑147 (MQ=255)                                                   gTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcg                                                                                               <  6:62722/148‑1 (MQ=255)                                                        aGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCAc                                                                                          >  2:235104/1‑148 (MQ=255)                                                          gCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAg                                                                                       <  7:214023/149‑1 (MQ=255)                                                                  cTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggcag                                                                               >  2:77053/1‑149 (MQ=255)                                                                         gggTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGaaa                                                                        <  1:77053/149‑1 (MQ=255)                                                                               aTCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCGCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTgcgc                                                                  <  6:119115/149‑1 (MQ=255)                                                                                tCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGc                                                                                                                                     <  8:80086/81‑1 (MQ=255)                                                                                tCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGc                                                                                                                                     >  7:80086/1‑81 (MQ=255)                                                                                tCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATc                                                                                                                           >  1:104124/1‑91 (MQ=255)                                                                                tCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATc                                                                                                                           <  2:104124/91‑1 (MQ=255)                                                                                  cACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCct                                                                <  7:88342/148‑1 (MQ=255)                                                                                       gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTgg                                                           >  7:236354/1‑148 (MQ=255)                                                                                       gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTAGCCTTCGGCATCCCGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGCAAGTGCGCCTGCTgg                                                           >  7:134618/1‑148 (MQ=255)                                                                                          cTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTAGGCAGGAAGAAAGTGCGACTGCTGGGAc                                                        >  7:115921/1‑148 (MQ=255)                                                                                            ggTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTgg                                                           >  8:291755/1‑143 (MQ=255)                                                                                            ggTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTgg                                                           <  7:291755/143‑1 (MQ=255)                                                                                            ggTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAg                                                     >  8:57223/1‑149 (MQ=255)                                                                                            ggTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAg                                                     >  7:323554/1‑149 (MQ=255)                                                                                             gTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGAc                                                        <  5:304829/145‑1 (MQ=255)                                                                                             gTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGAc                                                        >  6:304829/1‑145 (MQ=255)                                                                                                  cgcgAAGGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCGGCTGGGACCAGCCATgg                                               <  8:269015/149‑1 (MQ=255)                                                                                                    cgAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGt                                             >  8:322522/1‑149 (MQ=255)                                                                                                    cgAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGt                                             <  8:236354/149‑1 (MQ=255)                                                                                                       aTGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                               >  6:168874/1‑144 (MQ=255)                                                                                                       aTGGTTTCGCGCACCGCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                               <  5:168874/144‑1 (MQ=255)                                                                                                        tGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGcc                                         <  5:143621/149‑1 (MQ=255)                                                                                                           tttCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGa                                      <  8:254128/149‑1 (MQ=255)                                                                                                               gcgcACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgc                                                                                                >  2:335682/1‑87 (MQ=255)                                                                                                                           ggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGc                        <  7:322522/147‑1 (MQ=255)                                                                                                                               gCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTcga                               <  8:231406/136‑1 (MQ=255)                                                                                                                               gCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTcga                               >  7:231406/1‑136 (MQ=255)                                                                                                                                   cccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATc              >  2:45774/1‑149 (MQ=255)                                                                                                                                    gaactcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCg             >  8:65439/6‑149 (MQ=255)                                                                                                                                    ccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCg             >  2:2882/1‑149 (MQ=255)                                                                                                                                      gCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCggg           >  3:7342/1‑149 (MQ=255)                                                                                                                                      actcAACCCCGCAGGGGTTGGGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCggg           >  8:289657/4‑149 (MQ=255)                                                                                                                                        tcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGCCGGCCGCAAGCATCGGgtg         >  7:279893/2‑149 (MQ=255)                                                                                                                                        ccAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAACGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgtg         >  3:205061/1‑149 (MQ=255)                                                                                                                                        ccAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgtg         >  4:174812/1‑149 (MQ=255)                                                                                                                                        ccAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgtg         >  3:21059/1‑149 (MQ=255)                                                                                                                                         cAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg            <  6:331879/145‑1 (MQ=255)                                                                                                                                         cAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg            >  5:331879/1‑145 (MQ=255)                                                                                                                                          aaCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGtt       >  1:352722/1‑149 (MQ=255)                                                                                                                                             cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCggg           >  7:310012/1‑142 (MQ=255)                                                                                                                                             cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCggg           <  8:310012/142‑1 (MQ=255)                                                                                                                                              ccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAgg   <  1:45774/149‑1 (MQ=255)                                                                                                                                               cTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGcc                                         >  4:211074/1‑110 (MQ=255)                                                                                                                                               cTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGcc                                         <  3:211074/110‑1 (MQ=255)                                                                                                                                               cTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAggg  <  7:57223/149‑1 (MQ=255)                                                                                                                                               cTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAggg  <  3:174812/149‑1 (MQ=255)                                                                                                                                               cTGCGGGGTTGAGTTCGCATTCGGCATCACGGCTCCGGCACTCGTCTGCAAGTGCGCGCTCCAGGTCGGCAGGCAGAAAGTGGGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAggg  <  8:134618/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGG  >  NC_002947/5884837‑5885127

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