Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	19,117	C→T	6.3%	A637V (GCA→GTA)	PP_0014 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	19,117	0	C	T	6.3%	189.7 / 2.3	63	A637V (GCA→GTA)	PP_0014	transposase
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base T (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CCCGCCGTTCGCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCTCAAGTCGAGCGAGCGTTGAC  >  NC_002947/18970‑19263                                                                                                                                                     |                                                                                                                                                   cccGCCGTTCGCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCcagca                                                                                                                                                   <  1:155266/149‑1 (MQ=255)  ccgccgTTCGCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCGCGTGACTCTATTCCTAACATTCCGCcagca                                                                                                                                                   <  1:28132/149‑1 (MQ=255)  ccgccgTTCGCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGAACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCcagcag                                                                                                                                                  >  1:310645/1‑149 (MQ=255)     ccgTTCGCAGACGCTCAGAAAACTACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACAGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCt                                                                                                                                               <  5:4500/149‑1 (MQ=255)          cGCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAg                                                                                                                                          >  4:42801/1‑149 (MQ=255)           gCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCagcagc                                                                                                                                                 <  7:172055/141‑1 (MQ=255)           gCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGTAGCCTTCAAgg                                                                                                                                         >  8:261175/1‑149 (MQ=255)           gCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAgg                                                                                                                                         >  1:339235/1‑149 (MQ=255)           gCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATCGACCG‑GTGACTCTATTCCTAACATTCCGCCagcagc                                                                                                                                                 >  8:172055/1‑141 (MQ=255)            cAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGt                                                                                                                                        >  1:197904/1‑149 (MQ=255)                  cTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGt                                                                                                                                  >  1:204694/1‑149 (MQ=255)                     aGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGtt                                                                                                                                 <  2:197904/147‑1 (MQ=255)                     aGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTtc                                                                                                                                <  2:71950/148‑1 (MQ=255)                     aGAAAACGACGCGCAGCTGCAATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTtc                                                                                                                                <  6:254979/148‑1 (MQ=255)                       aaaaCGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTtcgtc                                                                                                                             >  1:324906/1‑149 (MQ=255)                          acgacgCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGTAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGcc                                                                                                                          <  4:313294/149‑1 (MQ=255)                             acgCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGTAGCCTTCa                                                                                                                                            <  2:168005/128‑1 (MQ=255)                             acgCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGTAGCCTTCa                                                                                                                                            >  1:168005/1‑128 (MQ=255)                                 gcAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTgg                                                                                                                   <  6:41011/149‑1 (MQ=255)                                          aTAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCCCCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTAGTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCGGGCTGGGAGTTCt                                                                                                            >  4:285629/1‑147 (MQ=255)                                           tAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAg                                                                                                         >  2:143258/1‑149 (MQ=255)                                            aGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCcagca                                                                                                                                                   >  4:197392/1‑106 (MQ=255)                                            aGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCcagca                                                                                                                                                   <  3:197392/106‑1 (MQ=255)                                               gAAAATCTACAAGTAGTGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTTACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTTCTGGGAGTTc                                                                                                             <  2:368794/141‑1 (MQ=255)                                               gAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCAAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTc                                                                                                             >  1:368794/1‑141 (MQ=255)                                                aaaaTCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATg                                                                                                    <  3:42801/149‑1 (MQ=255)                                                   aTCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAAc                                                                                                 >  6:26499/1‑149 (MQ=255)                                                          aggaggGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCa                                                                                          <  3:265321/149‑1 (MQ=255)                                                           ggaggGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCa                                                                                          >  4:12045/1‑148 (MQ=255)                                                            gaggGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAAt                                                                                        <  2:339235/149‑1 (MQ=255)                                                            gaggGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAAt                                                                                        >  2:343219/1‑149 (MQ=255)                                                              gggTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCt                                                                                      <  6:263417/149‑1 (MQ=255)                                                               ggTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTGTAAGCCATGAACGCCATCAATCTg                                                                                     >  6:191774/1‑149 (MQ=255)                                                                  cGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTAc                                                                                  <  1:343219/149‑1 (MQ=255)                                                                      tGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAAt                                                                              <  2:324906/149‑1 (MQ=255)                                                                       gACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAAt                                                                              <  2:310645/148‑1 (MQ=255)                                                                        aCCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGc                                                                            >  3:348499/1‑149 (MQ=255)                                                                           gAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTAc                                                                         >  2:90439/1‑149 (MQ=255)                                                                             aaGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGa                                                                       <  2:204694/149‑1 (MQ=255)                                                                                   gCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATgg                                                                 <  5:215005/149‑1 (MQ=255)                                                                                     cTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGt                                                                >  7:1814/1‑148 (MQ=255)                                                                                      tCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTg                                                               >  6:214028/1‑148 (MQ=255)                                                                                        cACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGaa                                                             <  8:1814/148‑1 (MQ=255)                                                                                        cACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAAt                                                            <  5:214028/149‑1 (MQ=255)                                                                                         aCCCAATGTACTCGCTGAAGTGGAATGGACCG‑TTGACTCTATTCCTAAAATTCCTCCAGCAGCCTTCAAGGTCTAGGTTCGTCGCCCTGCTGGGAGTTCTAATCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATa                                                           <  5:360150/149‑1 (MQ=255)                                                                                          cccAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATa                                                           >  3:12668/1‑148 (MQ=255)                                                                                          cccAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATAc                                                          <  5:26499/149‑1 (MQ=255)                                                                                           ccAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGAGAATACa                                                         >  4:348597/1‑149 (MQ=255)                                                                                                  aCTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGCACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTAc                                                  >  3:192817/1‑149 (MQ=255)                                                                                                   cTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACg                                                 <  4:348499/149‑1 (MQ=255)                                                                                                   cTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACg                                                 >  4:305361/1‑149 (MQ=255)                                                                                                    tCGCTGAAGTGGAATGGACCG‑GTGAGTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCctgctg                                                                                                                    <  6:42569/81‑1 (MQ=255)                                                                                                    tCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCctgctg                                                                                                                    >  5:42569/1‑81 (MQ=255)                                                                                                      gCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGTGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTc                                              >  8:291769/1‑149 (MQ=255)                                                                                                               gAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGt                                     >  3:232367/1‑149 (MQ=255)                                                                                                                    gACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTa                                 <  3:12045/148‑1 (MQ=255)                                                                                                                    gACCG‑GTGAATCTATTCCTAACATTCCGCCAGTAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTa                                 <  7:261175/148‑1 (MQ=255)                                                                                                                     aCCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTAtt                               <  1:143258/149‑1 (MQ=255)                                                                                                                          gTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTcaca                           <  1:90439/149‑1 (MQ=255)                                                                                                                              ctctATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTAACAGcc                        >  8:189226/1‑148 (MQ=255)                                                                                                                               tctATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGccc                       >  1:80372/1‑148 (MQ=255)                                                                                                                               tctATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCt                      >  1:250369/1‑149 (MQ=255)                                                                                                                               tctATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCt                      >  6:75128/1‑149 (MQ=255)                                                                                                                                   ttCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCTCAAg                  <  2:80372/149‑1 (MQ=255)                                                                                                                                      cTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCTCAAGTc                >  5:290336/1‑148 (MQ=255)                                                                                                                                       tAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGACATGAAAGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATAGTTTTATTCACAGCCATCAAGTc                >  7:309544/1‑147 (MQ=255)                                                                                                                                        aaCATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAg                          >  6:273698/1‑136 (MQ=255)                                                                                                                                        aaCATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAg                          <  5:273698/136‑1 (MQ=255)                                                                                                                                         aCATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCTCAAGTcgagc            >  4:227696/1‑149 (MQ=255)                                                                                                                                            ttCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTAACAGCCCTCAAGTcgagcgag         >  3:21484/1‑149 (MQ=255)                                                                                                                                                  cagcagCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCTCAAGTCGAGCGAGCGTTGa   >  8:322034/1‑149 (MQ=255)                                                                                                                                                   agcagcCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCTCAAGTCGAGCGAGCGTTGAc  <  3:227696/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCCGCCGTTCGCAGACGCTCAGAAAACGACGCGCAGCTGCCATAGCGAAAATCGACAAGGAGGGTCGAGTGACCGAAAGCACGCCTCCACCCAAGGTACTCGCTGAAGTGGAATGGACCG‑GTGACTCTATTCCTAACATTCCGCCAGCAGCCTTCAAGGTCTCGGTTCGTCGCCCTGCTGGGAGTTCTAAGCCATGAACGCCATCAATCTGTACGAATGCTACGAGCATGGTGAATACAGCTTTACGCTCCGCGATCGTTTTATTCACAGCCCTCAAGTCGAGCGAGCGTTGAC  >  NC_002947/18970‑19263

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