Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	715,959	C→T	7.0%	Q60* (CAG→TAG)	PP_0609 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	715,959	0	C	T	7.0%	165.3 / 2.4	57	Q60* (CAG→TAG)	PP_0609	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base T (2/2);  total (28/29)
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.59e-01

GCTGATGTTGGCGCTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGCCTGATCGTGGCGGAGTTGCCTGGGCAAG  >  NC_002947/715817‑716106                                                                                                                                                |                                                                                                                                                      gCTGATGTTGGCGCTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGata                                                                                                                                                 >  5:549/1‑149 (MQ=255)       gTTGGCGCTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCa                                                                                                                                           >  4:191606/1‑149 (MQ=255)            cgcTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCa                                                                                                                                                      >  3:305278/1‑133 (MQ=255)            cgcTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCa                                                                                                                                                      <  4:305278/133‑1 (MQ=255)             gcTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCTAGGATATTCGCATGGCCg                                                                                                                                     >  1:341368/1‑149 (MQ=255)             gcTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCTAGGATATTCGCATGGCCg                                                                                                                                     >  1:346699/1‑149 (MQ=255)                gACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCa                                                                                                                                  <  2:116311/149‑1 (MQ=255)                 aCTATCGGGTTGATGCTGGTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGGGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAt                                                                                                                                 <  3:8308/149‑1 (MQ=255)                   tATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGGCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑Gt                                                                                                                             >  6:238054/1‑149 (MQ=255)                   tATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑Gt                                                                                                                             >  6:219630/1‑149 (MQ=255)                     tCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑Gtt                                                                                                                            <  4:158242/148‑1 (MQ=255)                     tCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑Gt                                                                                                                             <  5:219630/147‑1 (MQ=255)                       gggTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTgg                                                                                                                         >  4:177325/1‑149 (MQ=255)                       gggTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTgg                                                                                                                         <  6:549/149‑1 (MQ=255)                        ggTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGc                                                                                                                        >  7:82105/1‑149 (MQ=255)                        ggTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGc                                                                                                                        >  7:159791/1‑149 (MQ=255)                              tGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCGCCGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGttt                                                                                                                   <  1:271876/149‑1 (MQ=255)                                      gcggcCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCTAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGaa                                                                                                          <  2:346699/149‑1 (MQ=255)                                      gcggcCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCTAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGaa                                                                                                          <  2:341368/149‑1 (MQ=255)                                        ggcCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAAc                                                                                                        >  5:286402/1‑149 (MQ=255)                                                 tCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGAtt                                                                                               <  3:191606/149‑1 (MQ=255)                                                     ctctTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGa                                                                                           >  7:55235/1‑149 (MQ=255)                                                           gCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTcc                                                                                                   >  8:316719/1‑135 (MQ=255)                                                           gCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTcc                                                                                                   <  7:316719/135‑1 (MQ=255)                                                            cATCGGCTCAACAGGCCTGGCGGCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCg                                                                                    <  5:238054/149‑1 (MQ=255)                                                                ggCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCa                                                                                 <  8:159791/148‑1 (MQ=255)                                                                ggCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCAt                                                                                >  3:194674/1‑149 (MQ=255)                                                                agCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCAt                                                                                <  1:144982/148‑1 (MQ=255)                                                                          gCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGt                                                                      <  4:28534/149‑1 (MQ=255)                                                                               gcgcCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTtcg                                                                 >  4:266929/1‑149 (MQ=255)                                                                                cgcCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTtcgc                                                                <  8:55235/149‑1 (MQ=255)                                                                                 gcCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGACG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTtcgct                                                               >  1:224989/1‑149 (MQ=255)                                                                                    tGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTtcgctcgc                                                            >  2:335165/1‑149 (MQ=255)                                                                                     gCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGAGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGcc                                                           >  5:28769/1‑149 (MQ=255)                                                                                      cAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTtcgctcgc                                                            <  3:177325/147‑1 (MQ=255)                                                                                      cAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGcc                                                           <  6:28769/148‑1 (MQ=255)                                                                                             ctgcgctgcgCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGt                                                    >  3:287443/1‑148 (MQ=255)                                                                                             ctgcgctgcgCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTg                                                   <  2:224989/149‑1 (MQ=255)                                                                                                cgctgcgcATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGAt                                                 <  4:287443/148‑1 (MQ=255)                                                                                                cgctgcgcATGCAGGATGAGGAG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGAATGGCCGGGAT‑‑GTTTGGCTGTGTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATc                                                >  2:37910/1‑149 (MQ=255)                                                                                                   tgcgcATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGAGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGaa                                             >  3:160440/1‑149 (MQ=255)                                                                                                     cgcATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGc                                           <  6:286402/149‑1 (MQ=255)                                                                                                      gcATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCATAGGTTTGGGTTTTTGCGCGTGAAACGTACCTTTTTCGATGCCCCGGCATAACAGAGGGTGTTGGCGCGCCCCGGGGTGATGGacg                                            >  1:276125/1‑147 (MQ=255)                                                                                                       cATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGa                                                                                           >  6:184551/1‑99 (MQ=255)                                                                                                       cATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGa                                                                                           <  5:184551/99‑1 (MQ=255)                                                                                                         tGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTc                                        >  8:44584/1‑148 (MQ=255)                                                                                                               aTGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTc                                 <  4:160440/149‑1 (MQ=255)                                                                                                                  aGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGc                              <  3:266929/149‑1 (MQ=255)                                                                                                                    gcg‑cGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAgg                                                                       <  4:309034/106‑1 (MQ=255)                                                                                                                    gcg‑cGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAgg                                                                       >  3:309034/1‑106 (MQ=255)                                                                                                                      g‑cgGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCa                                      <  3:98860/137‑1 (MQ=255)                                                                                                                      g‑cgGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCa                                      >  4:98860/1‑137 (MQ=255)                                                                                                                                    ctgcGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGCCTGATCGTGGCGGAGtt             <  2:276125/149‑1 (MQ=255)                                                                                                                                       cgcATGGACCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGCCTGATCGTGGCGGAGTTGcc          <  4:194674/149‑1 (MQ=255)                                                                                                                                         cATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGt                                                    <  8:90093/105‑1 (MQ=255)                                                                                                                                         cATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGt                                                    >  7:90093/1‑105 (MQ=255)                                                                                                                                             gCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGCCTGATCGTGGCGGAGTTGCCTGGGc     >  5:256388/1‑148 (MQ=255)                                                                                                                                              cccAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGc                              <  2:208693/122‑1 (MQ=255)                                                                                                                                              cccAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGc                              >  1:208693/1‑122 (MQ=255)                                                                                                                                               ccAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGCCTGATCGTGGCGGAGTTGCCTGGGCAAg  >  6:353101/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      GCTGATGTTGGCGCTGACTATCGGGTTGATGCTGTTGGCGGCCGCCAGTCAGCTCTTTGCATCGGCTCAACAGGCCTGGCGCCTGCAAAGTACTGCGCTGCGCATGCAGGATGAGGCG‑CGGCAGGCGTTGCTGCGCATGGCCCAGGATATTCGCATGGCCGGCAT‑‑GTTTGGCTGTTTGCGCCTGAAACCTTCCGATTTCGATGACCCGGCATTACAGCAGGTCTTCGCTCGCCCCCTGGTGATCGAAGCGTCCACCCTCAGCCTGATCGTGGCGGAGTTGCCTGGGCAAG  >  NC_002947/715817‑716106

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