Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,425,738	C→T	5.8%	Q92* (CAG→TAG)	PP_4765 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,425,738	0	C	T	5.8%	209.6 / 2.8	69	Q92* (CAG→TAG)	PP_4765	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/36);  new base T (2/2);  total (31/38)
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.93e-01

GTGGGCCGCAAGTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCCACACCCGCCATCGCTGGTT  >  NC_002947/5425594‑5425886                                                                                                                                                     |                                                                                                                                                      gTGGGCCGCAAGTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCAACCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAATGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTg                                                                                                                                                   <  2:4841/149‑1 (MQ=255)    ggCCGCAAGTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGa                                                                                                                                                 <  1:146521/148‑1 (MQ=255)     gCCGCAAGTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGAtt                                                                                                                                               <  2:290085/149‑1 (MQ=255)      ccGCAAGTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGGGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTc                                                                                                                                              <  4:95039/149‑1 (MQ=255)       cGCAAGTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCa                                                                                                                                             >  2:186510/1‑149 (MQ=255)            gTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGtt                                                                                                                                        >  6:286525/1‑149 (MQ=255)                ctgctgGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTgg                                                                                                                                     <  8:8995/148‑1 (MQ=255)                   ctgGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCa                                                                                                                                  <  6:45435/148‑1 (MQ=255)                      gCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATc                                                                                                                                <  2:105915/147‑1 (MQ=255)                      gCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑a                                                                                                                             <  5:255731/149‑1 (MQ=255)                      gCGGGTGTCGGAGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCg                                                                                                                               >  6:255731/1‑148 (MQ=255)                        gggTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑aaa                                                                                                                           >  2:95193/1‑149 (MQ=255)                         ggTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAAc                                                                                                                          <  7:52659/149‑1 (MQ=255)                             tCGGCGGAATGAACATCACCCATTCGGAGCCCAACCCGGGGTTCTTGGCAGGGTATGCGTAACGCCAAGGCGAATTCGTGGCCCTGCTACGCTA‑‑ATTCTATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTtc                                                                                                                      <  1:76527/149‑1 (MQ=255)                               ggcggcATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTtcgt                                                                                                                    <  1:186510/149‑1 (MQ=255)                                gcggcATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTtcgtc                                                                                                                   <  8:185962/149‑1 (MQ=255)                                   gcATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCTAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCgg                                                                                                                 >  1:280209/1‑148 (MQ=255)                                   gcATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCTAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCgg                                                                                                                 >  3:264145/1‑148 (MQ=255)                                   gcATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGc                                                                                                                <  6:137020/149‑1 (MQ=255)                                     aTGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCAc                                                                                                              <  3:100227/149‑1 (MQ=255)                                        aaCATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTc                                                                                                           >  7:69259/1‑149 (MQ=255)                                         aCATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCg                                                                                                          >  8:59719/1‑149 (MQ=255)                                           aTCACCCATTCCGAGCCCTACCCGGCGTTAGTGGCACGCTATGCATAACGCCAAGGCGCAATCGCGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCTGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCggg                                                                                                        >  8:265752/1‑149 (MQ=255)                                            tCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGc                                                                                                       >  4:211012/1‑149 (MQ=255)                                             cACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGAGAAATCGAGCACCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGcc                                                                                                      >  5:214082/1‑149 (MQ=255)                                              aCCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCg                                                                                                     >  6:22284/1‑149 (MQ=255)                                               cgcATTCGGATCCCTACCCGGGGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCgg                                                                                                    <  6:214082/147‑1 (MQ=255)                                               cccATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCgg                                                                                                    >  6:4495/1‑149 (MQ=255)                                                ccATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCGCAAACCTTCGTCGGCACCTCGGGCCgg                                                                                                    <  1:3069/149‑1 (MQ=255)                                                 cATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGgt                                                                                                  <  7:59719/149‑1 (MQ=255)                                                  aTTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGgtg                                                                                                 >  2:49260/1‑149 (MQ=255)                                                      cGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTcc                                                                                             >  7:53333/1‑149 (MQ=255)                                                      cGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTc                                                                                              <  3:50264/148‑1 (MQ=255)                                                      cGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTc                                                                                              <  8:69259/148‑1 (MQ=255)                                                       gAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTccc                                                                                            <  5:4495/149‑1 (MQ=255)                                                            cTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCAcc                                                                                        >  3:188564/1‑148 (MQ=255)                                                                   gCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGACTCTTCGATGGGGACGGCTACAGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACa                                                                                    <  1:95193/149‑1 (MQ=255)                                                                     gTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAg                                                                              <  3:259662/149‑1 (MQ=255)                                                                         gTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGc                                                                          <  4:188564/149‑1 (MQ=255)                                                                                    tGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTtgctgc                                                               <  1:49260/149‑1 (MQ=255)                                                                                      ccGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCg                                                              <  8:53333/148‑1 (MQ=255)                                                                                       cGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTg                                                            >  8:27280/1‑149 (MQ=255)                                                                                        gAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGc                                                           <  3:242358/149‑1 (MQ=255)                                                                                         aaCGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGaaa                                                                               <  7:170706/128‑1 (MQ=255)                                                                                         aaCGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGaaa                                                                               >  8:170706/1‑128 (MQ=255)                                                                                           cGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTg                                                        <  5:102288/149‑1 (MQ=255)                                                                                            gCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTgg                                                       >  7:79984/1‑149 (MQ=255)                                                                                                    gAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAgcg                                               >  8:40001/1‑149 (MQ=255)                                                                                                        tCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTg                                           >  3:267348/1‑149 (MQ=255)                                                                                                           atgCCCTTCTACGCGA‑‑CTTCGATGCGGCCGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCCTCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACCTGAAAGCCGCACCGTTGCTGCGTGCCTGGCTc                                                    <  6:97597/135‑1 (MQ=255)                                                                                                           aGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑ACACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGCCATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTc                                                    >  5:97597/1‑137 (MQ=255)                                                                                                             gCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGa                                      >  1:177822/1‑149 (MQ=255)                                                                                                               ccTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGAc                                     >  3:227739/1‑148 (MQ=255)                                                                                                                    tACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCTAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACggg                               <  2:280209/149‑1 (MQ=255)                                                                                                                      cgcgA‑‑CTTCGATGCGGACGCGT‑‑TGCGCTAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACggg                               <  4:264145/147‑1 (MQ=255)                                                                                                                      cgcgA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACgggg                              >  4:52630/1‑148 (MQ=255)                                                                                                                       gcgA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTtgct                                                                 <  5:138524/112‑1 (MQ=255)                                                                                                                       gcgA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTtgct                                                                 >  6:138524/1‑112 (MQ=255)                                                                                                                       gcgA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACgggg                              >  2:93861/1‑147 (MQ=255)                                                                                                                        cgA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAg                           <  5:22284/149‑1 (MQ=255)                                                                                                                             cTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTa                        <  5:286525/149‑1 (MQ=255)                                                                                                                              ttCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTAt                       >  2:55994/1‑149 (MQ=255)                                                                                                                                 gATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCca                    <  3:211012/149‑1 (MQ=255)                                                                                                                                 gATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCca                    <  7:54209/149‑1 (MQ=255)                                                                                                                                 gATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCca                    >  7:236274/1‑149 (MQ=255)                                                                                                                                  aTGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCcac                   <  7:161743/149‑1 (MQ=255)                                                                                                                                      ggACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTgcgc                                        <  3:285109/124‑1 (MQ=255)                                                                                                                                      ggACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTgcgc                                        >  4:285109/1‑124 (MQ=255)                                                                                                                                            gT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCCACACCCGCCATc         >  2:11171/1‑149 (MQ=255)                                                                                                                                                     cAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTtgct                                                                 <  2:46192/86‑1 (MQ=255)                                                                                                                                                     cAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTtgct                                                                 >  1:46192/1‑86 (MQ=255)                                                                                                                                                     cAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCCACACCCGCCATCGCTgggt  >  4:272091/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTGGGCCGCAAGTTCCTGCTGGCGGGTGTCGGCGGCATGAACATCACCCATTCCGAGCCCTACCCGGCGTTCGTGGCACGCTATGCCGAACGCCAAGGCGAAATCGAGGCCCTGCTACGCGA‑‑CTTCGATGCGGACGCGT‑‑TGCGCCAGTGGATTCATGGTTTGGGCATCG‑AAACCTTCGTCGGCACCTCGGGCCGGGTGTTCCCCACCGACATGAAAGCCGCACCGTTGCTGCGTGCCTGGCTCAAGCGCCTGCGCGACAACGGGGTAGTTATCCACACCCGCCATCGCTGGTT  >  NC_002947/5425594‑5425886

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑