Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,403,012	A→T	8.1%	F17Y (TTC→TAC)	PP_2105 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,403,012	0	A	T	8.1%	170.9 / 2.0	58	F17Y (TTC→TAC)	PP_2105	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base T (2/3);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTACCGAGTCCTGGGTCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCAAGTGCGGCTGCTGGTATGGGCCTTG  >  NC_002947/2402866‑2403145                                                                                                                                                   |                                                                                                                                      cACGTACCGAGTCCTGGGTCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGaa                                                                                                                                      <  1:31187/148‑1 (MQ=255)    gTACCGAGTCCTGGGTCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCa                                                                                                                                  <  1:245416/149‑1 (MQ=255)               tGGGTCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAata                                                                                                                        >  4:135641/1‑148 (MQ=255)                gggTCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGaa                                                                                                                                      <  4:127730/133‑1 (MQ=255)                gggTCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGaa                                                                                                                                      >  3:127730/1‑133 (MQ=255)                   tCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAat                                                                                                                         <  1:68818/143‑1 (MQ=255)                   tCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAat                                                                                                                         >  2:68818/1‑143 (MQ=255)                         gCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCa                                                                                                             >  1:333694/1‑149 (MQ=255)                          cAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCAt                                                                                                            >  2:7576/1‑149 (MQ=255)                             ggTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCt                                                                                                         <  7:400153/149‑1 (MQ=255)                             ggTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCAGCAAAATCTCTTGTGCATTCt                                                                                                         >  1:79909/1‑149 (MQ=255)                               tttCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCttt                                                                                                       >  7:140675/1‑149 (MQ=255)                                 tCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTg                                                                                                      <  4:139690/148‑1 (MQ=255)                                 tCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGACTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGc                                                                                                     >  6:76442/1‑149 (MQ=255)                                  cgcgGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGTAGTCAGCCATCAATATCTCTTGTGCATTCttt                                                                                                       <  6:309866/146‑1 (MQ=255)                                  cgcgGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGcc                                                                                                    >  8:408861/1‑149 (MQ=255)                                  cgcgGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGGGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGTAGTCAGCCATCAATATCTCTTGTGCATTCttt                                                                                                       >  5:309866/1‑146 (MQ=255)                                        ccTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACAAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTg                                                                                              >  6:174228/1‑149 (MQ=255)                                           gggCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTg                                                                                           >  4:225959/1‑149 (MQ=255)                                            ggCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTTGTCGCCGGCCACGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTgg                                                                                          <  1:357041/149‑1 (MQ=255)                                            ggCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTgg                                                                                          <  1:351577/149‑1 (MQ=255)                                             gCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGc                                                                                         >  7:402577/1‑149 (MQ=255)                                                  gcgcTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATgg                                                                                    <  4:36173/149‑1 (MQ=255)                                                    gcTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCATGTTGGCGGTGTCGGCCATCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCATCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGt                                                                                  <  7:272845/149‑1 (MQ=255)                                                       cGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCAt                                                                                                                             <  6:96479/103‑1 (MQ=255)                                                       cGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCAt                                                                                                                             >  5:96479/1‑103 (MQ=255)                                                        gTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGtctct                                                                              <  5:76442/149‑1 (MQ=255)                                                             tCTGATCGCGCAATTAGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCAAAAAGGCATGGAATTAAGAAATCAATATATCTTGTGGATTGTTTTACGATTTTCTGTCCATGGGtctctgcgtg                                                                         >  4:396626/1‑144 (MQ=255)                                                               tGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAATGCCTGGAAGTCAGCCCTCAATATCGCTTGTGCATTGTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTc                                                                       >  8:228243/1‑149 (MQ=255)                                                                gCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTc                                                                       >  2:160682/1‑148 (MQ=255)                                                                 cTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGt                                                                        <  5:414940/146‑1 (MQ=255)                                                                    gcgcAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGAt                                                                    <  3:131493/147‑1 (MQ=255)                                                                    gcgcAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGa                                                                  <  1:7576/149‑1 (MQ=255)                                                                    gcgcAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGa                                                                  <  3:315352/149‑1 (MQ=255)                                                                     cgcAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAg                                                                 >  7:108754/1‑149 (MQ=255)                                                                      gcAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAgt                                                                <  3:225959/149‑1 (MQ=255)                                                                      gcAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAgt                                                                <  6:160863/149‑1 (MQ=255)                                                                         aTTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGtt                                                             <  3:135641/149‑1 (MQ=255)                                                                                 gCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGTAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTc                                                     >  6:401031/1‑149 (MQ=255)                                                                                    tGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCTCCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAg                                                  >  3:50475/1‑149 (MQ=255)                                                                                      gTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGt                                                 >  4:268280/1‑148 (MQ=255)                                                                                      gTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGt                                                 <  8:140675/148‑1 (MQ=255)                                                                                      gTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGt                                                 <  3:268280/148‑1 (MQ=255)                                                                                      gTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGt                                                 <  6:395465/148‑1 (MQ=255)                                                                                      gTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGGCGCGCACCAAGGCCTGGTAGTCAGCCATCAATATCTCTTGTGAATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGt                                                 <  3:393116/148‑1 (MQ=255)                                                                                       tCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGa                                               <  8:108754/149‑1 (MQ=255)                                                                                       tCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTAGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGa                                               <  2:333694/149‑1 (MQ=255)                                                                                         gCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACa                                             <  8:402577/149‑1 (MQ=255)                                                                                                  gTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATg                                    >  5:117211/1‑149 (MQ=255)                                                                                                       cGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAg                                                  >  8:15224/1‑130 (MQ=255)                                                                                                       cGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAg                                                  <  7:15224/130‑1 (MQ=255)                                                                                                       cGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGAACTTCGAGTGACAGCACGCATGGAAgg                               >  8:329142/1‑149 (MQ=255)                                                                                                        ggTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGAGTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGt                              >  7:198349/1‑149 (MQ=255)                                                                                                            tCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCaa                          <  1:160682/149‑1 (MQ=255)                                                                                                                          ctcGGTGTCGCGCACCAAGGCATGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCAAGTGCGGctgctg              <  7:408861/147‑1 (MQ=255)                                                                                                                              gtgtCTCGCAACAAGGACTGGAAGTCATCCATCATTATCTGTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCAAGTGCGGCTGCTGGTATgg        <  3:364399/149‑1 (MQ=255)                                                                                                                              gtgtCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCAAGTGCGGCTGCTGGTATgg        <  6:117211/149‑1 (MQ=255)                                                                                                                                   gcgcACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTACCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCAAGTGCGGCTGCTGGTATGGGCCtt   >  1:319196/1‑149 (MQ=255)                                                                                                                                    cgcAGCAAGGCCTGGTAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCAAGTGCGGCTGCTGGTATGGGCCTTg  <  5:401031/149‑1 (MQ=255)                                                                                                                                         cAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACg                                        >  1:54535/1‑106 (MQ=255)                                                                                                                                         cAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACg                                        <  2:54535/106‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                      CACGTACCGAGTCCTGGGTCAGTTGCAGGGTTTCGCGGGCCTGGGCCAGGCGCTCGTGAATCTGCTCGCGCAATTCGTCGGCCTGGTCGCCGGCCAGGTTGGCGGTGTCGGCCAGCAGCTTCTCGGTGTCGCGCACCAAGGCCTGGAAGTCAGCCATCAATATCTCTTGTGCATTCTTTGCCGATTTGCTGGCCATGGGTCTCTCCCTGTCGATGAGTGTTGGTACTTCGAGTGACAGCACGCATGGAAGGTTCAAGTGCGGCTGCTGGTATGGGCCTTG  >  NC_002947/2402866‑2403145

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