Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,061,913	A→T	5.2%	Q262L (CAG→CTG)	PP_1838 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,061,913	0	A	T	5.2%	285.3 / 2.3	89	Q262L (CAG→CTG)	PP_1838	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/43);  new base T (2/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGCCGATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATG  >  NC_002947/2061769‑2062061                                                                                                                                                 |                                                                                                                                                     gCCCGCCGATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAgg                                                                                                                                                 >  1:605853/1‑150 (MQ=255)   ccgccgATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGAc                                                                                                                                               >  1:515731/1‑150 (MQ=255)    cgccgATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGAc                                                                                                                                               <  2:1513780/149‑1 (MQ=255)    cgccgATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGAc                                                                                                                                               >  1:1513780/1‑149 (MQ=255)       cgATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCTGCAGGActggc                                                                                                                                           <  1:1074706/150‑5 (MQ=255)        gATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGc                                                                                                                                          >  1:395992/1‑150 (MQ=255)         aTGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGcc                                                                                                                                         <  2:1811541/150‑1 (MQ=255)             gcTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGc                                                                                                                                     <  2:1872465/150‑1 (MQ=255)                ggggTGATGCCTTTATCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCaaa                                                                                                                                  <  2:1352682/150‑1 (MQ=255)                   gTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAgtg                                                                                                                               >  2:455453/1‑150 (MQ=255)                      aTGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGc                                                                                                                            >  1:1162726/1‑150 (MQ=255)                      aTGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGc                                                                                                                            <  1:47006/150‑1 (MQ=255)                        gCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGt                                                                                                                          >  2:279476/1‑150 (MQ=255)                          cTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGTCGCCAAGAGCCGCATTGGCGAGGACCGTGCCAGCATCTTGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGa                                                                                                                        <  1:1493258/150‑1 (MQ=255)                              aCCACGGACTCCAACCTCGCCATCCAGCCGCGGGTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACACAAAGTGTGCGTGAGCaa                                                                                                                    <  1:912335/150‑1 (MQ=255)                                 aCGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAActg                                                                                                                 <  1:279476/150‑1 (MQ=255)                                   tcaCTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAActgct                                                                                                               >  1:1341737/3‑150 (MQ=255)                                     aCTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGa                                                                                                             <  2:395992/150‑1 (MQ=255)                                     aCTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGCGACCAActttaca                                                                                                             <  2:1449919/150‑6 (MQ=255)                                         cAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGAc                                                                                                                                               <  2:1609180/112‑1 (MQ=255)                                         cAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGAc                                                                                                                                               >  1:1609180/1‑112 (MQ=255)                                         cAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGt                                                                                                                          <  2:550715/133‑1 (MQ=255)                                         cAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGt                                                                                                                          >  1:550715/1‑133 (MQ=255)                                             ttCGTCAACCACCTGGGCGTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCGGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCACAAGCACGTCGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCTCAc                                                                                                     <  1:452845/150‑1 (MQ=255)                                                 tCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGaca                                                                                                 >  2:1863545/1‑150 (MQ=255)                                                   aaCCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGCCGCCACCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACAcc                                                                                               <  2:375813/150‑1 (MQ=255)                                                   aaCCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACAcc                                                                                               <  1:1623944/150‑1 (MQ=255)                                                   aaCCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACAcc                                                                                               <  1:600656/150‑1 (MQ=255)                                                     ccAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCt                                                                                             >  2:1774179/1‑150 (MQ=255)                                                      cAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCaa                                                                                                                    >  1:483362/1‑126 (MQ=255)                                                      cAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCaa                                                                                                                    <  2:483362/126‑1 (MQ=255)                                                      cAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCACTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTg                                                                                            >  1:605934/1‑150 (MQ=255)                                                        gCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGt                                                                                          >  2:1544680/1‑150 (MQ=255)                                                             gCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGccc                                                                                                       >  1:324299/1‑132 (MQ=255)                                                             gCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGccc                                                                                                       <  2:324299/132‑1 (MQ=255)                                                             gCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACg                                                                                     >  2:1658356/1‑150 (MQ=255)                                                                       cACCCTGACCCTGATCGACGCCGCCAAGAGCCACTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGc                                                                           <  2:605934/150‑1 (MQ=255)                                                                        accctgaccctgATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGcc                                                                          >  2:1285415/1‑150 (MQ=255)                                                                          cctgaccctgaTCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGc                                                                        >  1:1339888/1‑150 (MQ=255)                                                                                cctgaTCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATccg                                                                  >  1:269727/1‑150 (MQ=255)                                                                                 ctgaTCGACGCCGCCAAGAGCCGCTTTGGTGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATccgc                                                                 >  1:1189286/1‑150 (MQ=255)                                                                                  tgaTCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATccgcc                                                                <  2:605853/150‑1 (MQ=255)                                                                                   gaTCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATccgccg                                                               >  1:786917/1‑150 (MQ=255)                                                                                   gaTCGACGCCGCCAAGAGCCGCTTTGGCAAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATccgccg                                                               <  2:65167/150‑1 (MQ=255)                                                                                        aCGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCACACTTTGCGTggc                                                                                                                       <  1:1314019/89‑3 (MQ=255)                                                                                            cgccAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAActgct                                                                                                               <  2:898734/93‑1 (MQ=255)                                                                                            cgccAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAActgct                                                                                                               >  1:898734/1‑93 (MQ=255)                                                                                             gccAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCTGACGAAGc                                                                           <  1:1752569/128‑1 (MQ=255)                                                                                             gccAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCTGACGAAGc                                                                           >  2:1752569/1‑128 (MQ=255)                                                                                                aaGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCaa                                                                                                                    <  2:1361891/84‑1 (MQ=255)                                                                                                aaGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCaa                                                                                                                    >  1:1361891/1‑84 (MQ=255)                                                                                                    gCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCgg                                              <  2:641705/150‑1 (MQ=255)                                                                                                      cGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATcc                                                                   >  2:1306315/1‑127 (MQ=255)                                                                                                      cGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATcc                                                                   <  1:1306315/127‑1 (MQ=255)                                                                                                        cTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGa                                          >  1:1749344/1‑150 (MQ=255)                                                                                                         tttGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTggc                                         <  1:1863545/150‑3 (MQ=255)                                                                                                               gAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGGAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTAcc                                   >  1:1640934/1‑150 (MQ=255)                                                                                                               gAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTAcc                                   <  2:515731/150‑1 (MQ=255)                                                                                                               gAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTAcc                                   >  2:1919827/1‑150 (MQ=255)                                                                                                                aGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTAccc                                  <  2:1178676/150‑1 (MQ=255)                                                                                                                 ggACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCt                                 <  2:1794476/150‑1 (MQ=255)                                                                                                                  gACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTg                                <  2:269727/150‑1 (MQ=255)                                                                                                                  gACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTg                                >  2:1001733/1‑150 (MQ=255)                                                                                                                   aCCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGc                               <  1:455453/150‑1 (MQ=255)                                                                                                                    ccGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGc                                                                                                                      >  1:186753/1‑62 (MQ=255)                                                                                                                    ccGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGc                                                                                                                      <  2:186753/62‑1 (MQ=255)                                                                                                                    ccGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGcc                              <  2:128070/150‑1 (MQ=255)                                                                                                                     cGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCg                             >  1:286308/1‑150 (MQ=255)                                                                                                                       tGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGAt                           <  2:1162726/150‑1 (MQ=255)                                                                                                                        gCCAACATCTGGAAGCCAGTCCTGCAGCAGGAAGTGGCCACGCAAAGTGTGCGTAAGTCCTTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATc                          <  1:1774179/150‑1 (MQ=255)                                                                                                                           aaCATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAg                       <  1:1919827/150‑1 (MQ=255)                                                                                                                              aTCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACAGCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAAc                    >  1:912931/1‑150 (MQ=255)                                                                                                                              aTCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAAc                    <  1:1285415/150‑1 (MQ=255)                                                                                                                                 tGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCATGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACgtg                 >  2:66519/1‑150 (MQ=255)                                                                                                                                  ggAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACgtgg                >  2:439827/1‑150 (MQ=255)                                                                                                                                  ggAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACgtgg                >  2:463335/1‑150 (MQ=255)                                                                                                                                   gAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGc                                             <  1:676045/120‑1 (MQ=255)                                                                                                                                   gAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGc                                             >  2:676045/1‑120 (MQ=255)                                                                                                                                   gAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACgtggt               >  2:1899163/1‑150 (MQ=255)                                                                                                                                     gccacgTCCTGCTGCAGGAC‑TGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGa                            >  2:246828/6‑134 (MQ=255)                                                                                                                                     gccacgTCCTGCTGCAGGAC‑TGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGa                            <  1:246828/129‑1 (MQ=255)                                                                                                                                       cacgTCCTGCTGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACTCCGACGAAGCCGCTa                                                                      >  2:522040/4‑91 (MQ=255)                                                                                                                                       cacgTCCTGCTGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACTCCGACGAAGCCGCTa                                                                      <  1:522040/88‑1 (MQ=255)                                                                                                                                         aGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGAt         >  2:1500871/1‑150 (MQ=255)                                                                                                                                         aGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGAt         >  2:1100034/1‑150 (MQ=255)                                                                                                                                            ccTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCt      <  1:463335/150‑1 (MQ=255)                                                                                                                                             cTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTg     >  2:836853/1‑150 (MQ=255)                                                                                                                                              tgctgcAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGa    >  1:1214800/1‑150 (MQ=255)                                                                                                                                              tGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGa    <  2:1749344/150‑1 (MQ=255)                                                                                                                                              tGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGa    <  2:786917/150‑1 (MQ=255)                                                                                                                                                cagcagGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGAGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATg  >  2:262554/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCCCGCCGATGCGCTGGGGTGATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATG  >  NC_002947/2061769‑2062061

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