Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,650,019	A→T	8.4%	intergenic (+1496/+195)	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,019	0	A	T	8.4%	179.7 / 2.6	63	intergenic (+1496/+195)	PP_3214/PP_3216	hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base T (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGGC  >  NC_002947/3649877‑3650166                                                                                                                                                |                                                                                                                                                    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                                                                                           <  8:369106/132‑1 (MQ=255)                                                                      aagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                           >  7:369106/1‑132 (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                                                          <  7:288137/148‑1 (MQ=255)                                                                                       gCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATgg                                                          <  5:133764/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                                              <  3:203777/149‑1 (MQ=255)                                                                                                  gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGaaa                                              <  1:19238/149‑1 (MQ=255)                                                                                                  gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCAGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTAt                                                                                <  8:88784/115‑1 (MQ=255)                                                                                                  gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCAGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTAt                                                                                >  7:88784/1‑115 (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                                            >  1:358575/1‑132 (MQ=255)                                                                                                                     tGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAgg                                            <  2:358575/132‑1 (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                   <  1:193654/144‑1 (MQ=255)                                                                                                                                  tCACCGTCTA‑CGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg                   >  2:193654/1‑144 (MQ=255)                                                                                                                                    aCCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTc                          <  6:109361/135‑1 (MQ=255)                                                                                                                                    aCCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTc                          >  5:109361/1‑135 (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                                   >  2:47446/1‑122 (MQ=255)                                                                                                                                        tCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCttcttc                                   <  1:47446/122‑1 (MQ=255)                                                                                                                                           t‑CGTGGCCGGTGAGTGCGGGGGCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCtt     <  8:357467/149‑1 (MQ=255)                                                                                                                                               tAGAAGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTAAGTCGATCAGGGCCAGCTCCTCCGTGGc  >  6:321845/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    CTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCAT‑CGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGGC  >  NC_002947/3649877‑3650166

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