Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,376,482	G→C	9.5%	intergenic (+14/‑134)	PP_2975 → / → PP_2976	insertion element ISR1 protein A3/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,376,482	0	G	C	9.5%	226.4 / 12.6	76	intergenic (+14/‑134)	PP_2975/PP_2976	insertion element ISR1 protein A3/transposase
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base C (3/4);  total (36/40)
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.64e-01

GCGGCATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCATATGAGCAA  >  NC_002947/3376338‑3376623                                                                                                                                                   |                                                                                                                                               gcggcATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑ga                                                                                                                                          <  1:968318/150‑1 (MQ=255)     cATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATg                                                                                                                                      >  2:1417467/1‑150 (MQ=255)      aTGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATgg                                                                                                                                     <  1:1342065/150‑1 (MQ=255)      aTGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATgg                                                                                                                                     <  1:281617/150‑1 (MQ=255)       tGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTC‑GCGCCGCTGAGATGgc                                                                                                                                    <  1:382591/150‑1 (MQ=255)           tGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCgcag                                                                                                                                >  2:2038429/1‑150 (MQ=255)               tgatTGATGACCAGTCTGCGTGAACGGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAgg                                                                                                                            >  2:643649/3‑150 (MQ=255)               tcatTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAgg                                                                                                                            >  2:1679153/1‑150 (MQ=255)               tcatTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAgg                                                                                                                            <  1:2038429/150‑1 (MQ=255)                catTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAATCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGc                                                                                                                           >  2:539422/1‑150 (MQ=255)                          cAGGCTGCGTGACCTGGGAGAGGAGAACCGTCGGCTCAAGAAGATGTAAGCC‑‑GCAACGCTTTGGCGACAAGGCGTGAGCGCGGAGGGCGGGGCAAAACAGTGCTACAGCTACCTCAGCGGGGT‑GGGATGGCGCAGCAGGCGGTTCGTTCa                                                                                                                 <  1:1749802/150‑1 (MQ=255)                          cAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCa                                                                                                                 >  2:2377574/1‑150 (MQ=255)                          cAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCa                                                                                                                 <  1:1431879/150‑1 (MQ=255)                          cAGGCTGCGTGAACTGGAAGAGGAGAACCGTCCGCTCAAGAAGATGTACGCCGAGGAGCGCCTAAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCa                                                                                                                 <  1:1079718/150‑1 (MQ=255)                                  gTGAACTTGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGc                                                                                                         <  2:475215/150‑1 (MQ=255)                                  gTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGc                                                                                                         >  2:849327/1‑150 (MQ=255)                                    gAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCa                                                                                                       <  2:1888503/150‑1 (MQ=255)                                     aaCTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAg                                                                                                      <  1:2377574/150‑1 (MQ=255)                                      aatGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGACTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGt                                                                                                     <  2:549859/148‑1 (MQ=255)                                        tGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTAt                                                                                                   >  1:2117700/1‑150 (MQ=255)                                        tGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTC‑GCGCCGCTGAGATgg                                                                                                                                     >  2:2270570/1‑116 (MQ=255)                                        tGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTC‑GCGCCGCTGAGATgg                                                                                                                                     >  2:1036742/1‑116 (MQ=255)                                        tGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTC‑GCGCCGCTGAGATgg                                                                                                                                     <  1:1036742/116‑1 (MQ=255)                                        tGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTC‑GCGCCGCTGAGATgg                                                                                                                                     <  1:2270570/116‑1 (MQ=255)                                         ggAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATc                                                                                                  <  2:2045467/150‑1 (MQ=255)                                          gAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCa                                                                                                 >  1:2242936/1‑150 (MQ=255)                                              aggagAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAAtt                                                                                             >  1:965464/1‑150 (MQ=255)                                                gagaACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAATGCCGAG‑‑ATCAGCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAg                                                                                           >  2:346353/1‑150 (MQ=255)                                                 agaACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAACAAACTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGc                                                                                          <  2:467044/150‑1 (MQ=255)                                                 agaACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGc                                                                                          <  2:2242936/150‑1 (MQ=255)                                                   aaCCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGt                                                                                        <  1:1471977/150‑1 (MQ=255)                                                      cGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCaa                                                                                                >  1:1394465/1‑139 (MQ=255)                                                      cGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCaa                                                                                                <  2:1394465/139‑1 (MQ=255)                                                      cGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGcc                                                                                     <  1:38044/150‑1 (MQ=255)                                                             tCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATgg                                                                                                                                     <  2:447293/95‑1 (MQ=255)                                                             tCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATgg                                                                                                                                     >  1:447293/1‑95 (MQ=255)                                                             tCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGtt                                                                              >  2:2257098/1‑150 (MQ=255)                                                             tCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGtt                                                                              >  1:417527/1‑150 (MQ=255)                                                             tCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGtt                                                                              <  1:1587951/150‑1 (MQ=255)                                                             tCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGtt                                                                              >  1:544267/1‑150 (MQ=255)                                                               aagaagATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTg                                                                            <  1:2257098/150‑1 (MQ=255)                                                                agaagaTGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTgg                                                                           >  1:1656116/1‑150 (MQ=255)                                                                agaagaTGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTgg                                                                           <  1:539422/150‑1 (MQ=255)                                                                 gaagaTGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGt                                                                          >  1:1047735/1‑150 (MQ=255)                                                                   agaTGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTAt                                                                        <  1:2221352/150‑1 (MQ=255)                                                                          cGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGGGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAACGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACg                                                                 >  1:1768700/1‑150 (MQ=255)                                                                           gCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGa                                                                <  2:1416067/150‑1 (MQ=255)                                                                               aggagCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGa                                                            <  2:1656116/150‑1 (MQ=255)                                                                                   gcgcCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAAc                                                          <  2:1847564/148‑1 (MQ=255)                                                                                   gcgcCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAAc                                                          >  1:1847564/1‑148 (MQ=255)                                                                                     gcCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGcc                                                      >  1:1009911/1‑150 (MQ=255)                                                                                       cTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGcccc                                                    >  2:1931269/1‑150 (MQ=255)                                                                                             gCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCa                                                                      <  2:1997020/126‑1 (MQ=255)                                                                                             gCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCa                                                                      >  1:1997020/1‑126 (MQ=255)                                                                                             gCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGttt                                              >  1:440075/1‑150 (MQ=255)                                                                                                     atcatcCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGa                                        <  1:1931269/150‑1 (MQ=255)                                                                                                         tcCAGGAGGCCATGGCAAAAAGGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCt                                    >  1:2173732/1‑150 (MQ=255)                                                                                                          ccAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTc                                   <  1:1417467/150‑1 (MQ=255)                                                                                                             ggaggCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCAc                                >  2:463701/1‑150 (MQ=255)                                                                                                                  ccATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCttt                           <  1:1679153/150‑1 (MQ=255)                                                                                                                  ccATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCttt                           <  2:1009911/150‑1 (MQ=255)                                                                                                                   cATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTg                          >  1:2377259/1‑150 (MQ=255)                                                                                                                    aTGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTga                         <  1:463701/150‑1 (MQ=255)                                                                                                                      ggCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTgaga                       <  1:1294585/150‑1 (MQ=255)                                                                                                                       gCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTgagag                      <  2:544267/150‑1 (MQ=255)                                                                                                                       gCAAAAAAGTGGTAAAGCCATCTCAGCGAGGT‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTgagag                      >  1:442345/1‑150 (MQ=255)                                                                                                                        cAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTgagaga                     >  1:1245699/1‑150 (MQ=255)                                                                                                                         aaaaaaGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCggg                                                 <  2:2382641/121‑1 (MQ=255)                                                                                                                         aaaaaaGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCggg                                                 >  1:2382641/1‑121 (MQ=255)                                                                                                                         aaaaaaGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGaa                    >  1:1206110/1‑150 (MQ=255)                                                                                                                           aaaaGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATg                  <  2:442345/150‑1 (MQ=255)                                                                                                                                     aaGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCatat        >  2:9986/1‑150 (MQ=255)                                                                                                                                       gCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCATATGa      <  2:2117700/150‑1 (MQ=255)                                                                                                                                       gCCATCTC‑GCGCCGCTGAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGaa                    >  2:2110335/1‑136 (MQ=255)                                                                                                                                       gCCATCTC‑GCGCCGCTGAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGaa                    <  1:2110335/136‑1 (MQ=255)                                                                                                                                           tctcAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCATATGagcaa  >  2:1385828/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                               GCGGCATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAG‑‑ATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCATATGAGCAA  >  NC_002947/3376338‑3376623

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