Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,686,122	T→A	5.7%	L177Q (CTG→CAG)	PP_3251 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,686,122	0	T	A	5.7%	221.9 / 2.2	74	L177Q (CTG→CAG)	PP_3251	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

CCACGGGCTGATTCGCACCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCTGGCCGAACTGGGCATCGACGGGCAGTATTT  >  NC_002947/3685976‑3686260                                                                                                                                                   |                                                                                                                                            ccACGGGCTGATTCGCACCTATCTGGAAGAGCTCGGCAACGGCGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTg                                                                                                                                            <  5:195418/148‑1 (MQ=255)   aCGGGCTGATTCGCACCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCaa                                                                                                                                         >  6:384701/1‑149 (MQ=255)   aCGGGCTGATTCGCACCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCa                                                                                                                                          >  3:240119/1‑148 (MQ=255)      ggCTGATTCGCACCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATTGGCTGTTTTTTTTGGGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACtt                                                                                                                                      <  1:120442/149‑1 (MQ=255)       gCTGATTCGCACCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTg                                                                                                                                     >  8:250236/1‑149 (MQ=255)              cGCACCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTgg                                                                                                                               <  5:252587/148‑1 (MQ=255)                 aCCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTgggg                                                                                                                             >  3:6930/1‑147 (MQ=255)                 aCCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTgggg                                                                                                                             >  3:179213/1‑147 (MQ=255)                       cTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑g                                                                                                                    >  2:112866/1‑149 (MQ=255)                         ggAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑Gc                                                                                                                   >  6:226634/1‑148 (MQ=255)                          gAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACt                                                                                                                      <  5:46367/145‑1 (MQ=255)                           aaGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGa                                                                                                                 >  1:172098/1‑148 (MQ=255)                              agCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGACCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACt                                                                                                             <  7:59101/149‑1 (MQ=255)                                cTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCa                                                                                                           >  4:303887/1‑149 (MQ=255)                                cTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCa                                                                                                           >  4:129319/1‑149 (MQ=255)                                          gTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCagag                                                                                                 <  5:384701/149‑1 (MQ=255)                                          gTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCagag                                                                                                 <  4:71270/149‑1 (MQ=255)                                            gATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAggt                                                                                               >  6:107033/1‑149 (MQ=255)                                             aTCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAggtg                                                                                              >  2:209577/1‑149 (MQ=255)                                             aTCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAggt                                                                                               <  5:107033/148‑1 (MQ=255)                                              tCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAggtgg                                                                                             <  3:118330/149‑1 (MQ=255)                                                cACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAggtggt                                                                                            <  4:240119/148‑1 (MQ=255)                                                cACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGtt                                                                                           >  6:151397/1‑149 (MQ=255)                                                    cTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCt                                                                                       <  8:76518/149‑1 (MQ=255)                                                       cAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTAc                                                                                     >  7:288661/1‑148 (MQ=255)                                                       cAACCATGTGCTGATATATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACCCTCCACTGCCAGAGGTGGTTGGCTAc                                                                                     >  2:212882/1‑148 (MQ=255)                                                            aTGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTc                                                                               >  5:137284/1‑149 (MQ=255)                                                             cttgCTGATCTATCAGCGCCTGCTTGCGCGCCTGGGTTTCGTAGAGGCCATGCCCGTGGATGACTCGCGGTATCTGCAAGGGGCCCTGCAAGTTGCCCTGGGGCGGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTc                                                                               <  7:306371/146‑1 (MQ=255)                                                               tgCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCTTAGAGGCCATGCCCCGGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGt                                                                            <  7:161081/149‑1 (MQ=255)                                                               tgCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGt                                                                            >  8:161081/1‑149 (MQ=255)                                                                gCTGATCTATCAGCGCCTGCTTGCTCGCCTGCGTTTCGTAGAGGCCATGCCCCTGGATGACACGCGCTATCTGCAACGGCCCCTGCAACTTGCCCTCCGGCGGCACT‑GCGACACTCACCTGCCAGAGCTGGTTGGCTACAACCTCgag                                                                            >  8:306371/1‑146 (MQ=255)                                                                   gATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCAGTTAt                                                                         >  1:128648/1‑148 (MQ=255)                                                                    aTCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTAt                                                                         <  8:189459/147‑1 (MQ=255)                                                                    aTCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTAt                                                                         <  4:179213/147‑1 (MQ=255)                                                                    aTCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATg                                                                        <  5:151397/148‑1 (MQ=255)                                                                       tATCAGCGCCTGCTTGCTCGTCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGagca                                                                    <  3:43382/149‑1 (MQ=255)                                                                        aTCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGCGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTCGTTGGCTACAACCTCGGTTATGagcag                                                                   >  7:300398/1‑149 (MQ=255)                                                                         tCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGagcagc                                                                  >  8:213621/1‑149 (MQ=255)                                                                         tCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGagca                                                                    <  3:129319/147‑1 (MQ=255)                                                                         tCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGagca                                                                    <  5:226634/147‑1 (MQ=255)                                                                         tCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGagca                                                                    <  2:172098/147‑1 (MQ=255)                                                                         tCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGag                                                                      <  3:303887/145‑1 (MQ=255)                                                                           aGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCGCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGtt                                                                                           <  4:192930/122‑1 (MQ=255)                                                                           aGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGtt                                                                                           >  3:192930/1‑122 (MQ=255)                                                                            gcgcCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCAc                                                               <  2:181796/149‑1 (MQ=255)                                                                                   cttgctCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGc                                                         <  4:6930/148‑1 (MQ=255)                                                                                   cttgctCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGc                                                         <  7:213621/148‑1 (MQ=255)                                                                                    ttgctCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCAt                                                       >  7:257779/1‑149 (MQ=255)                                                                                     tgcaagCCTGGGTTGCGTAGAGGCCATGCCCCTGGACGACACGCGGTATCTGCAAGGGGCCCTGCAACTAGCCATGGGGCACCACT‑GCGACACTCAACACCCAGAGGACGATGCCAACAACCACGGTTATGAGCAGCCACCGCTGCATc                                                      >  4:30166/6‑149 (MQ=255)                                                                                         cGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATctgct                                                  <  6:15071/149‑1 (MQ=255)                                                                                          gCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACTGGTGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGa                                                                       <  1:264182/128‑1 (MQ=255)                                                                                          gCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGa                                                                       >  2:264182/1‑127 (MQ=255)                                                                                          gCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGa                                                                       <  4:226032/127‑1 (MQ=255)                                                                                          gCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGa                                                                       >  3:226032/1‑127 (MQ=255)                                                                                          gCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATctgc                                                   <  8:257779/147‑1 (MQ=255)                                                                                              gggTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATca                                             >  3:208404/1‑149 (MQ=255)                                                                                                 ttGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATcacca                                          >  1:60183/1‑149 (MQ=255)                                                                                                       agagGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGc                                                            <  8:330956/125‑1 (MQ=255)                                                                                                       agagGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGc                                                            >  7:330956/1‑125 (MQ=255)                                                                                                       agagGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTAc                                     >  1:57235/1‑148 (MQ=255)                                                                                                        gagGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTAc                                     >  6:234545/1‑147 (MQ=255)                                                                                                         agGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCAGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTAc                                     >  2:227273/1‑146 (MQ=255)                                                                                                         agGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCAGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTAc                                     >  2:221891/1‑146 (MQ=255)                                                                                                         agGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCAGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGa                                   <  1:221891/148‑1 (MQ=255)                                                                                                         agGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCAGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGa                                   <  1:227273/148‑1 (MQ=255)                                                                                                           gCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCt                                >  7:353372/1‑149 (MQ=255)                                                                                                              aTGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCTGGc                             <  7:250236/149‑1 (MQ=255)                                                                                                                   ccTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCCACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCTGGCCGAAc                        <  1:112866/149‑1 (MQ=255)                                                                                                                       gATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCTGGCCGAACTgg                     <  5:133208/148‑1 (MQ=255)                                                                                                                       gATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCTGGCCGAACTg                      <  4:208404/147‑1 (MQ=255)                                                                                                                              cgcgGTCTCTGCACCGGGCCCCGCAACTTGCCCTGGGGCCGCACC‑GCGACACTCACCTCCCAGAGGTGGCTGGCTACAACCTCGGTTATCAGCCGCCACCGCTGCATCTGCTGAGCACCACGTACGAGCTGCCCGAACTGGCCATCGAc             <  8:288661/149‑1 (MQ=255)                                                                                                                               gcgGTATCTTCAAGGGGGCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATGTGCTGATCACCACGTACGAGCTGGCCGAACTGGGCATCGACg            <  3:30166/149‑1 (MQ=255)                                                                                                                                       tGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTCCAAATGCTGATCACAACGTACGAGCTGGCCGAACTGGGCATCGACGGCCAGTAt    >  3:156417/1‑149 (MQ=255)                                                                                                                                         cAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCTGGCCGAACTGGGCATCGACGGGCAGTAttt  >  5:214778/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                            CCACGGGCTGATTCGCACCTATCTGGAAGAGCTCGGCAACGGTGATCCACGCTGCAACCATGTGCTGATCTATCAGCGCCTGCTTGCTCGCCTGGGTTGCGTAGAGGCCATGCCCCTGGATGACACGCGGTATCTGCAAGGGGCCCTGCAACTTGCCCTGGGGCAGCACT‑GCGACACTCACCTGCCAGAGGTGGTTGGCTACAACCTCGGTTATGAGCAGCCACCGCTGCATCTGCTGATCACCACGTACGAGCTGGCCGAACTGGGCATCGACGGGCAGTATTT  >  NC_002947/3685976‑3686260

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