Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,088,803	T→C	6.7%	W16R (TGG→CGG)	PP_0943 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,088,803	0	T	C	6.7%	215.8 / 4.4	74	W16R (TGG→CGG)	PP_0943	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base C (2/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCCTGCTGCTGGA  >  NC_002947/1088657‑1088948                                                                                                                                                    |                                                                                                                                                  gCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTgg                                                                                                                                                 <  4:51522/149‑1 (MQ=255)   ccTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCAAGTAACCCGGCGATGATGTTCATGAAAT‑ACGTCCTCCGCGGACTGCCCTAACTGGAAAACTGGCGCCGGCt                                                                                                                                               >  4:138370/1‑149 (MQ=255)         gcagggTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATAATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCg                                                                                                                                         >  2:100382/4‑149 (MQ=255)           cGGGTTTTTATTTCATATCGTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGATTTTCATTATCCAGTAAGCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGgc                                                                                                                                       <  8:130342/149‑1 (MQ=255)              ggttttATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCCGGCTGAGCCAgcgcc                                                                                                                                    <  3:138370/147‑2 (MQ=255)               tttttATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCAt                                                                                                                                   >  5:153386/1‑149 (MQ=255)                ttttATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCCGGCTGAGCCAgcgcca                                                                                                                                   <  4:7529/148‑3 (MQ=255)                 tttATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTCGAAAACTGGCGCCGGCTGAGCCAgcgccag                                                                                                                                  <  1:26618/148‑4 (MQ=255)                         atgaCCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTgcg                                                                                                                         <  1:150335/146‑1 (MQ=255)                          tcaCCTTTTTCTTGAAATGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGAGCTGGCTCAGCCGGCGCATTCGCTGTgcgc                                                                                                                        >  6:209798/1‑149 (MQ=255)                                       gAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCcgac                                                                                                             <  2:59807/147‑1 (MQ=255)                                       gAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCcga                                                                                                              <  2:72700/146‑1 (MQ=255)                                        aaGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCcgacga                                                                                                           <  2:83050/148‑1 (MQ=255)                                         aGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGc                                                                                                         <  6:63090/149‑1 (MQ=255)                                          gTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGcc                                                                                                        >  7:231471/1‑149 (MQ=255)                                            gATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCcg                                                                                                       <  8:231471/148‑1 (MQ=255)                                             aTTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGc                                                                                                                               <  5:302827/123‑1 (MQ=255)                                             aTTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGc                                                                                                                               >  6:302827/1‑123 (MQ=255)                                             aTTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCcgcg                                                                                                     >  3:166501/1‑149 (MQ=255)                                              ttCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCgcgc                                                                                                    <  6:153386/149‑1 (MQ=255)                                                   tttATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCATCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGAt                                                                                               >  8:201733/1‑149 (MQ=255)                                                         aCTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAg                                                                                           <  1:100196/147‑1 (MQ=255)                                                         aCTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAg                                                                                           >  2:100196/1‑147 (MQ=255)                                                         aCTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCa                                                                                         >  2:150484/1‑149 (MQ=255)                                                         aCTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCa                                                                                         <  1:31636/149‑1 (MQ=255)                                                          cTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCAt                                                                                        >  6:233767/1‑149 (MQ=255)                                                              ataatGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCa                                                                                                                                             <  6:268071/92‑1 (MQ=255)                                                              ataatGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCa                                                                                                                                             >  5:268071/1‑92 (MQ=255)                                                                     atatCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTgcg                                                                                                                         >  7:81678/1‑105 (MQ=255)                                                                     atatCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTgcg                                                                                                                         <  8:81678/105‑1 (MQ=255)                                                                        tCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAgg                                                                          <  1:77411/149‑1 (MQ=255)                                                                         cATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGa                                                                                            >  3:93933/1‑130 (MQ=255)                                                                         cATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGa                                                                                            <  4:93933/130‑1 (MQ=255)                                                                          aTTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGAGCTGGCTCAGACGGCTCAATCGCTGTGAGCTCGACCACGACGAGCCGCGACTGATAGAGCATTACATGGCAGAAGGcc                                                                        >  5:88371/1‑149 (MQ=255)                                                                           ttATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTTCGCGTACTGCCCTACCTGGAAATCTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGATCCGCGCCTGATCGAGCATTAACTGGCCGAAGGCCg                                                                       <  1:100382/149‑1 (MQ=255)                                                                           ttATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCg                                                                       <  5:233767/149‑1 (MQ=255)                                                                                 aGTAACCCGGCGATGATGTTCATGAAATGCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATc                                                                  >  1:17296/1‑149 (MQ=255)                                                                                   tAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACATGGCCGAAGGCCGCTATCTgg                                                               >  6:170061/1‑149 (MQ=255)                                                                                     aCCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGgt                                                              <  4:166501/148‑1 (MQ=255)                                                                                      cccGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGgtgt                                                            <  3:239444/149‑1 (MQ=255)                                                                                        cGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGtgc                                                          >  6:61909/1‑149 (MQ=255)                                                                                          gCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGtgctg                                                        <  1:150484/149‑1 (MQ=255)                                                                                             atgatgTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCAc                                                     >  8:21973/1‑149 (MQ=255)                                                                                                  gTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGaaa                                                >  3:40860/1‑149 (MQ=255)                                                                                                  gTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGaaa                                                >  7:11737/1‑149 (MQ=255)                                                                                                     cATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAAcc                                              <  7:201733/148‑1 (MQ=255)                                                                                                     cATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCt                                             >  8:130788/1‑149 (MQ=255)                                                                                                      aTGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTc                                            >  6:151405/1‑149 (MQ=255)                                                                                                       tGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCg                                           <  7:21973/149‑1 (MQ=255)                                                                                                       tGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCACATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCg                                           >  5:48745/1‑149 (MQ=255)                                                                                                        gAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGc                                          >  5:261416/1‑149 (MQ=255)                                                                                                         aaaT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGc                                          <  6:261416/148‑1 (MQ=255)                                                                                                                    tccGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTg                               >  5:155213/1‑149 (MQ=255)                                                                                                                    tccGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTg                               >  2:181085/1‑149 (MQ=255)                                                                                                                       gcgAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGgcg                            <  6:155213/149‑1 (MQ=255)                                                                                                                        cgAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGgcgc                           <  1:181085/149‑1 (MQ=255)                                                                                                                               cccTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAg                     <  8:10981/148‑1 (MQ=255)                                                                                                                               cccTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGc                    >  6:313043/1‑149 (MQ=255)                                                                                                                                 cTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGa                                                                                            <  2:206432/75‑1 (MQ=255)                                                                                                                                 cTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGa                                                                                            >  1:206432/1‑75 (MQ=255)                                                                                                                                    ccTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGtt                 <  2:17296/147‑1 (MQ=255)                                                                                                                                    ccTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTc               >  5:38841/1‑149 (MQ=255)                                                                                                                                    ccTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTc               <  5:313043/149‑1 (MQ=255)                                                                                                                                      tGGAAAACTGGCGCTGGCTCAGCCGGCGAATTCGCTGTGCGCTCGACCCCGACGAGCCGAGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGc             >  6:141091/1‑149 (MQ=255)                                                                                                                                        gAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGcc            <  7:130788/148‑1 (MQ=255)                                                                                                                                        gAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGcc            <  4:40860/148‑1 (MQ=255)                                                                                                                                         agcgaatgcgcCGGCTGAGCCAGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCctg          >  1:273397/8‑149 (MQ=255)                                                                                                                                          aaaCTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCctgc         <  5:141091/149‑1 (MQ=255)                                                                                                                                              tGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCctgctgct     >  2:115344/1‑149 (MQ=255)                                                                                                                                               ggCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCctgctgctg    >  4:244392/1‑149 (MQ=255)                                                                                                                                                gcgcTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCctgctgct     <  5:151405/147‑1 (MQ=255)                                                                                                                                                 cgcTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCCTGCTGCTGGa  <  3:244392/149‑1 (MQ=255)                                                                                                                                                   cTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGaaa                                                >  1:195270/1‑101 (MQ=255)                                                                                                                                                   cTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGaaa                                                <  2:195270/101‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAT‑CCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGCTGGCTCAGCCGGCGCATTCGCTGTGCGCTCGACCCCGACGAGCCGCGCCTGATCGAGCATTACCTGGCCGAAGGCCGCTATCTGGTGTGCTGCACCGAAACCTCGCCATGGACGGTGGCGCTGACAGCGTTTCGCCTGCTGCTGGA  >  NC_002947/1088657‑1088948

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