Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,000,765	T→G	6.2%	V77V (GTA→GTC)	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,000,765	0	T	G	6.2%	201.2 / 4.6	64	V77V (GTA→GTC)	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base G (2/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAC  >  NC_002947/3000620‑3000911                                                                                                                                                  |                                                                                                                                                    gggCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTAcc                                                                                                                                                  >  7:66777/1‑149 (MQ=255)   gCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGt                                                                                                                                                <  6:193088/149‑1 (MQ=255)        cAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGt                                                                                                                                          >  4:359136/1‑149 (MQ=255)           tGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCa                                                                                                                                       <  4:212604/149‑1 (MQ=255)              ttCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGt                                                                                                                                                <  7:228434/138‑1 (MQ=255)              ttCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGt                                                                                                                                                >  8:228434/1‑138 (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                                                                                                    <  8:46975/149‑1 (MQ=255)                                              gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAat                                                                                                    <  1:240260/149‑1 (MQ=255)                                                  cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAaca                                                                                                               <  2:349711/134‑1 (MQ=255)                                                  cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAaca                                                                                                               >  1:349711/1‑134 (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                                                               <  7:70627/148‑1 (MQ=255)                                                                                    aTGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCaa                                                               <  2:65388/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                                                                 <  8:75879/117‑1 (MQ=255)                                                                                                                 cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGTCGGTCCTc                                                                 >  7:81592/1‑117 (MQ=255)                                                                                                                 cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGTCGGTCCTc                                                                 <  8:81592/117‑1 (MQ=255)                                                                                                                 cTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGTCGGTCCTc                                                                 >  7:75879/1‑117 (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                   <  2:148156/140‑1 (MQ=21)                                                                                                                                        gCCTGGACGTGACGGTA‑CCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATc                   >  1:148156/1‑140 (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)                                                                                                                                                  |                                                                                                                                                    GGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTC‑ACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAC  >  NC_002947/3000620‑3000911

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