Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,655,512	G→C	6.1%	L419L (CTC→CTG) ‡	PP_1449 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,655,512	0	G	C	6.1%	210.3 / 3.1	65	L419L (CTC→CTG) ‡	PP_1449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base C (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GGTTGTCGCTTGCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTG  >  NC_002947/1655364‑1655659                                                                                                                                                      |                                                                                                                                                    ggTTGTCGCTTGCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCg                                                                                                                                                     >  2:182360/1‑149 (MQ=255)  gTTGTCGCTTGCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGa                                                                                                                                                    <  7:171247/149‑1 (MQ=255)            gCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTcc                                                                                                                                            <  8:316966/147‑1 (MQ=255)              cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  6:228850/1‑149 (MQ=255)              cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  6:354502/1‑149 (MQ=255)              cTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGc                                                                                                                                        >  7:137367/1‑149 (MQ=255)               tttATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGcc                                                                                                                                       >  7:180302/1‑149 (MQ=255)               tttATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGcc                                                                                                                                       >  7:349325/1‑149 (MQ=255)                      aCATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGTTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCg                                                                                                                                <  7:59472/149‑1 (MQ=255)                        aTGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGc                                                                                                                                    >  2:227171/1‑143 (MQ=255)                        aTGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGc                                                                                                                                    <  1:227171/143‑1 (MQ=255)                            gCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTc                                                                                                                           >  4:322799/1‑148 (MQ=255)                             cccGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGa                                                                                                                         <  3:117242/149‑1 (MQ=255)                              ccGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGAc                                                                                                                        >  4:354048/1‑149 (MQ=255)                               cGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGAcc                                                                                                                       <  5:228850/149‑1 (MQ=255)                                  cgcTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGGCGCTATCGCTGTCACCGCGCCATACATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTc                                                                                                                    <  4:324106/149‑1 (MQ=255)                                    cTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAg                                                                                                                  >  5:243934/1‑149 (MQ=255)                                       aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCt                                                                                                                <  8:180302/148‑1 (MQ=255)                                       aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTc                                                                                                               <  3:322799/149‑1 (MQ=255)                                       aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTc                                                                                                               <  5:277381/149‑1 (MQ=255)                                       aGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCGGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGCCCGGCACCGCTGTCGACCGTCAGCTc                                                                                                               <  4:291008/149‑1 (MQ=255)                                             cTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTAc                                                                                                          <  5:354502/148‑1 (MQ=255)                                                   tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGtt                                                                                                            >  5:245960/1‑140 (MQ=255)                                                   tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGtt                                                                                                            <  6:245960/140‑1 (MQ=255)                                                   tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGCATCgaca                                                                                                                             <  7:344062/123‑4 (MQ=38)                                                   tttGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCCTGGACATCCTCGCCGGCATCgaca                                                                                                                             >  8:344062/1‑120 (MQ=38)                                                        ggTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCt                                                                                              >  2:40278/1‑149 (MQ=255)                                                          tttCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTc                                                                                            >  3:169326/1‑149 (MQ=255)                                                          tttCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGACGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTc                                                                                            >  8:137849/1‑149 (MQ=255)                                                              attGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCCCCTTCgcgg                                                                                        >  4:252368/2‑149 (MQ=255)                                                               ttGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCgcggc                                                                                       >  4:293891/1‑149 (MQ=255)                                                               ttGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCgcggc                                                                                       >  2:224805/1‑149 (MQ=255)                                                                       gTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATgc                                                                               <  1:224805/149‑1 (MQ=255)                                                                          gCTATCGCTGGCACCGCGCCATAATTC‑CTTGCGATGACGCTCGTGCTGTTTGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATgcgc                                                                             <  5:152533/148‑1 (MQ=255)                                                                           cTATCGCTGTCACCGCGCCATAAATCGCTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCa                                                                            >  1:192264/1‑149 (MQ=255)                                                                           cTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCAt                                                                           <  5:48062/149‑1 (MQ=255)                                                                              tCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATg                                                                          <  3:293891/147‑1 (MQ=255)                                                                                 cTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATgtcgtc                                                                     <  5:153512/149‑1 (MQ=255)                                                                                 cTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATgtcgtc                                                                     >  4:161201/1‑149 (MQ=255)                                                                                  tGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGtcgtcg                                                                    >  6:220458/1‑149 (MQ=255)                                                                                      aCCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCg                                                                 <  8:137367/148‑1 (MQ=255)                                                                                      aCCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCg                                                                 >  6:311776/1‑148 (MQ=255)                                                                                       ccGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGa                                                               <  3:46938/149‑1 (MQ=255)                                                                                       ccGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGa                                                               <  6:243934/149‑1 (MQ=255)                                                                                        cgcgCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGAc                                                              >  6:258667/1‑149 (MQ=255)                                                                                             cATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCaa                                                         <  3:354048/149‑1 (MQ=255)                                                                                               tAAATC‑CTTGCGATGACGCACCTGCTGTGCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGTGTCCACCTTCGCGGCGGTCATGCGCATGGCGTCGCCGGACTGCAACa                                                       <  1:267468/149‑1 (MQ=255)                                                                                                 aaTC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGc                                                     <  3:161201/149‑1 (MQ=255)                                                                                                          cGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGctct                                              <  5:311776/148‑1 (MQ=255)                                                                                                          cGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCtctc                                             <  8:349325/149‑1 (MQ=255)                                                                                                           gATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCa                                            <  7:204204/149‑1 (MQ=255)                                                                                                                  cACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCa                                     <  2:81645/149‑1 (MQ=255)                                                                                                                    ccTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAg                                    >  3:242865/1‑148 (MQ=255)                                                                                                                       gctgTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAAt                                >  6:374428/1‑149 (MQ=255)                                                                                                                       gctgTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAAt                                >  2:34999/1‑149 (MQ=255)                                                                                                                                  cGTTTGCTGTCGATGCCGACGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCtg                     <  7:137849/149‑1 (MQ=255)                                                                                                                                        cTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCTGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGttgttg               <  5:258667/149‑1 (MQ=255)                                                                                                                                        cTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGttgttg               <  1:182360/149‑1 (MQ=255)                                                                                                                                         tGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGc              >  4:331939/1‑149 (MQ=255)                                                                                                                                             gATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgt          <  4:169326/149‑1 (MQ=255)                                                                                                                                              aTGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgtg         <  4:242865/149‑1 (MQ=255)                                                                                                                                              aTGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgt          <  7:57594/148‑1 (MQ=255)                                                                                                                                              aTGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTtgt          >  8:57594/1‑148 (MQ=255)                                                                                                                                                  cGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTgtgggtg     >  6:351647/1‑149 (MQ=255)                                                                                                                                                   ggCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTgtgggtgg    >  2:315235/1‑149 (MQ=255)                                                                                                                                                     cGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGggtggtg  <  5:374428/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GGTTGTCGCTTGCCTTTATCGACATGTGCCCGGCGCTCAGGGTACTGCCCTTTGCGGTTTCCTTGAGCTGGTCGCTATCGCTGTCACCGCGCCATAAATC‑CTTGCGATGACGCACCTGCTGTTCGGTGCGTTTGCTGTCGATGCCGGCGAGGATGTCCAGGCCGGCACCGCTGTCGACCGTCAGCTCGTTACCGGCGTCCACCTTCGCGGCGGTCATGCGCATGTCGTCGCCGGACTGCAACAGCACGCTCTCAGTGCCCAGCAATTGACTGCCCTGTTGTTGCTTGTGGGTGGTG  >  NC_002947/1655364‑1655659

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