Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,563,123	C→T	88.0%	Q289* (CAA→TAA)	PP_1372 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,563,123	0	C	T	88.0%	134.5 / 9.3	58	Q289* (CAA→TAA)	PP_1372	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (26/25);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

TTGTGGCTGCTGGCCGAACTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAGGAGCAGAAGTCGCT  >  NC_002947/1562988‑1563271                                                                                                                                         |                                                                                                                                                      ttgtGGCTGCTGGCCGAACTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGc                                                                                                                                           <  6:140560/148‑1 (MQ=255)   gtgGCTGCTGGCCGAACTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGa                                                                                                                                        >  4:55971/1‑149 (MQ=255)    tgGCTGCTGGCCGAACTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGAt                                                                                                                                       <  1:248132/149‑1 (MQ=255)      gctgctGGCCGAACTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATc                                                                                                                                      <  6:140974/148‑1 (MQ=255)                 aaCTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAgggg                                                                                                                          >  5:132524/1‑149 (MQ=255)                 aaCTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTAA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAgggg                                                                                                                          >  4:77572/1‑149 (MQ=255)                  aCTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCACTCGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAgggg                                                                                                                          <  1:289773/149‑1 (MQ=255)                  aCTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGc                                                                                                                                           <  3:65318/131‑1 (MQ=255)                  aCTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGc                                                                                                                                           >  4:65318/1‑131 (MQ=255)                  aCTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGATCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCGGATCCGCTCCAAggggg                                                                                                                         >  7:331918/1‑149 (MQ=255)                    tGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACATGCTCTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATACGCTCCAAgggggtg                                                                                                                       >  2:89439/1‑149 (MQ=255)                      ttCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGATCCGCTCCAAgggggtggg                                                                                                                     >  3:37435/1‑149 (MQ=255)                         gAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCt                                                                                                                                          <  4:322570/125‑1 (MQ=255)                         gAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCt                                                                                                                                          >  3:322570/1‑125 (MQ=255)                              gTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTAtt                                                                                                             >  7:24149/1‑149 (MQ=255)                               ttGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTAttt                                                                                                            <  3:55971/149‑1 (MQ=255)                                 gCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGc                                                                                                                                           >  1:14264/1‑116 (MQ=255)                                 gCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGc                                                                                                                                           <  2:14264/116‑1 (MQ=255)                                 gCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCAGATCCGCTCCAAGGGGGTGGGGGTCTATTTCg                                                                                                          >  2:118189/1‑149 (MQ=255)                                  cccGAGCGCGGGGATGCCGACTTGCGGGTGGTGGCGGTGATCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGGCCAAGGGGTTGCAGGAGCGTCTTGAGTAAGGGGTCCGGCTGATCCGCTCCTAGGGGGTGGGGGTCTATTTCg                                                                                                          <  2:32577/148‑1 (MQ=255)                                   ccGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGt                                                                                                         >  2:145700/1‑148 (MQ=255)                                          cgGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTcc                                                                                                                                              >  8:47410/1‑104 (MQ=255)                                           ggTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTc                                                                                               >  4:150820/1‑149 (MQ=255)                                               aTGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCg                                                                                           >  4:309782/1‑149 (MQ=255)                                               aTGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCg                                                                                           <  4:37435/149‑1 (MQ=255)                                                tGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTc                                                                                                           >  3:205512/1‑133 (MQ=255)                                                tGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTc                                                                                                           <  4:205512/133‑1 (MQ=255)                                                  ccGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGt                                                                                        <  6:100639/149‑1 (MQ=255)                                                   cGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGt                                                                                        >  6:170358/1‑148 (MQ=255)                                                         tccGGTGCAGCCGCTGTGCGTCGATGAAGCACACCTGCTGGTCA‑ACGCCACGCCCAAGGCGTTGCAGGACCGTCTTGCGCAAGAGGTCCGGCTGATCCGCTCCAAGGGGGTGCGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGc                                                                                 <  4:220763/148‑1 (MQ=255)                                                          ccGGTGCTGGCGCTGAGCGTCGCTGCAGGACACCTGCGGTTCA‑ACGGCACGCCCATGGCGTTGCAGGACCGTCGTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGGGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGc                                                                                 <  8:331918/148‑1 (MQ=255)                                                                cTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGAcgc                                                                          <  5:170358/149‑1 (MQ=255)                                                                  ggCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCgg                                                                        <  4:179732/149‑1 (MQ=255)                                                                  ctcgcTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTATGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGAcgcg                                                                         <  7:265831/146‑1 (MQ=255)                                                                     gcTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCccgccgag                                                                         >  8:265831/1‑143 (MQ=255)                                                                         ttcttcGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCGGGCTGATCCGCTGCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTgg                                                                  <  3:77572/148‑1 (MQ=255)                                                                           cttcGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGcc                                                                >  8:148751/1‑148 (MQ=255)                                                                           cttcGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGGCTCTTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGcc                                                                <  6:110639/148‑1 (MQ=255)                                                                            ttcGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTc                                                                                                                                               <  5:266722/69‑1 (MQ=255)                                                                            ttcGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTc                                                                                                                                               >  6:266722/1‑69 (MQ=255)                                                                               gATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCCGtt                                                           >  8:184060/1‑149 (MQ=255)                                                                                   aaGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCTACGCGGTGCTGGCCCAGTTggg                                                        <  8:24149/148‑1 (MQ=255)                                                                                   aaGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGc                                                       <  1:145700/149‑1 (MQ=255)                                                                                       acacCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTgc                                                   >  5:281136/1‑149 (MQ=255)                                                                                        cacCTGCTGTTCA‑ACGGCACGCCCTAGGCGTTGCAGGACGGTCTTGTGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTGTATTTCGTTACCGCAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTgc                                                   <  1:89439/149‑1 (MQ=255)                                                                                           ctgctgTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTTTTTTGTAAGTGGTCCGGCTGATCCGCTCCATGGGGGTGGGGGTCTATTTCGTTACG‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCAt                                               <  1:118189/149‑1 (MQ=255)                                                                                                   cA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCgg                                                                        >  2:145751/1‑116 (MQ=255)                                                                                                   cA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCgg                                                                        <  1:145751/116‑1 (MQ=255)                                                                                                   cA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACg                                       >  1:364194/1‑149 (MQ=255)                                                                                                    a‑aCGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACgg                                      >  2:69649/1‑149 (MQ=255)                                                                                                              cccAAGGCGTTGCAGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGc                             >  4:337719/1‑149 (MQ=255)                                                                                                                     cGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACg                      <  6:281136/149‑1 (MQ=255)                                                                                                                           aGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAg                <  3:309782/149‑1 (MQ=255)                                                                                                                           aGGACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAg                <  6:132524/149‑1 (MQ=255)                                                                                                                            ggACCGTCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAgg               <  7:148751/149‑1 (MQ=255)                                                                                                                                  tCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAGGAGCAg          >  5:113922/1‑148 (MQ=255)                                                                                                                                  tCTTGAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAGGAGCAg          >  2:124435/1‑148 (MQ=255)                                                                                                                                      gAGTAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAGGAGCAGaa        <  1:69649/146‑1 (MQ=255)                                                                                                                                         cAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAGGAGCAGAAGTCGCt  >  5:23290/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                      TTGTGGCTGCTGGCCGAACTGTTCGAGCAGTTGCCCGAGCGCGGTGATGCCGACAAGCCGGTGCTGGCGCTGTTCTTCGATGAAGCACACCTGCTGTTCA‑ACGGCACGCCCAAGGCGTTGCAGGACCGTCTTGAGCAAGTGGTCCGGCTGATCCGCTCCAAGGGGGTGGGGGTCTATTTCGTTACC‑CAGTCGCCGGGTGACTTGCCCGACGCGGTGCTGGCCCAGTTGGGCTTGCGCATTCAGCACGGCTTGCGGGCGTTCACGGCCAAGGAGCAGAAGTCGCT  >  NC_002947/1562988‑1563271

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