Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,436,199	T→C	6.5%	L45P (CTG→CCG) ‡	gpI →	tail formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,436,199	0	T	C	6.5%	170.8 / 2.8	61	L45P (CTG→CCG) ‡	gpI	tail formation protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base C (2/2);  total (30/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 = 4.36e-01

GAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCG  >  NC_002947/3436054‑3436347                                                                                                                                                   |                                                                                                                                                      gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctggctgg                                                                                                                                                     <  4:11314/148‑1 (MQ=255) gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCtggctggc                                                                                                                                                    <  3:298356/149‑1 (MQ=255)           aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCAGCCC‑TCGACCTTCCATGCCCAGGCCAGCCAGGGCagga                                                                                                                                          <  5:319179/149‑4 (MQ=37)             aTGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATgg                                                                                                                                          >  7:132172/1‑147 (MQ=255)                     gTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctggctgg                                                                                                                                                     >  7:319054/1‑128 (MQ=255)                     gTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctggctgg                                                                                                                                                     <  8:319054/128‑1 (MQ=255)                        aCCGCACAACACCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                                             >  6:355225/1‑149 (MQ=255)                          cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                           <  2:256691/149‑1 (MQ=255)                          cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTCGGCATGGAAGGTCGAGGGCTgg                                                                                                                           >  7:152688/1‑149 (MQ=255)                               aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                                      >  8:92122/1‑149 (MQ=255)                               aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                                      <  7:201739/149‑1 (MQ=255)                                  gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                                   <  2:331140/149‑1 (MQ=255)                                    cACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCt                                                                                                                  >  6:92390/1‑148 (MQ=255)                                      cccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTa                                                                                                                >  6:197777/1‑148 (MQ=255)                                       ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                              >  1:338598/1‑149 (MQ=255)                                       ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                              <  5:92390/149‑1 (MQ=255)                                       ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                              <  8:132172/149‑1 (MQ=255)                                              gAGCGCGCCCTGGAGGTGGCCGACGACCTTGGGCTTGAGCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGGGCTGGCTGGGGTGGGCATGGAAGGTCGTGGGCTGGGAGGCGGCt                                                                                                                  <  8:184519/138‑1 (MQ=255)                                              gAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATACGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAAGGAAGGTCGAGGGCTGGGAGGCGGCt                                                                                                                  >  7:184519/1‑138 (MQ=255)                                                gcgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGAcc                                                                                                     >  2:133000/1‑149 (MQ=255)                                                      ccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgcc                                                                                               <  5:355225/149‑1 (MQ=255)                                                             gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCt                                                                                        >  6:70162/1‑149 (MQ=255)                                                                gccgccGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTg                                                                                                                                                 >  3:289747/1‑89 (MQ=255)                                                                gccgccGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTg                                                                                                                                                 <  4:289747/89‑1 (MQ=255)                                                                  cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATcc                                                                                   >  6:121314/1‑149 (MQ=255)                                                                     cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                                 >  6:232090/1‑148 (MQ=255)                                                                      gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                               <  1:136321/149‑1 (MQ=255)                                                                       aCCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                              <  7:92122/149‑1 (MQ=255)                                                                              ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTc                                                                       >  5:246400/1‑149 (MQ=255)                                                                                  gACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                                     <  5:232090/147‑1 (MQ=255)                                                                                   aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                  <  1:133000/149‑1 (MQ=255)                                                                                   aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGt                                                                   <  5:70162/148‑1 (MQ=255)                                                                                      cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                               <  5:197777/149‑1 (MQ=255)                                                                                      cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                               <  6:246400/149‑1 (MQ=255)                                                                                      cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                               <  5:196016/149‑1 (MQ=255)                                                                                      cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                                >  6:337837/1‑148 (MQ=255)                                                                                           aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                          >  8:49438/1‑149 (MQ=255)                                                                                           aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                          >  7:105605/1‑149 (MQ=255)                                                                                               tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                                         <  1:82390/114‑1 (MQ=255)                                                                                               tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                                         >  2:82390/1‑114 (MQ=255)                                                                                                      cGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACgg                                               <  5:121314/149‑1 (MQ=255)                                                                                                          gCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCg                                            <  7:49438/148‑1 (MQ=255)                                                                                                             gACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcgc                                        <  8:105605/149‑1 (MQ=255)                                                                                                                   tCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGTGCGA‑TTcg                                 >  7:163851/1‑149 (MQ=255)                                                                                                                   tCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTcg                                 >  1:35385/1‑149 (MQ=255)                                                                                                                     aCGCTGCCCGCCGAACTTCC‑TGCCATGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAtt                                                                        <  6:97226/109‑1 (MQ=255)                                                                                                                     aCGCTGCCCGCCGAACTTCC‑TGCCATGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAtt                                                                        >  5:97226/1‑109 (MQ=255)                                                                                                                       gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                              <  5:356191/101‑1 (MQ=255)                                                                                                                       gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                              >  6:356191/1‑101 (MQ=255)                                                                                                                       gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTcggcgg                             >  2:315238/1‑149 (MQ=255)                                                                                                                         tGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTcggcgg                             >  3:131811/1‑147 (MQ=255)                                                                                                                         tGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTcggcgg                             <  4:131811/147‑1 (MQ=255)                                                                                                                             cctcgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                      <  6:345476/132‑1 (MQ=25)                                                                                                                             cctcgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                      >  5:345476/4‑135 (MQ=25)                                                                                                                              ctcgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACAAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTc                       >  7:252213/3‑148 (MQ=37)                                                                                                                                 gAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAgg                   <  2:35385/149‑1 (MQ=255)                                                                                                                                  aaCTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGc                  <  1:315238/149‑1 (MQ=255)                                                                                                                                     ttCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGt                           >  4:311584/1‑137 (MQ=255)                                                                                                                                     ttCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGt                           <  3:311584/137‑1 (MQ=255)                                                                                                                                          tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCgg           >  3:247390/1‑149 (MQ=255)                                                                                                                                          tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCgg           >  2:185798/1‑149 (MQ=255)                                                                                                                                          tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCGTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAGGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCgg           <  2:338598/149‑1 (MQ=255)                                                                                                                                            cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGACTTCGGCGGGTGCTCAAGGCGGTTCgg           <  1:185798/148‑1 (MQ=255)                                                                                                                                            cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCGGGt         >  6:43673/1‑149 (MQ=255)                                                                                                                                            cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGATTCgg           <  5:337837/147‑1 (MQ=255)                                                                                                                                              ctggctggCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCGGGTc        >  3:265477/1‑148 (MQ=255)                                                                                                                                                   tggcCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCGGGTCACcgcg  >  1:312605/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGA‑TTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCG  >  NC_002947/3436054‑3436347

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