Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,650,015	T→A	6.3%	intergenic (+1492/+199)	PP_3214 → / ← PP_3216	hypothetical protein/cysteine dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,650,015	0	T	A	6.3%	185.7 / 2.9	63	intergenic (+1492/+199)	PP_3214/PP_3216	hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base A (2/2);  total (31/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 = 8.64e-01

AAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTT  >  NC_002947/3649868‑3650163                                                                                                                                                    |                                                                                                                                                      aaaGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑c                                                                                                                                                     <  5:3075/149‑1 (MQ=255)   aGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑Cg                                                                                                                                                    >  2:168161/1‑148 (MQ=255)          cTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACgg                                                                                                                                             <  3:149243/148‑1 (MQ=255)           ttGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGCCGGt                                                                                                                                            >  6:2210/1‑148 (MQ=255)           ttGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGACCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGt                                                                                                                                            >  3:64086/1‑148 (MQ=255)                aCGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCTA‑CGATGACGGTGAgggg                                                                                                                                      <  3:341905/149‑4 (MQ=255)                     gACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGgtg                                                                                                                                  <  3:5656/148‑1 (MQ=255)                      aCATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAACCGCTCCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTc                                                                                                                                >  3:77326/1‑149 (MQ=255)                      aCATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTc                                                                                                                                >  6:30243/1‑149 (MQ=255)                       cATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTcc                                                                                                                               >  7:346511/1‑149 (MQ=255)                           aCGAAACAACTCGCCGCCCGCAGCGCGCCACCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCccgcc                                                                                                                           >  4:203777/1‑149 (MQ=255)                             gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCCGTAGACGGTGAGTGCGGTGTCccgccg                                                                                                                          <  1:23947/149‑1 (MQ=255)                             gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                         >  3:70878/1‑149 (MQ=255)                              aaaCAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCggg                                                                                                                        <  5:9565/149‑1 (MQ=255)                                     tCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGAGGTGTCCCGCCGGGGAAcgcg                                                                                                                 >  4:361791/1‑149 (MQ=255)                                          gcccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGa                                                                                                            <  3:370587/149‑1 (MQ=255)                                             cgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGCCTGCGGCGAGCCCCTCACCGTCAT‑CGTGGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCg                                                                                                             >  7:357467/1‑145 (MQ=255)                                                 gcgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTc                                                                                                     <  2:261406/149‑1 (MQ=255)                                                  cgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTcc                                                                                                    <  1:243807/149‑1 (MQ=255)                                                         gCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGc                                                                                                              >  4:141933/1‑132 (MQ=255)                                                         gCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGc                                                                                                              <  3:141933/132‑1 (MQ=255)                                                          ccGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCg                                                                                                                                     >  1:110430/1‑108 (MQ=255)                                                          ccGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCg                                                                                                                                     <  2:110430/108‑1 (MQ=255)                                                          ccGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAg                                                                                            >  2:19238/1‑149 (MQ=255)                                                            ggCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGtt                                                                                          >  1:230241/1‑149 (MQ=255)                                                              cAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                        <  3:279086/149‑1 (MQ=255)                                                                   ccgccATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAg                                                                                   <  1:168161/149‑1 (MQ=255)                                                                     gccATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAggg                                                                                 <  4:70878/149‑1 (MQ=255)                                                                         taaagcaagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGGGCGGGGACTGCGGCTAGCCCCTCACCGTCAT‑CTTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTAt                                                                             <  1:318851/147‑1 (MQ=255)                                                                               aagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                        >  7:369106/1‑132 (MQ=255)                                                                               aagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                        <  8:369106/132‑1 (MQ=255)                                                                               aagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGc                                                                       >  6:133764/1‑149 (MQ=255)                                                                                   cAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCcgc                                                                    >  3:7261/1‑148 (MQ=255)                                                                                      cATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGaa                                                                 >  2:68960/1‑148 (MQ=255)                                                                                             gTAGCCCAATGCCAGTAAAAGCTGCGCTGTGACTGCGGCGAGCCCCTAACCGTCAT‑CGTTGACGGTGAGTGCGGTGTACCGCCGGGGAACGCGGGCGAAGATTTAAACCTCCAGTATCGACAGGCGAAAGTTGGCCGCGAACAGAGGAt                                                         >  6:92033/1‑149 (MQ=255)                                                                                               aGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATg                                                        <  8:346511/148‑1 (MQ=255)                                                                                                gCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATgg                                                       <  5:133764/148‑1 (MQ=255)                                                                                                gCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATgg                                                       <  7:288137/148‑1 (MQ=255)                                                                                                  ccAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGcc                                                    >  6:25916/1‑149 (MQ=255)                                                                                                         cAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCg                                              >  8:302444/1‑148 (MQ=255)                                                                                                          aGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGa                                             >  3:92402/1‑148 (MQ=255)                                                                                                           gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGaaa                                           <  1:19238/149‑1 (MQ=255)                                                                                                           gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGaaa                                           <  3:203777/149‑1 (MQ=255)                                                                                                           gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCAGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTAt                                                                             >  7:88784/1‑115 (MQ=255)                                                                                                           gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCAGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTAt                                                                             <  8:88784/115‑1 (MQ=255)                                                                                                                cGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGTGtc                                      <  4:64086/149‑1 (MQ=255)                                                                                                                        ggcgACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTcgc                              <  3:361791/149‑1 (MQ=255)                                                                                                                           gACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa                           >  7:141511/1‑149 (MQ=255)                                                                                                                              tGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAgg                                         <  2:358575/132‑1 (MQ=255)                                                                                                                              tGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAgg                                         >  1:358575/1‑132 (MQ=255)                                                                                                                                 tgcgAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGGCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGGTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGa                     <  4:77326/148‑1 (MQ=255)                                                                                                                                   caccgcacTCACCGTCTA‑CGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTAGGTCGAt                    >  7:328939/8‑148 (MQ=255)                                                                                                                                      gCCCCTCACCGTCAT‑CGTAGAGGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGTGTATGTTTGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg                <  5:2210/149‑1 (MQ=255)                                                                                                                                      gCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg                <  7:302444/149‑1 (MQ=255)                                                                                                                                           tCACCGTCTA‑CGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg                >  2:193654/1‑144 (MQ=255)                                                                                                                                           tCACCGTCTA‑CGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg                <  1:193654/144‑1 (MQ=255)                                                                                                                                             aCCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTc                       >  5:109361/1‑135 (MQ=255)                                                                                                                                             aCCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTc                       <  6:109361/135‑1 (MQ=255)                                                                                                                                               cGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATg                                                        <  7:51920/100‑1 (MQ=255)                                                                                                                                               cGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATg                                                        >  8:51920/1‑100 (MQ=255)                                                                                                                                               cGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGctcc       >  1:82484/1‑149 (MQ=255)                                                                                                                                                 tCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCttcttc                                <  1:47446/122‑1 (MQ=255)                                                                                                                                                 tCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCttcttc                                >  2:47446/1‑122 (MQ=255)                                                                                                                                                    t‑CGTGGCCGGTGAGTGCGGGGGCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCtt  <  8:357467/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTT  >  NC_002947/3649868‑3650163

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