Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	944,831	G→A	8.3%	P264L (CCC→CTC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	944,831	0	G	A	8.3%	171.2 / 4.5	60	P264L (CCC→CTC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (25/30);  new base A (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCG  >  NC_002947/944687‑944977                                                                                                                                                 |                                                                                                                                                    ccAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTGCCGAGcc                                                                                                                                                  <  2:239705/148‑2 (MQ=255)  cAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCg                                                                                                                                                >  1:4146/1‑149 (MQ=255)            tGCGGCGCTGATCTCCACGGCCAGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTGCCGAGCCCGGCTCGAcg                                                                                                                                        <  5:249403/147‑2 (MQ=255)                          ccAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTAccc                                                                                                                        >  7:59884/1‑149 (MQ=255)                           cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCgg                                                                                                                                                     <  6:57482/119‑1 (MQ=255)                           cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCgg                                                                                                                                                     >  5:57482/1‑119 (MQ=255)                           cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGtt                                                                                                                            <  2:86103/144‑1 (MQ=255)                           cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGtt                                                                                                                            <  2:139108/144‑1 (MQ=255)                            aGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAATCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTg                                                                                                                      <  1:263183/149‑1 (MQ=255)                            aGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAgg                                                                                                                                 >  8:218975/1‑138 (MQ=255)                            aGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAgg                                                                                                                                 <  7:218975/138‑1 (MQ=255)                              ggCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCt                                                                                                                    <  2:4146/149‑1 (MQ=255)                                      cGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATa                                                                                                            <  3:206460/149‑1 (MQ=255)                                      cGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATa                                                                                                            <  8:99576/149‑1 (MQ=255)                                              tgtgGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAggtc                                                                                                    >  5:148399/1‑149 (MQ=255)                                              tgtgGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAggt                                                                                                     <  8:59884/148‑1 (MQ=255)                                                 ggCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGc                                                                                                                                              >  4:119153/1‑104 (MQ=255)                                                 ggCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGc                                                                                                                                              <  3:119153/104‑1 (MQ=255)                                                    ggTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGtcggtcgg                                                                                              >  5:53387/1‑149 (MQ=255)                                                         cAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGgcgc                                                                                           >  2:165565/1‑147 (MQ=255)                                                          aGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTGCCGAGCCCGGCTCGAcg                                                                                                                                        >  4:93437/1‑100 (MQ=255)                                                            ggTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCAt                                                                                                                 <  2:356423/122‑1 (MQ=255)                                                            ggTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCAt                                                                                                                 >  1:356423/1‑122 (MQ=255)                                                                 tGCTGAACTGGCCGTCCGCACCGGCCACCGCAGGGCCCAGCAATACGCCACCTGCGCCTCGCACGGAACCCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCa                                                                                 <  1:190366/149‑1 (MQ=255)                                                                  gCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAg                                                                                >  1:293566/1‑149 (MQ=255)                                                                    tgcggtGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCAGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGtt                                                                              >  3:73560/6‑149 (MQ=255)                                                                      aaCTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGtt                                                                              <  1:355461/147‑1 (MQ=255)                                                                           gCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCg                                                                       <  2:293566/149‑1 (MQ=255)                                                                             cGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGt                                                                     <  7:306625/149‑1 (MQ=255)                                                                              gTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTg                                                                    >  2:90834/1‑149 (MQ=255)                                                                              gTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTg                                                                    <  2:289764/149‑1 (MQ=255)                                                                                ccGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTg                                                                    >  7:142072/1‑147 (MQ=255)                                                                                  gCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGt                                                                <  4:73560/149‑1 (MQ=255)                                                                                   cACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtg                                                                                      >  6:374288/1‑126 (MQ=255)                                                                                   cACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtg                                                                                      <  5:374288/126‑1 (MQ=255)                                                                                   cACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGt                                                                >  5:310653/1‑148 (MQ=255)                                                                                       ggCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑g                                                          <  8:142072/149‑1 (MQ=255)                                                                                        gCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑gc                                                         >  8:227506/1‑149 (MQ=255)                                                                                           aCCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑gcgc                                                       >  2:18272/1‑148 (MQ=255)                                                                                                     ccAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGgtg                                                                                                                              >  6:71404/1‑68 (MQ=255)                                                                                                     ccAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGgtg                                                                                                                              <  5:71404/68‑1 (MQ=255)                                                                                                     ccAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGcc                                               >  1:198229/1‑146 (MQ=255)                                                                                                      cAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCgaga                                           >  4:212071/1‑149 (MQ=255)                                                                                                       aGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGAc                                          <  1:165565/149‑1 (MQ=255)                                                                                                       aGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGAc                                          <  6:148399/149‑1 (MQ=255)                                                                                                            tACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTAc                                     <  1:90834/149‑1 (MQ=255)                                                                                                              cGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTAcc                                    <  5:328352/148‑1 (MQ=255)                                                                                                                ccACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACcggc                                 >  6:242733/1‑149 (MQ=255)                                                                                                                       cgcCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGGGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCgtc                           <  7:138286/148‑1 (MQ=255)                                                                                                                        gcCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCgtcag                         >  7:322445/1‑149 (MQ=255)                                                                                                                         ccTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACcg                                   >  3:364588/1‑138 (MQ=255)                                                                                                                         ccTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACcg                                   <  4:364588/138‑1 (MQ=255)                                                                                                                           tCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCgtcagt                        >  7:260467/1‑147 (MQ=255)                                                                                                                             gCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCgtcagtca                      <  2:198229/147‑1 (MQ=255)                                                                                                                               ctggAATCCGTCGAGCCGGGCTCGGCACTACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGTCGCCGAGACATTACCGGCGGAGTCAGTCAggg                   <  2:68704/146‑1 (MQ=255)                                                                                                                               aCGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGt                  <  3:212071/149‑1 (MQ=255)                                                                                                                                 ggAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGa                >  3:162625/1‑149 (MQ=255)                                                                                                                                  ttaccggTCGTGCCGAGCCCGGCTCGACGGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTTGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGa                >  5:367499/7‑148 (MQ=255)                                                                                                                                  gAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGAGTTACCCTGCTCATTGATAGACAGGTAGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTAAGTCAggg                   >  1:68704/1‑145 (MQ=255)                                                                                                                                      ccGTCGAGCCGGGCGCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGTCATGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATc            <  1:273947/149‑1 (MQ=255)                                                                                                                                        gTCGTGCCGAGCCCGGCTCGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAg         >  5:147998/1‑149 (MQ=255)                                                                                                                                               cGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCg  <  5:242733/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACA‑GCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCG  >  NC_002947/944687‑944977

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