Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	626,572	C→A	6.3%	L208L (CTG→CTT)	dedA ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	626,572	0	C	A	6.3%	247.0 / 3.0	79	L208L (CTG→CTT)	dedA	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base A (2/3);  total (38/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.43e-01

CCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAGCGAAC  >  NC_002947/626424‑626720                                                                                                                                                     |                                                                                                                                                     ccGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                                     <  2:1862240/150‑1 (MQ=255)  cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAg                                                                                                                                                    <  1:1678787/150‑1 (MQ=255)  cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGGTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAg                                                                                                                                                    <  1:1883347/150‑1 (MQ=255)   gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                                   <  1:529710/150‑1 (MQ=255)   gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                                   <  2:481886/150‑1 (MQ=255)    tCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTCCCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgc                                                                                                                                                  <  2:859127/150‑1 (MQ=255)      cAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcgc                                                                                                                                                >  1:131261/1‑150 (MQ=255)      cAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcgc                                                                                                                                                >  2:547339/1‑150 (MQ=255)          aTTGGCAGCCGCTAGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGTCGTACGGCCCAGGCGCCCAc                                                                                                                                            <  2:368673/150‑1 (MQ=255)             ggCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCANgca                                                                                                                                         <  2:1192536/150‑1 (MQ=255)             ggCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACgca                                                                                                                                         <  1:952149/150‑1 (MQ=255)                    gCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGcc                                                                                                                                  <  1:925372/150‑1 (MQ=255)                      tCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCga                                                                                                                                <  2:1868597/150‑1 (MQ=255)                         cgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGaa                                                                                                                             <  2:393983/150‑1 (MQ=255)                               gCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                       <  2:1145232/150‑1 (MQ=255)                                  tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCa                                                                                                                    >  1:110678/1‑150 (MQ=255)                                  tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGATTGCGTCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCa                                                                                                                    <  1:408522/150‑1 (MQ=255)                                   gCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGCGCTTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                                   <  1:1482696/150‑1 (MQ=255)                                     cATTGCTGTGGCGGGATCGGGTAACGCGCCAGAGGACGATTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTAGCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGACGCCCACGCAGCAAGCCGAGAATCAACGGCACCa                                                                                                                 >  2:1479518/1‑150 (MQ=255)                                       ttGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAACAAGCTGTGAAGCATCGGCACCAGc                                                                                                               <  2:169098/150‑1 (MQ=255)                                           tgtgGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACa                                                                                                           <  2:773957/150‑1 (MQ=255)                                            gtgGCGTGATCGGGTAACGCGCCAGAGTGCGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAg                                                                                                          >  1:1793540/1‑150 (MQ=255)                                                cGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatg                                                                                                      >  2:1248905/1‑150 (MQ=255)                                                 gTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCCGCACCAGCGACAGgatga                                                                                                     <  1:65442/150‑1 (MQ=255)                                                    aTCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgatga                                                                                                  >  2:559848/1‑150 (MQ=255)                                                     tCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgatgat                                                                                                 <  2:908211/150‑1 (MQ=255)                                                      cGGGTAACGCGGCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATGATCGGCACCAGCGACAGGatgatgatg                                                                                                <  1:627882/150‑1 (MQ=255)                                                       gggTAACGCTCCAGAGTACGCTTGCCCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                               <  1:1833618/150‑1 (MQ=255)                                                       gggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                               >  2:1427355/1‑150 (MQ=255)                                                                gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCTCATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccaccacca                                                                                      >  1:1190792/1‑150 (MQ=255)                                                                gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccaccacca                                                                                      >  1:1442745/1‑150 (MQ=255)                                                                    gagTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAg                                                                                  >  2:1837967/1‑150 (MQ=255)                                                                    gagTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAg                                                                                  <  2:1479851/150‑1 (MQ=255)                                                                    gagTACGCTTGCACCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCCGCCAGCCGAGAATCATCGGCACCAGCTACAGGATGATGATGACCACCACCATCAg                                                                                  <  1:1185384/150‑1 (MQ=255)                                                                     agTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTAACTCTTCAGTGAGCCTTGGAAGGTGTTGGCAGTCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGc                                                                                 >  1:1671454/1‑150 (MQ=255)                                                                                  cgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCCGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGGGAATGATGCCCACCACCATCAGCGACAGGTgctgct                                                                    >  2:1095494/1‑150 (MQ=255)                                                                                     cAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtga                                                                 <  1:1248789/150‑1 (MQ=255)                                                                                        gcgcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAACCGAGAATCATCGGCACCAGCGACAGGAAGATGATGCCCACCACCATCAGCGACAGGTGCTGCAtgatga                                                              >  1:466254/1‑150 (MQ=255)                                                                                         cgcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCTGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCCCCACCAGCGACAGGTGCTGCTTGATGaa                                                             <  2:213777/150‑1 (MQ=255)                                                                                          gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            >  2:31682/1‑150 (MQ=255)                                                                                          gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTGCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            <  1:26523/150‑1 (MQ=255)                                                                                                tGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATTGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACg                                                      >  1:1090245/1‑150 (MQ=255)                                                                                                tGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACg                                                      >  1:1775674/1‑150 (MQ=255)                                                                                                 gACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGt                                                     >  2:413274/1‑150 (MQ=255)                                                                                                  aCCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGtt                                                    >  1:981103/1‑150 (MQ=255)                                                                                                    cacaTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGc                                                  >  2:621369/1‑150 (MQ=255)                                                                                                      caTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCg                                                <  2:981103/150‑1 (MQ=255)                                                                                                       aTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGACCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGAAGAATGGCACGTTGCCga                                               >  2:192851/1‑150 (MQ=255)                                                                                                         gCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaag                                             <  1:559848/150‑1 (MQ=255)                                                                                                           cgttcgttCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaagaa                                           <  1:621369/150‑1 (MQ=255)                                                                                                             ttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                                >  1:1622514/1‑127 (MQ=255)                                                                                                             ttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                                <  2:1622514/127‑1 (MQ=255)                                                                                                                gttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                                      <  1:1427355/150‑1 (MQ=255)                                                                                                                 ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGcc                                     >  1:1747441/1‑150 (MQ=255)                                                                                                                 ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCAAAACCGACAGTATCTGCGTGAGGAATGGCACGTTGCCGAAGAAGTAGcc                                     >  2:1754079/1‑150 (MQ=255)                                                                                                                 ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGAAAGGATGATGATTCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGcc                                     >  1:1771772/1‑150 (MQ=255)                                                                                                                  tcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGccc                                    >  2:1272077/1‑150 (MQ=255)                                                                                                                       aCCTGTTCAGTGAGCCTTGGCCGTCAGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGt                               <  2:331395/150‑1 (MQ=255)                                                                                                                       aCCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCATCGATAGGATGATGACGCCCACCCCCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGt                               >  1:1077217/1‑150 (MQ=255)                                                                                                                         cTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGa                             <  1:1015525/150‑1 (MQ=255)                                                                                                                            ttCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATAAACGATGGCGACGGCGGGATGAAGGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCa                          >  2:989884/1‑150 (MQ=255)                                                                                                                              cAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAgg                        <  1:1837967/150‑1 (MQ=255)                                                                                                                               aGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGGTGTGCGGGATGAGCGACAGTTGCTGCTTGATGGCTGGCCCCACGGTGACCAGGTATTCCAAGGTGACCAGgc                       >  1:35420/1‑150 (MQ=255)                                                                                                                                gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                                      >  2:209047/1‑134 (MQ=255)                                                                                                                                gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                                      <  1:209047/134‑1 (MQ=255)                                                                                                                                 tGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCCGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGGGACCAGgccg                     >  1:1271154/1‑150 (MQ=255)                                                                                                                                    gCCTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgat                                                                                                       <  2:1498589/65‑1 (MQ=38)                                                                                                                                    gCCTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgat                                                                                                       >  1:1498589/1‑65 (MQ=38)                                                                                                                                     ccTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGa                 <  2:1873560/150‑1 (MQ=255)                                                                                                                                      ctgggCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc                >  2:1721219/4‑150 (MQ=255)                                                                                                                                      ctgggCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc                >  1:598253/4‑150 (MQ=255)                                                                                                                                      cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc                >  1:1766145/1‑150 (MQ=255)                                                                                                                                      cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc                <  2:1442745/150‑1 (MQ=255)                                                                                                                                       ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAcc               >  1:478020/1‑150 (MQ=255)                                                                                                                                           ccGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCaaa           >  1:1658460/1‑150 (MQ=255)                                                                                                                                             gTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAagc         >  1:1284101/1‑150 (MQ=255)                                                                                                                                              tACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCCCCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGACCGAGCCCAAGGTGACCAGGCCGCCGACCCAAagca        <  1:1095494/150‑1 (MQ=255)                                                                                                                                                cGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAagcagc      <  1:1248905/150‑1 (MQ=255)                                                                                                                                                 ggCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCTAGTGGACCAGGCCGCCGACCCAAAGCAGCg     <  1:989884/150‑1 (MQ=255)                                                                                                                                                   cccAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAGCGaa   >  1:1266120/1‑150 (MQ=255)                                                                                                                                                    ccAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCCGCGAAc  <  2:1271154/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAGCGAAC  >  NC_002947/626424‑626720

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