Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	819,241:1	+C	14.5%	coding (395/627 nt)	PP_0705 →	DNA‑3‑methyladenine glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	819,241	1	.	C	14.5%	191.3 / 36.7	62	coding (395/627 nt)	PP_0705	DNA‑3‑methyladenine glycosylase
Reads supporting (aligned to +/- strand):  ref base . (26/27);  new base C (4/5);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

GCAGTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGG  >  NC_002947/819095‑819387                                                                                                                                                    |                                                                                                                                                   gCAGTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑Ga                                                                                                                                                  <  3:236644/149‑1 (MQ=255) gCAGTAGCTGGCGGTATCACAGGAAACACAGCGTGCCTGCGGCTTTTCGGCGAGCAAGCAGGTGACGGTCCAGGCCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCCACGCGGGAGCAGGCGTTAACAATGGCTGAAGAGGCGTT‑Ga                                                                                                                                                  >  2:215666/1‑149 (MQ=255)   aGTTGCTGGCGTTATCGCCGGAAACCCTGCGTGCCTGGGGCTTTTCGGCGAGCAAGCTGTTGACGGTCCAGGTCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGGGGGAGCTGGCGTTAACCATGGCTGACGAGGCGTT‑GAt                                                                                                                                                 <  1:215666/148‑1 (MQ=255)   aGTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATc                                                                                                                                                >  3:109462/1‑149 (MQ=255)            cGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCATTGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtc                                                                                                                                       <  5:112824/149‑2 (MQ=255)             gatATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcg                                                                                                                                      <  5:136024/147‑3 (MQ=255)              gTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCtggt                                                                                                                                     <  6:170360/149‑1 (MQ=255)               tATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCtggtg                                                                                                                                    <  3:5225/149‑1 (MQ=255)                   cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtca                                                                                                                                   >  6:142553/1‑144 (MQ=255)                   cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtca                                                                                                                                   <  5:142553/146‑3 (MQ=255)                   cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGgcgc                                                                                                                                >  3:194784/1‑149 (MQ=255)                    ccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtcag                                                                                                                                  >  4:340388/1‑143 (MQ=255)                    ccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtcag                                                                                                                                  <  3:340388/146‑4 (MQ=255)                      ggAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGc                                                                                                                             >  6:296827/1‑149 (MQ=255)                      ggAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGc                                                                                                                             >  6:217926/1‑149 (MQ=255)                       gAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCg                                                                                                                            >  5:100781/1‑149 (MQ=255)                         aaCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTg                                                                                                                          <  5:19795/149‑1 (MQ=255)                           cccTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTggg                                                                                                                        >  8:201332/1‑149 (MQ=255)                             cTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTgggg                                                                                                                       >  6:142281/1‑148 (MQ=255)                             cTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTgggg                                                                                                                       >  8:92753/1‑148 (MQ=255)                             cTGCGGGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGt                                                                                                                      <  7:92753/149‑1 (MQ=255)                                     cTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggt                                                                                                              <  6:40588/149‑1 (MQ=255)                                      tGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtg                                                                                                             >  3:335484/1‑149 (MQ=255)                                        cGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtgga                                                                                                           <  4:109462/149‑1 (MQ=255)                                          gCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacg                                                                                                         <  2:111218/149‑1 (MQ=255)                                          gCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacg                                                                                                         >  2:286513/1‑149 (MQ=255)                                              ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                      <  4:335484/148‑1 (MQ=255)                                              ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                      <  6:349534/148‑1 (MQ=255)                                                cGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtgga                                                                                                    <  6:100781/148‑1 (MQ=255)                                                    gAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtc                                                                                                                                    <  4:348262/112‑2 (MQ=255)                                                    gAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtc                                                                                                                                    >  3:348262/1‑111 (MQ=255)                                                           cTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGAt                                                                                        <  5:296827/149‑1 (MQ=255)                                                           cTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGAt                                                                                        <  5:217926/149‑1 (MQ=255)                                                                gACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGAt                                                                                        <  3:221671/144‑1 (MQ=255)                                                                gACGGTACAGGGCATTGCTCAAGCCAACCTGGAGTGCCTGGTACCGACGCGGGAGCAGTCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGAGTGGGGTAGGGCGGTGGAAGGTGGAGATGCTGCGGATcagc                                                                                    >  4:221671/1‑145 (MQ=255)                                                                       cAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATgctgct                                                                                           >  4:302665/1‑134 (MQ=255)                                                                       cAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATgctgct                                                                                           <  3:302665/134‑1 (MQ=255)                                                                       cAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg                                                                            >  6:270772/1‑149 (MQ=255)                                                                        aGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTgggg                                                                           >  7:43068/1‑149 (MQ=255)                                                                               tGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCt                                                                     >  5:206942/1‑148 (MQ=255)                                                                                gCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGa                                                                   >  1:322161/1‑149 (MQ=255)                                                                                      aatgcccTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCt                                                             <  6:3908/145‑1 (MQ=255)                                                                                       ccACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgc                                                                                               <  7:289947/114‑1 (MQ=255)                                                                                       ccACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgc                                                                                               >  8:289947/1‑114 (MQ=255)                                                                                           ccTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                      <  7:245095/103‑1 (MQ=255)                                                                                           ccTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                      >  8:245095/1‑103 (MQ=255)                                                                                                ctggCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGTGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGAc                                                   <  4:194784/147‑1 (MQ=255)                                                                                                    ccTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGtt                                                                       <  1:91622/125‑1 (MQ=255)                                                                                                    ccTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGtt                                                                       >  2:91622/1‑125 (MQ=255)                                                                                                      tGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGAtttt                                             >  7:165412/1‑149 (MQ=255)                                                                                                        gTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTg                                            <  5:169891/148‑1 (MQ=255)                                                                                                               cgcgGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTgcgc                                    <  1:286513/149‑1 (MQ=255)                                                                                                                    gAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAggg                               <  5:270772/149‑1 (MQ=255)                                                                                                                            gTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATcggc                       >  5:346951/1‑149 (MQ=255)                                                                                                                             ttAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATcggcg                      <  1:122338/149‑1 (MQ=255)                                                                                                                                 ccATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCt                   >  3:268891/1‑148 (MQ=255)                                                                                                                                      gCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAggg             >  3:295311/1‑149 (MQ=255)                                                                                                                                        ccaccaGGCGTTCGATC‑AACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGTCGATTTTGGCGTGAGCGAGGGGTATCGGCGGCTCAAgggg            >  5:313357/6‑148 (MQ=255)                                                                                                                                            gAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGga       <  4:268891/149‑1 (MQ=255)                                                                                                                                             aGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTg                                                            <  7:99245/95‑1 (MQ=255)                                                                                                                                             aGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTg                                                            >  8:99245/1‑95 (MQ=255)                                                                                                                                                 gTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACgg  >  1:244104/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCAGTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGG  >  NC_002947/819095‑819387

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