Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,386,817	T→C	6.4%	D237G (GAT→GGT)	PP_3863 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,386,817	0	T	C	6.4%	182.2 / 2.2	62	D237G (GAT→GGT)	PP_3863	tail protein
Reads supporting (aligned to +/- strand):  ref base T (25/33);  new base C (2/2);  total (27/35)
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.91e-01

CCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTGCCCAGCTCCAGC  >  NC_002947/4386672‑4386963                                                                                                                                                   |                                                                                                                                                   ccgccACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATcc                                                                                                                                                  <  4:27621/148‑1 (MQ=255)  cgccACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCCCCCGCGCGCGCGCCATGGCTGGGGTCCTATGACCCGCCGGGTTGATCACCTTGCCCCGCCGTCGAGCCACGC‑GGTCATCCGc                                                                                                                                                >  7:307421/1‑149 (MQ=255)                cGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGT‑GGTCATCCGCCAGTTTGGCCTCGa                                                                                                                                  <  7:234204/149‑1 (MQ=255)                    tAGTTCAGGGGCATGGCCCGACTGATCGCGTGGGCCCGCTCCCACTCCCCCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTc                                                                                                                              <  6:48569/149‑1 (MQ=255)                         cAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTg                                                                                                                          <  8:238816/148‑1 (MQ=255)                                   gCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGc                                                                                                                <  6:73002/148‑1 (MQ=255)                                            aTCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGgccgc                                                                                                                     <  6:339205/134‑1 (MQ=255)                                            aTCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGgccgc                                                                                                                     >  5:339205/1‑134 (MQ=255)                                                cgTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCATTGGCCGCCGTGCCGAACATCGTGTCg                                                                                                  <  3:197954/149‑1 (MQ=255)                                                cgTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCg                                                                                                  >  5:348198/1‑149 (MQ=255)                                                cgTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCg                                                                                                  <  5:145818/149‑1 (MQ=255)                                                  ttGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGt                                                                                                <  3:94858/149‑1 (MQ=255)                                                  ttGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGt                                                                                                >  4:94858/1‑149 (MQ=255)                                                     gCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGc                                                                                              <  7:94454/148‑1 (MQ=255)                                                       ccGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgc                                                                                           >  4:37051/1‑149 (MQ=255)                                                        cGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgct                                                                                          >  5:118592/1‑149 (MQ=255)                                                           tCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgctgcg                                                                                       >  1:178235/1‑149 (MQ=255)                                                            cccACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCACCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgctgcgc                                                                                      >  7:166200/1‑149 (MQ=255)                                                                      ccGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCACCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgctgc                                                                                        >  3:288954/1‑137 (MQ=255)                                                                      ccGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCACCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgctgc                                                                                        <  4:288954/137‑1 (MQ=255)                                                                        gcgcgcgcgcCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCt                                                                           >  4:124525/1‑148 (MQ=255)                                                                           ctcgggcCATGGCTGGGTTCATCTGACCCGACGGGTTGATCACCTTGCCCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATg                                                                       <  8:307421/144‑1 (MQ=255)                                                                              gcgcCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATc                                                                                                        <  7:118009/113‑1 (MQ=255)                                                                              gcgcCATGGCTGGGGTCATCTGACCCGACTCGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATc                                                                                                        >  8:118009/1‑113 (MQ=255)                                                                                gcCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTc                                                                  <  4:185041/149‑1 (MQ=255)                                                                                    ccagtcgggTCAGATGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGa                                                               >  2:203805/7‑148 (MQ=255)                                                                                      gCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGaaaa                                                            <  8:38697/149‑1 (MQ=255)                                                                                            gTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGtt                                                      <  3:37051/149‑1 (MQ=255)                                                                                               aTCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGTCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGga                                                   <  7:106306/149‑1 (MQ=255)                                                                                                tCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCAGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTg                                                     >  1:263997/1‑147 (MQ=255)                                                                                                tCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGc                                                                                                                <  7:117482/87‑1 (MQ=255)                                                                                                tCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGc                                                                                                                >  8:117482/1‑87 (MQ=255)                                                                                                  tGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGgaga                                                 >  3:252974/1‑148 (MQ=255)                                                                                                    aCCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAAt                                              >  7:275594/1‑149 (MQ=255)                                                                                                      ccGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTc                                                                  >  5:235345/1‑127 (MQ=255)                                                                                                      ccGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTc                                                                  <  6:235345/127‑1 (MQ=255)                                                                                                             gTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTgg                                     >  4:85855/1‑149 (MQ=255)                                                                                                             gTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTgg                                     >  5:283436/1‑149 (MQ=255)                                                                                                             gTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTg                                      <  3:79512/148‑1 (MQ=255)                                                                                                             gTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTg                                      <  1:203805/148‑1 (MQ=255)                                                                                                             gTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTg                                      <  8:275594/148‑1 (MQ=255)                                                                                                              ttGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGgt                                    >  1:85493/1‑149 (MQ=255)                                                                                                              ttGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGgt                                    <  2:178235/149‑1 (MQ=255)                                                                                                              ttGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGgt                                    <  2:85493/149‑1 (MQ=255)                                                                                                              ttGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGgt                                    <  4:252974/149‑1 (MQ=255)                                                                                                                gATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGgtgt                                  >  6:145677/1‑149 (MQ=255)                                                                                                                  tCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATc                                             <  4:305495/136‑1 (MQ=255)                                                                                                                  tCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATc                                             >  3:305495/1‑136 (MQ=255)                                                                                                                    aCCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCAc                               <  6:335487/148‑1 (MQ=255)                                                                                                                      cTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTg                                                                          >  8:321596/1‑103 (MQ=255)                                                                                                                      cTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTg                                                                          <  7:321596/103‑1 (MQ=255)                                                                                                                      cTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGa                                                                         <  6:191159/104‑1 (MQ=255)                                                                                                                      cTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGa                                                                         >  5:191159/1‑104 (MQ=255)                                                                                                                      cTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCg                            >  3:39972/1‑149 (MQ=255)                                                                                                                           cccGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAgg                       <  4:39972/149‑1 (MQ=255)                                                                                                                                    gAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTg              <  6:118592/149‑1 (MQ=255)                                                                                                                                    gAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTg              <  6:283436/149‑1 (MQ=255)                                                                                                                                     aGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTGc             >  3:125509/1‑149 (MQ=255)                                                                                                                                       ccACGC‑GGTCACCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTGccc           <  8:166200/149‑1 (MQ=255)                                                                                                                                              ggTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTGCCCAGCTcc     <  2:263997/149‑1 (MQ=255)                                                                                                                                                tCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCACTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTGCCCAGCTCCAg   >  4:308301/1‑149 (MQ=255)                                                                                                                                                  aTCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTGCCCAGCTCCAGc  <  6:348198/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGC‑GGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATCCAGATTGGTGTCACCCGTCAGGATGTTCTTGCCCAGCTCCAGC  >  NC_002947/4386672‑4386963

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