Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,224,832	C→G	6.7%	E601Q (GAG→CAG) ‡	PP_1963 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,224,832	0	C	G	6.7%	206.9 / 4.3	71	E601Q (GAG→CAG) ‡	PP_1963	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base G (2/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATT  >  NC_002947/2224686‑2224979                                                                                                                                                    |                                                                                                                                                    tAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCaa                                                                                                                                                   <  3:94789/149‑1 (MQ=255)      gggCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGa                                                                                                                                              <  5:325072/149‑1 (MQ=255)        gCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGc                                                                                                                                            <  3:261308/149‑1 (MQ=255)         cGTTCCGCCGTCTATATTGCGAAGCCTTCTTGTTTGTGTCTTCATCATGATTGTCGTTTATGAGAAACTCCAAAATTTTAGCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCt                                                                                                                                           <  7:50290/149‑1 (MQ=255)          gTTCCGCCGTCTATATCGCGAAGCATGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCtg                                                                                                                                          <  7:268116/149‑2 (MQ=255)             ccgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGAc                                                                                                                                       >  2:15505/1‑149 (MQ=255)              cgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAg                                                                                                                                             >  3:22400/1‑142 (MQ=255)              cgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAg                                                                                                                                             <  4:22400/142‑1 (MQ=255)              cgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACg                                                                                                                                      >  1:146981/1‑149 (MQ=255)                     tataTCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCAt                                                                                                                                >  8:299168/1‑148 (MQ=255)                       tatCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGc                                                                                                                                            <  8:296353/134‑1 (MQ=255)                       tatCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATc                                                                                                                               <  8:287758/147‑1 (MQ=255)                         tCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                            >  2:362092/1‑148 (MQ=255)                            cgAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAgg                                                                                                                        <  1:241836/149‑1 (MQ=255)                            cgAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAgg                                                                                                                        >  4:299974/1‑149 (MQ=255)                               aGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTa                                                                                                                     <  1:362092/149‑1 (MQ=255)                                        gATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCtt                                                                                                            >  2:151659/1‑149 (MQ=255)                                          tGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                           <  2:125892/148‑1 (MQ=255)                                            gtttCTTCATCATGATTGTCGTTGAAGATAAACACCAAAATTTTACCTTCCCACCGCAGGTCAGAGGCCCGTCCAA‑CTACGCTTTCATTTTCCACCTCATACTCAAAATGAGCTAAACGCTCCATCAAATTGTGTAGTGAGTCttaga                                                                                                         >  8:96095/4‑145 (MQ=255)                                            gtgtCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGt                                                                                                         <  1:15505/148‑1 (MQ=255)                                                  tcatcaTGATTGTCGTTGGTGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                   >  8:260274/1‑148 (MQ=255)                                                   catcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                   <  5:168120/147‑1 (MQ=255)                                                         gATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                            >  3:267403/1‑116 (MQ=255)                                                         gATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                            <  4:267403/116‑1 (MQ=255)                                                          aTTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTggg                                                                                          >  8:242717/1‑149 (MQ=255)                                                              tCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGagca                                                                                      >  5:254631/1‑149 (MQ=255)                                                                     tgaGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGAt                                                                                <  7:299168/148‑1 (MQ=255)                                                                          aaCTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACtt                                                                          <  7:242717/149‑1 (MQ=255)                                                                          aaCTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACtt                                                                          <  2:107428/149‑1 (MQ=255)                                                                             tCCAAAATTTTACCTTCCGCGCGCAGGTCTGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCTTCAAATAGGGTAGTGAGTCTTCGTAAATTCGaa                                                                                                <  8:248098/124‑1 (MQ=255)                                                                             tCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATa                                                                                                                          <  4:228413/98‑1 (MQ=255)                                                                             tCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATa                                                                                                                          >  3:228413/1‑98 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGCCGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                           >  5:68688/1‑127 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                           <  5:329428/111‑1 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                           >  6:329428/1‑111 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                           <  6:68688/127‑1 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                     >  2:12248/1‑149 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                     >  1:290570/1‑149 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                     >  7:129495/1‑149 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                     >  4:103752/1‑149 (MQ=255)                                                                               cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                     >  4:339966/1‑149 (MQ=255)                                                                                aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGa                                                                                                                 <  5:91150/104‑1 (MQ=255)                                                                                aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGa                                                                                                                 >  6:91150/1‑104 (MQ=255)                                                                                aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                           >  1:203524/1‑126 (MQ=255)                                                                                aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                           <  2:203524/126‑1 (MQ=255)                                                                                aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGagcagc                                                                                    >  7:129397/1‑133 (MQ=255)                                                                                aaaaTTTTACCTTACCAGCGTAGGTCAGAGGACAGTGCAA‑CTTCGTTTTCATTTTGCAGCTCTTACTCAATTTCTGCTAGACGCTCCATCATATTGGGTAGTGAGTCTTCGTAAATTCGAAGCTGGGagcagc                                                                                    <  8:129397/133‑1 (MQ=255)                                                                                           ttCCCAGCGCAGGTCAGAGGACCGTGCACACTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCAt                                                          >  1:211263/1‑149 (MQ=255)                                                                                             cccAGCGCAGGTCGTAGGACCGTCCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGATCCATCAAATAGGGTAGTGAGTCTTCGTAAATTAGAACCAGGGAGCAGCCGATGCACTTGTCTAAGCTGTAGCAtttc                                                       >  6:7380/1‑149 (MQ=255)                                                                                                 gcgcAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCAtttctttc                                                   >  8:283291/1‑149 (MQ=255)                                                                                                   gcAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCa                                                 <  7:260274/149‑1 (MQ=255)                                                                                                      ggTCATAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCTTCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              <  7:96095/149‑1 (MQ=255)                                                                                                      ggTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              <  1:12248/149‑1 (MQ=255)                                                                                                      ggTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              <  3:339966/149‑1 (MQ=255)                                                                                                        tCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAgg                                            <  3:299974/149‑1 (MQ=255)                                                                                                          agagGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              <  4:188538/145‑1 (MQ=255)                                                                                                          agagGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              >  3:188538/1‑145 (MQ=255)                                                                                                            agGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAAcc                                        >  8:84101/1‑149 (MQ=255)                                                                                                            agGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAAcc                                        <  1:151659/149‑1 (MQ=255)                                                                                                              gACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGa                                      <  7:283291/149‑1 (MQ=255)                                                                                                                ccGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGAtt                                    <  2:290570/149‑1 (MQ=255)                                                                                                                       c‑cTTCGCTTTCATTTTCCATCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGaa                             <  7:224661/149‑1 (MQ=255)                                                                                                                       c‑cTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGaa                             >  3:24166/1‑149 (MQ=255)                                                                                                                              cTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAATTTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCtt                       <  3:103752/149‑1 (MQ=255)                                                                                                                               tttCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCCAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTc                      >  7:334002/1‑149 (MQ=255)                                                                                                                                   aTTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGACCCTGGGAGCAGCCGATGCACTTGTCTATGGTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGt                  >  8:105298/1‑149 (MQ=255)                                                                                                                                        ctaGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGt             >  5:357829/3‑149 (MQ=255)                                                                                                                                         taGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCTTTTCTTTCCATTAGGACCCGATTGCTCGAATCACTTCCAGTGTGGtg            >  8:85525/2‑148 (MQ=255)                                                                                                                                         cAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGtt            <  6:254631/149‑1 (MQ=255)                                                                                                                                          aGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGtt            >  6:32029/1‑148 (MQ=255)                                                                                                                                                   tCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGAtt  >  3:25173/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATT  >  NC_002947/2224686‑2224979

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