Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,547,156	A→G	6.4%	L633P (CTC→CCC)	rnr ←	exoribonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,547,156	0	A	G	6.4%	184.1 / 2.4	62	L633P (CTC→CCC)	rnr	exoribonuclease R
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (2/2);  total (31/31)
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.23e-01

GACGAACAGGCCGAAACCGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGATGCACCAGCAGGTCCGGGTAAC  >  NC_002947/5547008‑5547297                                                                                                                                                     |                                                                                                                                               gacgaACAGGCCGAAACCGGTCACCGCAGGGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGa                                                                                                                                                <  7:182602/149‑1 (MQ=255)  acgaACAGGCCGAAACCGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAg                                                                                                                                               <  1:54322/149‑1 (MQ=255)    gaACAGGCCGAAACCGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGtt                                                                                                                                             >  6:169288/1‑149 (MQ=255)      aCAGGCCGAAACCGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGc                                                                                                                                           >  2:223655/1‑149 (MQ=255)       cAGGCCGAAACCGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCt                                                                                                                                          >  3:166851/1‑149 (MQ=255)              aaaCCGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGgt                                                                                                                                  <  1:226459/149‑1 (MQ=255)                 cccgTCACCGCAGTGATCACACCCGGGAAGCTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTCAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGTCCCGGCGCTCGGTCATCGAGCACTTCTCGCCCAGATGATC‑CAGGctgt                                                                                                                                >  8:218393/4‑148 (MQ=255)                 ccGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGgtgt                                                                                                                                >  6:41506/1‑148 (MQ=255)                            aGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACg                                                                                                                    >  5:68030/1‑149 (MQ=255)                                   acacCCCGGCAGGTCGCGCCCACGCGGTCTTTCATGACCTCGCACTTGCGCCAGTTGACCACGGCACGGGTGGCGTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTGCTCGTCGTACGGGTAGAt                                                                                                             <  4:105256/149‑1 (MQ=255)                                     acCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGATCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAggg                                                                                                                                   <  6:37432/125‑1 (MQ=255)                                     acCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAggg                                                                                                                                   >  5:37432/1‑125 (MQ=255)                                     acCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAgg                                                                                                                                    <  7:239705/124‑1 (MQ=255)                                     acCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAgg                                                                                                                                    >  8:239705/1‑124 (MQ=255)                                       ccGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGAt                                                                                                             >  7:247477/1‑145 (MQ=255)                                       ccGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGAt                                                                                                             <  8:247477/145‑1 (MQ=255)                                          ggAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGcc                                                                                                      >  4:45827/1‑149 (MQ=255)                                             aGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTg                                                                                                   <  3:45827/149‑1 (MQ=255)                                               gTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTggg                                                                                                 <  4:262237/149‑1 (MQ=255)                                                 ctcGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATAcgcg                                                                                                        >  6:261272/1‑140 (MQ=255)                                                 ctcGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATAcgcg                                                                                                        <  5:261272/140‑1 (MQ=255)                                                 ctcGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGaa                                                                                               >  2:68814/1‑149 (MQ=255)                                                    gCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGaa                                                                                               <  7:206441/146‑1 (MQ=255)                                                    gCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGc                                                                                            >  1:115792/1‑149 (MQ=255)                                                        aCGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCAt                                                                                        <  1:68814/149‑1 (MQ=255)                                                                tttCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTGGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCgg                                                                                >  2:157396/1‑149 (MQ=255)                                                                  tCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCa                                                                              >  6:128963/1‑149 (MQ=255)                                                                    aTGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTg                                                                                                   >  1:119022/1‑126 (MQ=255)                                                                    aTGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTg                                                                                                   <  2:119022/126‑1 (MQ=255)                                                                    aTGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCAc                                                                             <  2:15761/148‑1 (MQ=255)                                                                    aTGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCCGCACg                                                                            >  7:18779/1‑149 (MQ=255)                                                                     tGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGc                                                                           <  5:169288/149‑1 (MQ=255)                                                                          tCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGGTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGAc                                                                      >  4:200509/1‑149 (MQ=255)                                                                          tCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATg                                                                                       <  1:758/132‑1 (MQ=255)                                                                          tCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATg                                                                                       >  2:758/1‑132 (MQ=255)                                                                            gCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGAc                                                                      <  3:264843/147‑1 (MQ=255)                                                                                  gAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCggg                                                               <  5:128963/148‑1 (MQ=255)                                                                                    gCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCtt                                                                                                    >  4:184130/1‑109 (MQ=255)                                                                                    gCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCtt                                                                                                    <  3:184130/109‑1 (MQ=255)                                                                                      cAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATc                                                          <  4:166851/149‑1 (MQ=255)                                                                                       aGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGGGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATc                                                          >  2:52206/1‑148 (MQ=255)                                                                                       aGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGGGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCg                                                         >  2:224462/1‑149 (MQ=255)                                                                                         ttGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGGGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGAc                                                       <  1:52206/149‑1 (MQ=255)                                                                                         ttGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGGGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGAc                                                       <  1:224462/149‑1 (MQ=255)                                                                                           gACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCgg                                                                <  3:17976/138‑1 (MQ=255)                                                                                           gACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCgg                                                                >  4:17976/1‑138 (MQ=255)                                                                                               aCGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGa                                                 <  6:68030/149‑1 (MQ=255)                                                                                                cGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCg                                                         >  1:164984/1‑140 (MQ=255)                                                                                                cGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCg                                                         <  2:164984/140‑1 (MQ=255)                                                                                                cGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAg                                                <  2:115792/149‑1 (MQ=255)                                                                                                 gTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGc                                               >  3:204967/1‑149 (MQ=255)                                                                                                      gggTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGgatg                                          <  5:41506/149‑1 (MQ=255)                                                                                                        gTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGgatgat                                        >  6:112925/1‑149 (MQ=255)                                                                                                               cgtcgGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGAtgctgc                                   <  1:23013/147‑1 (MQ=255)                                                                                                                gtcgGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGa                                >  3:158395/1‑149 (MQ=255)                                                                                                                      ccGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCg                                  >  7:104748/1‑141 (MQ=255)                                                                                                                      ccGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCg                                  <  8:104748/141‑1 (MQ=255)                                                                                                                         gcgcTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGAt                       <  4:158395/149‑1 (MQ=255)                                                                                                                                 tCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGATGCACcagc               <  7:106761/149‑1 (MQ=255)                                                                                                                                   aTCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGATGCACcagca              <  4:5539/148‑1 (MQ=255)                                                                                                                                     cGAGCACTGCTCGCCGAGTTGCTCTCAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGATGCACCAGCAgg            >  7:167282/1‑149 (MQ=255)                                                                                                                                       aGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGATGCACCAGCAGGTcc         >  7:113941/1‑149 (MQ=255)                                                                                                                                              cTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGATGCACCAGCAGGTCCGGGTAac  >  4:151229/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               GACGAACAGGCCGAAACCGGTCACCGCAGTGATCACACCCGGGAAGGTCTCGCCCACGCGGTCTTTCATGAACTCGCACTTGAGCCAGTTGACCACGTCACGGGTGGCTTCGTCGGCCCGGCGCTCGGTCATCGAGCACTGCTCGCCGAGTTGCTC‑CAGGGTGTTCTCGTCGTACGGGTAGATACGCGCCTTGGGAATGCTCATGGCGCCGGCACGCTTGACGTGCGGGGTATCGACCTTGGAGCGGATGATGCTGCGGATGGCGCGATGCACCAGCAGGTCCGGGTAAC  >  NC_002947/5547008‑5547297

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