Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,170,187	C→T	6.8%	A139V (GCC→GTC)	PP_4548 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,170,187	0	C	T	6.8%	172.3 / 2.0	59	A139V (GCC→GTC)	PP_4548	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base T (2/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

GCCATGCGCTCGCAGATGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGTGCGTTGAAGAACCT  >  NC_002947/5170039‑5170334                                                                                                                                                     |                                                                                                                                                       gCCATGCGCTCGCAGATGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGc                                                                                                                                                        >  1:361166/1‑149 (MQ=255)     tGCGCTCGCAGATGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGGCCTCAAGCGTGGCCACCTGCATGCAGCCAt                                                                                                                                                     >  1:387669/1‑148 (MQ=255)            gCAGATGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCATGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTGAAGCGTGGGCACCTGCATTCAGCCATGAAGCCt                                                                                                                                              <  4:65360/148‑1 (MQ=255)               gATGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGc                                                                                                                                          <  4:180088/149‑1 (MQ=255)                aTGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGca                                                                                                                                         >  3:193102/1‑149 (MQ=255)                 tGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGcac                                                                                                                                        <  1:122556/149‑1 (MQ=255)                 tGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGcac                                                                                                                                        <  4:37286/149‑1 (MQ=255)                  gcgcAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCCCCTGCATGCAGCCATGAAGCCTTCGcaca                                                                                                                                       >  2:39639/1‑149 (MQ=255)                   cgcAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCTACCTGCATGCAGCCATGAAGCCTTCGCACAt                                                                                                                                      <  2:372588/149‑1 (MQ=255)                       aCCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATgaac                                                                                                                                  <  5:95764/149‑1 (MQ=255)                        ccATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATgaac                                                                                                                                  >  1:361480/1‑148 (MQ=255)                               ggCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACtgc                                                                                                                           <  4:335921/148‑1 (MQ=255)                                  gcCGAAGTTGACGACTTCATATGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACtgcgt‑‑‑gc                                                                                                                    >  6:209715/1‑149 (MQ=255)                                     gAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGtt                                                                                                                 >  1:204261/1‑149 (MQ=255)                                      aaGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTc                                                                                                                <  1:115493/149‑1 (MQ=255)                                         ttGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCATCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGAGTTGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAg                                                                                                             >  5:205469/1‑149 (MQ=255)                                         gtGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAg                                                                                                             <  7:273236/148‑1 (MQ=255)                                          tGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAgg                                                                                                            >  8:83079/1‑149 (MQ=255)                                            acgacTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGc                                                                                                           <  1:88493/148‑1 (MQ=255)                                                 ttCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAgc                                                                                                     >  2:324776/1‑149 (MQ=255)                                                      cTGGCACTTGCGCTGGCGCGGGCACGAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTAGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGcgcg                                                                                                  <  4:255373/147‑1 (MQ=255)                                                      cTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTAGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCgcgga                                                                                                >  3:255373/1‑148 (MQ=255)                                                           aattGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GAGTTCTCGGCCGAAGAGCGGCGCCGt                                                                                           >  7:361773/3‑149 (MQ=255)                                                           aCTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGt                                                                                           >  3:160001/1‑149 (MQ=255)                                                            cTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTg                                                                                          <  2:361166/149‑1 (MQ=255)                                                            cTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTg                                                                                          <  1:377939/149‑1 (MQ=255)                                                                 gcTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCAtt                                                                                     <  5:209715/149‑2 (MQ=255)                                                                  cTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATgg                                                                                    >  7:42574/1‑149 (MQ=255)                                                                   tGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATggg                                                                                   >  8:233578/1‑149 (MQ=255)                                                                     gcgcGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGc                                                                                  >  6:125116/1‑148 (MQ=255)                                                                           gCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGCAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAgg                                                                           <  5:275371/149‑1 (MQ=255)                                                                            cACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGt                                                                          <  4:193102/149‑1 (MQ=255)                                                                                aGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCaa                                                                      <  7:233578/149‑1 (MQ=255)                                                                                 ggTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAAt                                                                     >  4:350798/1‑149 (MQ=255)                                                                                     aTCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGGCCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCt                                                                 <  2:387669/149‑1 (MQ=255)                                                                                      tCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTc                                                                <  7:225123/149‑1 (MQ=255)                                                                                       cGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTc                                                                >  7:23059/1‑148 (MQ=255)                                                                                                    gcgcgcTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAg                                                   <  3:352200/148‑1 (MQ=255)                                                                                                     cgcgcTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTg                                                 <  7:83079/149‑1 (MQ=255)                                                                                                              gCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGTCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCAt                                        >  1:14719/1‑149 (MQ=255)                                                                                                              gCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCAt                                        <  3:263638/149‑1 (MQ=255)                                                                                                                aTCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCt                                      <  8:355017/149‑1 (MQ=255)                                                                                                                aTCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCt                                      >  5:252183/1‑149 (MQ=255)                                                                                                                     ttGTGACCTCAAGCGTGGCCACCTGCATGCAGTCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCt                                                                 >  5:162525/1‑117 (MQ=255)                                                                                                                     ttGTGACCTCAAGCGTGGCCACCTGCATGCAGTCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCt                                                                 <  6:162525/117‑1 (MQ=255)                                                                                                                                aGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGACGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTAt                      <  4:235299/149‑1 (MQ=255)                                                                                                                                aGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTg                                                                         <  8:173017/98‑1 (MQ=255)                                                                                                                                aGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTg                                                                         >  7:173017/1‑98 (MQ=255)                                                                                                                                 gCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCg                            <  4:362413/142‑1 (MQ=255)                                                                                                                                 gCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCg                            >  3:362413/1‑142 (MQ=255)                                                                                                                                 gCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATc                     >  3:49892/1‑149 (MQ=255)                                                                                                                                    tGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTgg                  >  6:66821/1‑149 (MQ=255)                                                                                                                                      gCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGt                <  7:77219/149‑1 (MQ=255)                                                                                                                                          ccTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGTTCCCCCATCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGTg               <  1:324776/149‑1 (MQ=255)                                                                                                                                              catgcaGCCATGAAGCCTTCGCACCTGAACGAACTGCGT‑‑‑GCGGTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGGGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGTGCGTTgaa        <  2:361480/149‑1 (MQ=255)                                                                                                                                               atgcaGTCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGTGCGTTgaag       <  2:14719/149‑1 (MQ=255)                                                                                                                                                 gcaGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGt                >  3:129911/1‑138 (MQ=255)                                                                                                                                                 gcaGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGt                <  4:129911/138‑1 (MQ=255)                                                                                                                                                   aGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGTTGCGTTGAAGAAcc   >  3:41/1‑149 (MQ=255)                                                                                                                                                    gCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGTGCGTTGAAGAACCt  <  4:170082/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GCCATGCGCTCGCAGATGCGCAACCATCTGGGCGCCGAAGTTGACGACTTCATCTGGCACTTGCGCTGGCGCGGGCACCAGGTGATCGAGCAGCGCGTGGCGCGCTATGGCATCGATTGTGACCTCAAGCGTGGCCACCTGCATGCAGCCATGAAGCCTTCGCACATGAACGAACTGCGT‑‑‑GCGTTCGAGGCCGAAGCGCGGCGCCGTGGCATGGGCGATCAGGTGCAATTGCTCGATCGCGCTGCAGTGGCCGAGCATCTGCAAAGCCCGCTTTATCTGGGTGCGTTGAAGAACCT  >  NC_002947/5170039‑5170334

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