Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,612,097	C→A	6.6%	L346I (CTC→ATC)	PP_2287 →	internal core protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,612,097	0	C	A	6.6%	193.5 / 3.5	61	L346I (CTC→ATC)	PP_2287	internal core protein
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base A (2/2);  total (34/27)
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.44e-01

GGCTGGCTCACGGATTGCCAACGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCGGGCCTCATGCGCCCCGGTGAAGGC  >  NC_002947/2611955‑2612243                                                                                                                                               |                                                                                                                                                   ggctggctCACGGATTGCCAACGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCgggg                                                                                                                                              <  1:260268/149‑1 (MQ=255)        tCACGGATTGCCAACGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTgg                                                                                                                                       >  2:324371/1‑149 (MQ=255)               ttGCCAACGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGc                                                                                                                                >  8:79505/1‑149 (MQ=255)               ttGCCAACGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGc                                                                                                                                <  1:324371/149‑1 (MQ=255)                     aCGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCggg                                                                                                                                               >  6:129521/1‑128 (MQ=255)                     aCGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCggg                                                                                                                                               <  5:129521/128‑1 (MQ=255)                     aCGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCt                                                                                                                          >  3:268863/1‑149 (MQ=255)                              gCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGCGCCTCGGGGCCATTGGCGGTAGCCACGCTCACGTagag                                                                                                                 >  8:30933/1‑149 (MQ=255)                               cATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTag                                                                                                                   >  1:148215/1‑146 (MQ=255)                               cATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTag                                                                                                                   <  2:148215/146‑1 (MQ=255)                               cATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTa                                                                                                                    <  1:280502/145‑1 (MQ=255)                               cATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTa                                                                                                                    >  2:280502/1‑145 (MQ=255)                                        ggTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGCTGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCgctgg                                                                                                       <  5:51285/149‑1 (MQ=255)                                        ggTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCgctgg                                                                                                       >  8:43781/1‑149 (MQ=255)                                          tCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACGGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTGGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCgctggct                                                                                                     <  7:30933/149‑1 (MQ=255)                                            ggggCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCgctggct                                                                                                     >  5:107814/1‑147 (MQ=255)                                            ggggCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGa                                                                                                   <  6:107814/149‑1 (MQ=255)                                                cTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGa                                                                                                >  7:327993/1‑148 (MQ=255)                                                        aaCGTGGGCTATGACGCGCTGATCAACCCATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTACAGGGGCCGCTCGCTGAGGAAAAcgcccg                                                                                       >  5:13782/1‑149 (MQ=255)                                                               gCTATGACGCGCTGCTGAACCCATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGcc                                                                                <  6:83615/149‑1 (MQ=255)                                                                 tATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGc                                                                                                                           <  6:66616/104‑1 (MQ=255)                                                                     aCGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTgg                                                                          <  8:327993/149‑1 (MQ=255)                                                                         gcTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGgcc                                                                      >  7:111271/1‑149 (MQ=255)                                                                             atgaACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGgccgc                                                                    >  4:323476/1‑147 (MQ=255)                                                                             atgaACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCa                                                                  <  7:43781/149‑1 (MQ=255)                                                                               gaACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAg                                                                 >  8:329618/1‑148 (MQ=255)                                                                                aaCCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGa                                                               <  3:323476/149‑1 (MQ=255)                                                                                aaCCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGa                                                               <  4:268863/149‑1 (MQ=255)                                                                                 aCCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAgg                                                                                                                >  2:230183/1‑99 (MQ=255)                                                                                 aCCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAgg                                                                                                                <  1:230183/99‑1 (MQ=255)                                                                                  ccGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGa                                                               <  2:197552/147‑1 (MQ=255)                                                                                  ccGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGa                                                               >  1:197552/1‑147 (MQ=255)                                                                                   cGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTc                                                            >  3:39245/1‑149 (MQ=255)                                                                                    gATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGa                                                               >  3:57694/1‑145 (MQ=255)                                                                                    gATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGa                                                               <  4:57694/145‑1 (MQ=255)                                                                                      tACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTagag                                                                                                                 >  1:147088/1‑93 (MQ=255)                                                                                      tACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTagag                                                                                                                 <  2:147088/93‑1 (MQ=255)                                                                                           cccGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCgggg                                                    >  2:183750/1‑149 (MQ=255)                                                                                               cTCCGCACGTATGACGTCCTGTACATCTCCGGTGCGATGGGGCTGGGCCTCGGGGCGATTGGCGGTAGCGTCGCTCAGGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCGAGTGGGGCCGCCAGGAGTCGAGCGGGGCGCa                                                <  6:13782/149‑1 (MQ=255)                                                                                                tCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCa                                                >  4:138964/1‑148 (MQ=255)                                                                                                 ccGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCaaa                                              <  7:79505/149‑1 (MQ=255)                                                                                                   gCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCATCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAAt                                             >  6:91449/1‑148 (MQ=255)                                                                                                   gCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCATCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAAt                                             >  6:86352/1‑148 (MQ=255)                                                                                                    cACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATc                                            >  5:202359/1‑148 (MQ=255)                                                                                                     aCGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAg                                                                 >  8:140586/1‑126 (MQ=255)                                                                                                      cGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAg                                         <  8:111271/149‑1 (MQ=255)                                                                                                        gATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGa                                       <  1:183750/149‑1 (MQ=255)                                                                                                              gacCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTAc                                 >  1:308355/1‑149 (MQ=255)                                                                                                              gacCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTAc                                 >  8:68472/1‑149 (MQ=255)                                                                                                              gacCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTAc                                 >  2:150017/1‑149 (MQ=255)                                                                                                                ccTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTAcc                                >  2:8938/1‑148 (MQ=255)                                                                                                                    tACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTc                           >  5:195856/1‑149 (MQ=255)                                                                                                                      cATCTCCGCTGCGATGGGCCTGGGCATCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCg                          <  5:91449/148‑1 (MQ=255)                                                                                                                      cATCTCCGCTGCGATGGGCCTGGGCATCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCg                          <  5:86352/148‑1 (MQ=255)                                                                                                                       aTCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCggg                        >  6:220233/1‑149 (MQ=255)                                                                                                                          tcCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCaa                                          >  4:256888/1‑128 (MQ=255)                                                                                                                          tcCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCaa                                          <  3:256888/128‑1 (MQ=255)                                                                                                                            cGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAg                                         <  1:200880/127‑1 (MQ=255)                                                                                                                            cGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAg                                         >  2:200880/1‑127 (MQ=255)                                                                                                                             gCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCGGGCCTCAt                  >  5:323127/1‑149 (MQ=255)                                                                                                                                  gATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCGGGCCTCATGCGcc             <  1:8938/149‑1 (MQ=255)                                                                                                                                       gcctgggcctCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAg                                         >  1:207100/1‑116 (MQ=255)                                                                                                                                       gcctgggcctCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAg                                         <  2:207100/116‑1 (MQ=255)                                                                                                                                         ctgggcctCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCGGGCCTCATGCGCCCCGGTGa      <  3:138964/149‑1 (MQ=255)                                                                                                                                               ctCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCGGGCCTCATGCGCCCCGGTGAAgga  >  5:185005/1‑146 (MQ=255)                                                                                                                                               ctCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCGGGCCTCATGCGCCCCGGTGAAg    <  6:185005/145‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   GGCTGGCTCACGGATTGCCAACGCAGTCCGCATGGGCGTGGTCGGGGCTGGGACCAACGTGGGCTATGACGCGCTGATGAACCGATACAACCCGCTCCGCACGGATGACGACCTGTACATCTCCGCTGCGATGGGCCTGGGCCTCGGGGCCATTGGCGGTAGCCTCGCTCACGTAGAGGGGCCGCTGGCTGAGGAAAACGCCCGCCTGGCCCAGTGGGGCCGCCAGGAGTCCAGCGGGGCGCAAATCAAGGAGCTTACCGATTCGGGCCTCATGCGCCCCGGTGAAGGC  >  NC_002947/2611955‑2612243

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