Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,000,772	A→C	6.6%	V75G (GTG→GGG) ‡	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,000,772	0	A	C	6.6%	189.0 / 5.2	61	V75G (GTG→GGG) ‡	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base C (2/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTG  >  NC_002947/3000627‑3000920                                                                                                                                                   |                                                                                                                                                     cAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGt                                                                                                                                                   >  4:359136/1‑149 (MQ=255)    tGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCa                                                                                                                                                <  4:212604/149‑1 (MQ=255)               gcgcCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTTGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTg                                                                                                                                     <  2:229649/149‑1 (MQ=255)                 gcCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGc                                                                                                                                   <  5:185274/149‑1 (MQ=255)                     gcagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGcgtcg                                                                                                                               >  7:212083/1‑149 (MQ=255)                        gcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTc                                                                                                                            <  4:28827/149‑1 (MQ=255)                            gTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAaca                                                                                                                        >  7:46975/1‑149 (MQ=255)                            gTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAac                                                                                                                         <  8:212083/148‑1 (MQ=255)                                  cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGt                                                                                                                  >  6:404123/1‑149 (MQ=255)                                   tCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCTACGTCCAGGCTCACATTGGCGTCGTTc                                                                                                                            <  1:12483/139‑1 (MQ=255)                                   tCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTc                                                                                                                            >  2:12483/1‑138 (MQ=255)                                       gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAat                                                                                                             <  5:404123/149‑1 (MQ=255)                                       gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAat                                                                                                             <  1:240260/149‑1 (MQ=255)                                       gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAat                                                                                                             <  8:46975/149‑1 (MQ=255)                                           cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAaca                                                                                                                        >  1:349711/1‑134 (MQ=255)                                           cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAaca                                                                                                                        <  2:349711/134‑1 (MQ=255)                                                     aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                               <  3:359136/149‑1 (MQ=255)                                                      tccggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATc                                                                                                  >  6:70611/1‑145 (MQ=255)                                                      tccggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATc                                                                                                  <  5:70611/145‑1 (MQ=255)                                                      tccggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATg                                                                                              >  7:154184/1‑149 (MQ=255)                                                      tccgggccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATg                                                                                              <  8:173955/149‑1 (MQ=255)                                                         ggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgcgc                                                                                           <  8:129757/149‑1 (MQ=255)                                                                 ccGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCt                                                                                    >  2:369420/1‑148 (MQ=255)                                                                  cGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGa                                                                                  >  8:70627/1‑149 (MQ=255)                                                                  cGAAATCGTTGATGGA‑CTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGAc                                                                                 >  5:162080/1‑149 (MQ=255)                                                                         gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCa                                                                                                                           >  5:337894/1‑101 (MQ=255)                                                                         gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCa                                                                                                                           <  6:337894/101‑1 (MQ=255)                                                                             aTGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCaa                                                                        <  2:65388/148‑1 (MQ=255)                                                                             aTGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCaa                                                                        <  7:70627/148‑1 (MQ=255)                                                                                 gacTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAgcgc                                                                    <  6:162080/146‑1 (MQ=255)                                                                                   cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                                 >  6:330898/1‑149 (MQ=255)                                                                                   cTGAACTCCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                                 >  7:367599/1‑149 (MQ=255)                                                                                        cTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGt                                                            >  6:405417/1‑149 (MQ=255)                                                                                              gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATc                                                        >  4:52553/1‑147 (MQ=255)                                                                                                ggACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAgg                                                    <  3:52553/149‑1 (MQ=255)                                                                                                 gACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGc                                                   >  4:303497/1‑149 (MQ=255)                                                                                                   cGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGcc                                                  >  6:188738/1‑148 (MQ=255)                                                                                                     gCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACATGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                               <  8:66777/149‑1 (MQ=255)                                                                                                      ccAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATc                                                        <  6:376603/139‑1 (MQ=255)                                                                                                      ccAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATc                                                        >  5:376603/1‑139 (MQ=255)                                                                                                       cAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaacgac                                                                                                          >  8:376150/1‑84 (MQ=255)                                                                                                          cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGTCGGTCCTc                                                                          >  7:75879/1‑117 (MQ=255)                                                                                                          cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGTCGGTCCTc                                                                          <  8:81592/117‑1 (MQ=255)                                                                                                          cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGTCGGTCCTc                                                                          >  7:81592/1‑117 (MQ=255)                                                                                                          cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGTCGGTCCTc                                                                          <  8:75879/117‑1 (MQ=255)                                                                                                          cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGc                                          >  3:189027/1‑149 (MQ=255)                                                                                                             ggTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGccc                                        <  3:303497/148‑1 (MQ=255)                                                                                                             ggTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGcc                                         <  8:154184/147‑1 (MQ=255)                                                                                                             ggTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGcc                                         <  5:188738/147‑1 (MQ=255)                                                                                                                 tCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcac                                   <  4:189027/149‑1 (MQ=255)                                                                                                                   cGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCCTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATTACCTGACGCTACTCAAGCGCACCGCAGTAATCCAGGCCGAAGTCGCccagcacc                                  >  8:254083/1‑148 (MQ=25)                                                                                                                              cGAGCCTGGACGGTACCGTC‑ACGTCCAGTCGCACTTTGGCGTCGTTCAACAGGCTGTTGAATGTATGTTTATATAAGCTCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCTAAGTCGCCCAGAACCAGCATCTTGTGc                      <  8:367599/149‑1 (MQ=255)                                                                                                                                 gCCTGGACGTGACGGTA‑CCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATc                            >  1:148156/1‑140 (MQ=21)                                                                                                                                 gCCTGGACGTGACGGTA‑CCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATc                            <  2:148156/140‑1 (MQ=21)                                                                                                                                   cTGGACGTGACGGTA‑CCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg                 >  4:372343/1‑149 (MQ=18)                                                                                                                                   cTGGACGTGACGGTA‑CCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg                 <  3:372343/149‑1 (MQ=18)                                                                                                                                   cTGGACGGTACCGTC‑ACGTCCAGGCTCGCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCgg                    >  8:189409/1‑146 (MQ=255)                                                                                                                                   cTGGACGGTACCGTC‑ACGTCCAGGCTCGCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCa                  <  7:189409/148‑1 (MQ=255)                                                                                                                                         ggTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAc           <  5:330898/149‑1 (MQ=255)                                                                                                                                         ggTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCAACAGCATCTTGTGCGGCAGCTCCAc           <  7:254083/149‑1 (MQ=255)                                                                                                                                               gTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGc     <  5:405417/149‑1 (MQ=255)                                                                                                                                                   aCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTg  >  6:190498/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTG  >  NC_002947/3000627‑3000920

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