Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,812,702	A→G	6.6%	D328G (GAC→GGC)	PP_3368 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,812,702	0	A	G	6.6%	247.8 / 4.4	85	D328G (GAC→GGC)	PP_3368	MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (41/38);  new base G (3/3);  total (44/41)
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.70e-01

GGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCG  >  NC_002947/3812568‑3812842                                                                                                                                       |                                                                                                                                             ggTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCAATGTCCATGAGgc                                                                                                                               <  2:1249368/150‑1 (MQ=255) ggTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGgc                                                                                                                               >  2:1322925/1‑150 (MQ=255)      tGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGTGGCGCCTa                                                                                                                          >  1:1732190/1‑150 (MQ=255)      tGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTa                                                                                                                          >  1:767336/1‑150 (MQ=255)      tGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTa                                                                                                                          >  1:1488326/1‑150 (MQ=255)        aTGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGc                                                                                                                        <  1:446453/150‑1 (MQ=255)          gtgGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCAGACCAGGAGCTGTGGAGCCTGCTGCTAGGCATCGGCCTGGTGGGCACCATGAGCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCt                                                                                                                      >  2:216143/1‑150 (MQ=255)           tgGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTg                                                                                                                     <  2:709156/150‑1 (MQ=255)             gtcggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGGGGAACCGCTACAGATCCGCCCAAAGGATTTTGCCAACGAGGCGCCAAGCCTggg                                                                                                                   >  2:713624/1‑149 (MQ=255)             gtcggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGc                                                                                                                   >  1:740049/1‑150 (MQ=255)                ggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTgg                                                                                                                    <  1:562340/146‑1 (MQ=255)                ggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTgg                                                                                                                    >  2:562340/1‑146 (MQ=255)                         gctgctGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCaa                                                                                                       <  2:1303405/150‑1 (MQ=255)                          ctgctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                      <  1:102347/150‑1 (MQ=255)                          ctgctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCAACTTCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                      >  1:1482570/1‑150 (MQ=255)                           tgctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAcc                                                                                                     >  1:1604514/1‑150 (MQ=255)                             ctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCAt                                                                                                   >  2:933792/1‑150 (MQ=255)                              tgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATg                                                                                                  <  2:303279/150‑1 (MQ=255)                               gTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgc                                                                                                 >  2:349741/1‑150 (MQ=255)                                ttCTTTCCCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcg                                                                                                >  2:500178/1‑150 (MQ=255)                                  cTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcggc                                                                                              >  2:1220268/1‑150 (MQ=255)                                    ttANCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGt                                                                                            >  2:1199118/1‑150 (MQ=255)                                    tgtccgtCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGt                                                                                            <  2:1644614/144‑1 (MQ=255)                                     tACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGtt                                                                                           <  1:349741/150‑1 (MQ=255)                                       ccTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCa                                                                                         >  1:1700313/1‑150 (MQ=255)                                       ccTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCa                                                                                         >  1:262720/1‑150 (MQ=255)                                        cTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCaa                                                                                        >  1:1083562/1‑150 (MQ=255)                                        cTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACAATGATCGCCTTGGCGGCATGAGGACGGATCTTCTGCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCaa                                                                                        >  1:1670741/1‑150 (MQ=255)                                         ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAc                                                                                       >  1:1560987/1‑150 (MQ=255)                                         ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAc                                                                                       >  2:321566/1‑150 (MQ=255)                                         ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAc                                                                                       >  2:1488150/1‑150 (MQ=255)                                             ccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTgg                                                                                   >  1:1510939/1‑150 (MQ=255)                                                  cAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAAt                                                                              <  2:1207378/150‑1 (MQ=255)                                                    ggCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGc                                                                            >  1:783594/1‑150 (MQ=255)                                                               gCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGt                                                                 <  1:1220268/150‑1 (MQ=255)                                                                ccTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTg                                                                >  2:546330/1‑150 (MQ=255)                                                                 ctgctgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCAATGTCCATGAGgc                                                                                                                               >  2:767650/1‑86 (MQ=38)                                                                 ctgctgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCAATGTCCATGAGgc                                                                                                                               <  1:767650/86‑1 (MQ=38)                                                                     tgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTgg                                                           >  2:1227742/1‑150 (MQ=255)                                                                       ctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGa                                                         <  1:191927/150‑1 (MQ=255)                                                                       ctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGa                                                         <  1:500178/150‑1 (MQ=255)                                                                        tgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATg                                                                                                                                   >  2:525810/1‑75 (MQ=255)                                                                        tgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATg                                                                                                                                   <  1:525810/75‑1 (MQ=255)                                                                            cATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCCTGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCAt                                                    <  2:1215340/150‑1 (MQ=255)                                                                            cATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCAt                                                    <  1:1190890/150‑1 (MQ=255)                                                                             aTCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATg                                                   >  2:1283947/1‑150 (MQ=255)                                                                              tCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCCCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCCCGTAAGCGGGGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATgg                                                  <  1:1789162/150‑1 (MQ=255)                                                                                 gCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGa                                                                                                                                  >  1:1370234/1‑67 (MQ=255)                                                                                 gCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGa                                                                                                                                  <  2:1370234/67‑1 (MQ=255)                                                                                 gCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGa                                               <  1:1322925/150‑1 (MQ=255)                                                                                 gCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGa                                               <  2:1353308/150‑1 (MQ=255)                                                                                  ccTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGAt                                              <  1:1227742/150‑1 (MQ=255)                                                                                      gtggGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACt                                          <  1:1283947/150‑1 (MQ=255)                                                                                      gtggGCACCATGATCGACTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCTTGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACt                                          <  1:1891564/150‑1 (MQ=255)                                                                                       tggGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTg                                         >  2:1390163/1‑150 (MQ=255)                                                                                            aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGc                                    >  1:1749481/1‑150 (MQ=255)                                                                                            aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGc                                    >  1:563931/1‑150 (MQ=255)                                                                                                 gATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATggg                               <  2:1560987/150‑1 (MQ=255)                                                                                                   tCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCt                             <  1:1199118/150‑1 (MQ=255)                                                                                                     gCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCACTGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTgg                           <  2:740049/150‑1 (MQ=255)                                                                                                       cTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAc                         <  1:216143/150‑1 (MQ=255)                                                                                                             ggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGc                              >  1:1822645/1‑139 (MQ=255)                                                                                                             ggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGc                              <  2:1822645/139‑1 (MQ=255)                                                                                                                 ccGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa               <  2:1488326/150‑1 (MQ=255)                                                                                                                 ccGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa               <  1:1488150/150‑1 (MQ=255)                                                                                                                 ccGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa               <  2:1876445/150‑1 (MQ=255)                                                                                                                 ccGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGAGTCGAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCCTGGGCTGGACCAGCACCGGa               <  1:962116/150‑1 (MQ=255)                                                                                                                  cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGAt              >  2:413077/1‑150 (MQ=255)                                                                                                                      aCAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACAt          <  2:262720/150‑1 (MQ=255)                                                                                                                       cAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATc         <  2:1604514/150‑1 (MQ=255)                                                                                                                            ccGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGgcgc    >  2:1722328/1‑150 (MQ=255)                                                                                                                             cGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgcc   >  2:1066800/1‑150 (MQ=255)                                                                                                                              gCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGt                                                                                                           <  2:672942/45‑1 (MQ=255)                                                                                                                              gCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGt                                                                                                           >  1:672942/1‑45 (MQ=255)                                                                                                                              gCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGc                                                 >  1:1117370/1‑103 (MQ=255)                                                                                                                              gCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGc                                                 <  2:1117370/103‑1 (MQ=255)                                                                                                                              gCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccg  >  2:1223760/1‑150 (MQ=255)                                                                                                                                 tCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGAt              >  1:1887715/1‑135 (MQ=255)                                                                                                                                 tCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGAt              <  2:1887715/135‑1 (MQ=255)                                                                                                                                 tCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                      <  2:1255319/79‑1 (MQ=255)                                                                                                                                 tCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                      >  1:1255319/1‑79 (MQ=255)                                                                                                                                 tCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGa                          >  2:760370/1‑123 (MQ=255)                                                                                                                                 tCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGa                          <  1:760370/123‑1 (MQ=255)                                                                                                                                   aTGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCAt                                  <  2:443835/113‑1 (MQ=255)                                                                                                                                   aTGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCAt                                  >  1:443835/1‑113 (MQ=255)                                                                                                                                       |                                                                                                                                             GGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCG  >  NC_002947/3812568‑3812842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑