Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,508,373	A→C	7.0%	L105R (CTT→CGT)	PP_2202 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,508,373	0	A	C	7.0%	176.3 / 4.2	57	L105R (CTT→CGT)	PP_2202	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base C (2/2);  total (30/27)
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.09e-01

CACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCAGCAGC  >  NC_002947/2508229‑2508519                                                                                                                                                 |                                                                                                                                                     cACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCgg                                                                                                                                                  <  1:328598/149‑1 (MQ=255)  aCGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCggg                                                                                                                                                 <  5:225693/149‑1 (MQ=255)   cgcgAAACCTGAAAGGTGTTTCCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGTGTCACCATCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGt                                                                                                                                                <  5:27573/149‑1 (MQ=255)    gcgAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTc                                                                                                                                               <  1:8884/149‑1 (MQ=255)      gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTaa                                                                                                                                                      <  1:202711/140‑1 (MQ=255)      gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTaa                                                                                                                                                      >  2:202711/1‑140 (MQ=255)      gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGc                                                                                                                                             >  2:249086/1‑149 (MQ=255)           cTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATc                                                                                                                                        >  4:247811/1‑149 (MQ=255)              aaaGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGTATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCAc                                                                                                                                      <  1:85226/148‑1 (MQ=255)                  gtgtTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGc                                                                                                                                  >  8:232551/1‑148 (MQ=255)                       aCCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAg                                                                                                                            >  1:104823/1‑149 (MQ=255)                       aCCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAg                                                                                                                            >  1:286836/1‑149 (MQ=255)                         cAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCgg                                                                                                                                          >  5:156643/1‑132 (MQ=255)                         cAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCgg                                                                                                                                          <  6:156643/133‑2 (MQ=255)                             cGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGagc                                                                                                                      >  4:212015/1‑149 (MQ=255)                              gTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGCCATCCAGCTGagcc                                                                                                                     >  6:351523/1‑148 (MQ=255)                                   aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGAt                                                                                                                                         <  7:163384/124‑1 (MQ=255)                                   aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGAt                                                                                                                                         >  8:163384/1‑124 (MQ=255)                                   aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCg                                                                                                                                           <  2:223963/122‑1 (MQ=255)                                   aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCg                                                                                                                                           >  1:223963/1‑122 (MQ=255)                                      aGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GAc                                                                                                            <  2:319869/149‑1 (MQ=255)                                       gCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GAc                                                                                                            >  8:4353/1‑148 (MQ=255)                                                      cGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGtt                                                                                            <  3:212015/149‑1 (MQ=255)                                                       gCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTg                                                                                           >  2:79863/1‑149 (MQ=255)                                                           gCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAg                                                                                       >  3:121517/1‑149 (MQ=255)                                                                    aGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCtgctgc                                                                              >  7:177944/1‑149 (MQ=255)                                                                            gCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCAcgc                                                                        <  7:81815/147‑1 (MQ=255)                                                                            gCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCAcgc                                                                        >  8:81815/1‑147 (MQ=255)                                                                              cTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCa                                                                    <  2:286836/149‑1 (MQ=255)                                                                                   ccAGGGTCACCCGCAACCGGAACCGCCAGAATGGCTCAACCATTAACCCGGCGATCACCGTAAGCGGGGCGCCGATCACCGGCATCCAGCTGACCATCATCAGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTg                                                                <  1:249086/149‑1 (MQ=255)                                                                                         tCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAgcg                                                         >  8:17633/1‑149 (MQ=255)                                                                                           aCCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCt                                                       >  7:154719/1‑149 (MQ=255)                                                                                               gCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑Gaaa                                                  <  1:79863/149‑1 (MQ=255)                                                                                                    aGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGaa                                             >  7:271211/1‑149 (MQ=255)                                                                                                          cgcCAGAATGGCTAACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAgtg                                       >  4:249445/1‑147 (MQ=255)                                                                                                           gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCgctg                                                                                  >  3:198746/1‑106 (MQ=255)                                                                                                           gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCgctg                                                                                  <  4:198746/106‑1 (MQ=255)                                                                                                           gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAgcgc                                      >  4:38614/1‑149 (MQ=255)                                                                                                            ccAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAg                                                                                       >  8:202512/1‑100 (MQ=255)                                                                                                             cAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCtt                                    >  5:127870/1‑149 (MQ=255)                                                                                                                 aTGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTcg                                <  8:154719/149‑1 (MQ=255)                                                                                                                    gCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCgcgc                             <  8:271211/149‑1 (MQ=255)                                                                                                                        aCGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGTGCTTGTCGCGCAg                           <  3:249445/147‑1 (MQ=255)                                                                                                                        aCGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAg                           <  8:177944/147‑1 (MQ=255)                                                                                                                                tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGAt                  <  4:121517/148‑1 (MQ=255)                                                                                                                                tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATg                 >  1:225881/1‑149 (MQ=255)                                                                                                                                tCCCGGCGATCAGCGTAAGCGGGGCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATg                 <  7:232551/149‑1 (MQ=255)                                                                                                                                   cGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcc              >  6:135487/1‑149 (MQ=255)                                                                                                                                   cGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTAGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcc              >  6:50883/1‑149 (MQ=255)                                                                                                                                     gCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcccg            <  5:135487/149‑1 (MQ=255)                                                                                                                                         tCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAg                           <  2:3835/130‑1 (MQ=255)                                                                                                                                         tCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAg                           >  1:3835/1‑130 (MQ=255)                                                                                                                                           ctcGTAATCGGGGCGCCGATCACCCGCATCCAGCCGAGCAGCACC‑GACCATTGCCCCCACCGTTGGTAGCGCTGCTGGGCACGCGCCAGTTGGGCATCGCTCGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCca       <  1:76461/147‑1 (MQ=255)                                                                                                                                           aGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTCGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCcag      <  5:50883/149‑1 (MQ=255)                                                                                                                                            gCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCcagc     <  2:104823/149‑1 (MQ=255)                                                                                                                                            gCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCcag      >  2:57456/1‑148 (MQ=255)                                                                                                                                            gCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGATGCGCACGCGCCAGTTGGGCAGCGCT‑GAATGGGAACCAGCGCTTGTCGCGCAGGTGATCGATggcccggcccgg      >  5:205081/1‑146 (MQ=255)                                                                                                                                             cGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCcagca    <  7:4353/149‑1 (MQ=255)                                                                                                                                             cGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCcag      >  7:290650/1‑147 (MQ=255)                                                                                                                                                 aGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCagcagc  >  8:222458/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGC‑GACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCT‑GAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCAGCAGC  >  NC_002947/2508229‑2508519

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