Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,000,456	A→C	5.8%	K118N (AAA→AAC)	PP_1786 →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,000,456	0	A	C	5.8%	213.5 / 3.9	78	K118N (AAA→AAC)	PP_1786	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base C (2/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GGCCGTTATGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCA  >  NC_002947/2000307‑2000603                                                                                                                                                      |                                                                                                                                                       ggCCGTTATGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCaaa                                                                                                                                                        <  2:668529/150‑1 (MQ=255)    cGTTATGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑cc                                                                                                                                                   <  1:1024129/150‑1 (MQ=255)         tGCCTCGGCCTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑Ccggcgg                                                                                                                                              <  1:763320/150‑1 (MQ=255)           ccTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGgcggcg                                                                                                                                            >  2:1570366/1‑150 (MQ=255)            cTCGGCTTCGCACGCAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACACCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCaaa                                                                                                                                                        <  2:948402/139‑1 (MQ=255)            cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCaaa                                                                                                                                                        >  1:948402/1‑139 (MQ=255)            cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGt                                                                                                                                           <  2:826369/150‑1 (MQ=255)            cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGt                                                                                                                                           <  2:500057/150‑1 (MQ=255)                     gCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGCAAGTTGACGTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCa                                                                                                                                  >  2:1698416/1‑150 (MQ=255)                     gCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCa                                                                                                                                  >  2:389117/1‑150 (MQ=255)                     gCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCGTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCa                                                                                                                                  >  1:1258612/1‑150 (MQ=255)                      cACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTGTATGCGGCCGGCAGCGGCGTTATGGTCGAGCCGGATGTTGCGACAGTTGACCTGGATTTTTGCTTGTACCTTTAGCCTGCCGTCAACG‑‑CTGGCGGCGTTGAAACTCag                                                                                                                                 <  2:202932/150‑2 (MQ=255)                      cACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCaa                                                                                                                                 >  2:250758/1‑150 (MQ=255)                       aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGaggtcag                                                                                                                                  <  2:1082375/148‑7 (MQ=255)                       aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGaggtcag                                                                                                                                  >  1:1082375/1‑142 (MQ=255)                       aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTc                                                                                                                                   >  2:766755/1‑147 (MQ=255)                       aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTc                                                                                                                                   <  1:766755/147‑1 (MQ=255)                       aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAg                                                                                                                                >  1:1543641/1‑150 (MQ=255)                         gAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCt                                                                                                                              >  2:804844/1‑150 (MQ=255)                         gAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTGTACCCTGACCTCAAAT‑‑CCGGCTGCGTTGAAACTCAAGCt                                                                                                                              >  2:1753365/1‑150 (MQ=255)                              tggtgATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGc                                                                                                                         <  1:1301535/150‑1 (MQ=255)                                  gATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                     >  1:199415/1‑150 (MQ=255)                                  gATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                     <  1:1538126/150‑1 (MQ=255)                                   aTGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGTAAg                                                                                                                    <  1:804844/150‑1 (MQ=255)                                   aTGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCCACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAg                                                                                                                    <  1:577765/150‑1 (MQ=255)                                       ggAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGGTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGc‑aca                                                                                                               >  2:896432/1‑150 (MQ=255)                                           cTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTc                                                                                                           >  1:1063748/1‑150 (MQ=255)                                              aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCaa                                                                                                        <  1:948283/150‑1 (MQ=255)                                              aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCaa                                                                                                        <  2:630573/150‑1 (MQ=255)                                                   tGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAACCTTTGGCGGAAGC‑ACATTTCCAACGCtt                                                                                                   <  1:1698416/150‑1 (MQ=255)                                                       gCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGa                                                                                               >  1:303458/1‑150 (MQ=255)                                                        cAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGag                                                                                                                                       >  1:1095362/1‑109 (MQ=255)                                                        cAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGag                                                                                                                                       <  2:1095362/110‑2 (MQ=255)                                                           ttCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACggg                                                                                           <  2:424492/150‑1 (MQ=255)                                                              tcAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCgg                                                                                        >  2:1447537/1‑150 (MQ=255)                                                                       ccGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGAc                                                                               >  1:1053103/1‑150 (MQ=255)                                                                         gACAGCGCCGTTATGGTCAAGCCGGATGTTGCGACAGGTTAGCTGGAGTGTTACCTGACCCTTTACCCTGAGATAAAAGTACCTTTACCGTT‑‑ATCTGAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGAAGGGCGGGGGGGCGACtt                                                                             >  1:896375/1‑150 (MQ=255)                                                                                ccGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTActc                                                                      >  1:383277/1‑150 (MQ=255)                                                                                 cGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTActct                                                                     <  1:1406723/150‑1 (MQ=255)                                                                                 cGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTActct                                                                     <  1:1447537/150‑1 (MQ=255)                                                                                     aTGGTCAAGCCGGATGTTGCTAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTAt                                                                                                                            >  2:1510617/1‑92 (MQ=255)                                                                                     aTGGTCAAGCCGGATGTTGCTAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTAt                                                                                                                            <  1:1510617/92‑1 (MQ=255)                                                                                         tCAAGCCGGATGTTGCGACAGTTGGCCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCACGCTATTGCGGACGC‑ACACTTCCAACGCTTTGGACGGGCGGAGGGGCGAGTGCCTACTCTCGTTGCCt                                                             >  2:1648582/1‑150 (MQ=255)                                                                                              ccGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCt                                                        <  2:896375/150‑1 (MQ=255)                                                                                                  aTGTTGCGAAAGTTGACGTGGATTTTTACCTGACCTTTTACCCTGTCCTCAAAG‑‑CCGGCGGAGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGTTaaga                                                    <  2:759681/150‑1 (MQ=255)                                                                                                   tGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTaagaa                                                   <  1:250758/150‑1 (MQ=255)                                                                                                       gCGAAAGTTGACCTGGATTTTTACCTGACCCTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCAc                                               <  2:1258612/150‑1 (MQ=255)                                                                                                        cGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCg                                                                                                 <  2:39661/99‑1 (MQ=255)                                                                                                        cGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCg                                                                                                 >  1:39661/1‑99 (MQ=255)                                                                                                         gAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTCGTGTAGC‑ACATTGCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGa                                             >  1:1670465/1‑150 (MQ=255)                                                                                                         gAAAGTTAACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGa                                             <  2:300867/150‑1 (MQ=255)                                                                                                              ttGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTc                                        <  2:1063748/150‑1 (MQ=255)                                                                                                                   cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGaaa                                                                                                                                      >  2:1214106/1‑52 (MQ=255)                                                                                                                   cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGaaa                                                                                                                                      <  1:1214106/52‑1 (MQ=255)                                                                                                                    tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGc                                           >  2:344216/1‑141 (MQ=255)                                                                                                                    tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGc                                           <  1:344216/141‑1 (MQ=255)                                                                                                                     ggATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTaga                                 <  1:389117/150‑1 (MQ=255)                                                                                                                     ggATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTaga                                 <  2:1663441/150‑1 (MQ=255)                                                                                                                        tttttCCCTGACCTTCTGCCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAACCTCGT‑CGGCTTCGCCATCTGCAACTCGTCGGACGGTCGGAGCGGCGCCTCCCTACTATCGGTTGATGGGTTAAGAAGCACGAGCTTCCACTAGAGCt                              <  1:1876933/150‑1 (MQ=255)                                                                                                                              ccTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAGACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTAcc                        >  1:1256951/1‑150 (MQ=255)                                                                                                                               cTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCg                       >  1:1913495/1‑150 (MQ=255)                                                                                                                                tGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                      >  1:1559330/1‑150 (MQ=255)                                                                                                                                tGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                      >  1:215827/1‑150 (MQ=255)                                                                                                                                 gACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCAc                                     <  2:1293995/134‑1 (MQ=255)                                                                                                                                 gACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCAc                                     >  1:1293995/1‑134 (MQ=255)                                                                                                                                   ccTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATCTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACg                   >  2:299367/1‑150 (MQ=255)                                                                                                                                    cgtttACCCGGACCTCAAAT‑‑CCGGCGGCGGGGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACgg                  >  2:504209/3‑150 (MQ=255)                                                                                                                                     tgttACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACggg                 >  2:1268431/3‑150 (MQ=255)                                                                                                                                          cccTAACCTACAAG‑‑CCGGGGGCGGTGAAACGCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAg            <  2:1543641/150‑1 (MQ=255)                                                                                                                                           ccTGCCCTCATAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGAGCGTTCCTACTCTCGATGCCTGGGCTCAGACGCACTAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           <  1:1570366/150‑1 (MQ=255)                                                                                                                                           ccTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACTTTTCGAACGCGTCGGACGTTCTCAGGGGCAAAATCCTGCTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACTGGTTCAGc           >  1:1875333/1‑150 (MQ=255)                                                                                                                                           ccTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGtgcggc           >  1:1765051/1‑145 (MQ=255)                                                                                                                                            cTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGcc          >  2:1474199/1‑150 (MQ=255)                                                                                                                                             tGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCt         <  2:1913495/150‑1 (MQ=255)                                                                                                                                                 cTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTa                          >  2:317317/1‑129 (MQ=255)                                                                                                                                                 cTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTa                          <  1:317317/129‑1 (MQ=255)                                                                                                                                                  tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACgg                  <  1:1168505/136‑1 (MQ=255)                                                                                                                                                  tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACgg                  >  2:1168505/1‑136 (MQ=255)                                                                                                                                                  tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           >  1:1009591/1‑143 (MQ=255)                                                                                                                                                  tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           <  2:1009591/143‑1 (MQ=255)                                                                                                                                                   cAAAG‑‑CCGGCGGTGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTTGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTcagc   <  1:827683/150‑1 (MQ=255)                                                                                                                                                    aaaG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTcagca  <  1:1753365/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       GGCCGTTATGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCA  >  NC_002947/2000307‑2000603

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