Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,323	0	C	T	35.7%	116.3 / 93.3	98	pseudogene (54/279 nt)	exoD	pseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base T (19/16);  total (55/43)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACC  >  NC_000913/566185‑566459                                                                                                                                            |                                                                                                                                         aCTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGcc                                                                                                                                          <  2:25041/139‑1 (MQ=255)         tGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAt                                                                                                                                  <  2:1066239/139‑1 (MQ=255)            aCTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATAAGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAtggc                                                                                                                               <  2:435473/139‑4 (MQ=255)                   atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGaaa                                                                                                                        >  2:950162/4‑139 (MQ=255)                   atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGaaa                                                                                                                        >  2:1108714/4‑139 (MQ=255)                     cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAacc                                                                                                                      >  2:38969/2‑139 (MQ=255)                        gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                                                                                   <  2:19326/139‑1 (MQ=255)                          tCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAt                                                                                                                                  >  1:797016/1‑122 (MQ=255)                                 cGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGaa                                                                                                          >  2:678664/1‑139 (MQ=255)                                   aCATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAAtt                                                                                                        <  2:900455/139‑1 (MQ=255)                                    cATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAActt                                                                                                       >  1:881887/1‑139 (MQ=255)                                              gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                                                                                    >  1:363169/1‑116 (MQ=255)                                              gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                                                                                    <  2:363169/116‑1 (MQ=255)                                              gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATcc                                                                                             >  2:1094974/1‑139 (MQ=255)                                               ccTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATTTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                            <  1:835964/139‑1 (MQ=255)                                               ccTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                            <  1:963424/139‑1 (MQ=255)                                                  tCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                         <  1:945549/139‑1 (MQ=255)                                                   cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGc                                                                                                                               >  2:51055/1‑100 (MQ=255)                                                   cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGc                                                                                                                               <  1:51055/100‑1 (MQ=255)                                                   cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                                                                        >  1:712908/1‑139 (MQ=255)                                                        ttCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaa                                                                                   >  2:329498/1‑138 (MQ=255)                                                             acgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                                                              <  2:1032599/138‑7 (MQ=255)                                                                cgGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTAtt                                                                            >  1:112457/1‑139 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                         >  2:1210932/1‑128 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                         >  2:423621/1‑128 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                         <  1:38969/139‑12 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                         <  2:799207/139‑12 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgg                                                                               <  1:1024223/133‑6 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgg                                                                               >  2:1024223/1‑128 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                          >  2:1205831/1‑122 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                          <  1:1205831/122‑1 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                                                                         >  1:62567/1‑139 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                                                                         >  1:419293/1‑139 (MQ=255)                                                                  ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                                                                         >  2:818210/1‑139 (MQ=255)                                                                    tCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcatttca                                                                       >  1:180466/1‑126 (MQ=255)                                                                             gCAAAGTACCAGGCGTCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTACCTCGAAACCACCGAGCCGGAATTTCGTGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              >  1:751061/1‑139 (MQ=255)                                                                             gCAAAGTACCAGGCATTTTTTCGCGTCTCCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAtg                                                                                                                                 <  2:877931/72‑2 (MQ=255)                                                                             gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTa                                                                                                                                   >  1:921310/1‑70 (MQ=255)                                                                             gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTa                                                                                                                                   <  2:921310/70‑1 (MQ=255)                                                                             gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAtg                                                                                                                                 >  1:877931/1‑71 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                   >  1:881257/1‑102 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                   <  2:881257/102‑1 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                         <  2:850844/112‑1 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                         >  2:106105/1‑112 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                         <  1:106105/112‑1 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                         >  1:850844/1‑112 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              >  1:886009/1‑139 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              >  1:1191139/1‑139 (MQ=255)                                                                             gCAAAGTACCAGGCATCTTTGCGCGTCACCCACATGCTGTACTGCA‑CCTGTGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAAGCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              <  2:608703/139‑1 (MQ=255)                                                                                  gTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                                                                                   <  2:1183496/97‑1 (MQ=255)                                                                                  gTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                                                                                   >  1:1183496/1‑97 (MQ=255)                                                                                             cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCg                                                                        <  2:316408/113‑1 (MQ=255)                                                                                             cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCg                                                                        >  1:316408/1‑113 (MQ=255)                                                                                             cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              >  2:38175/1‑123 (MQ=255)                                                                                             cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              <  1:38175/123‑1 (MQ=255)                                                                                                   gcgTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTa                                                                                                                                   >  1:332917/1‑48 (MQ=255)                                                                                                   gcgTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTa                                                                                                                                   <  2:332917/48‑1 (MQ=255)                                                                                                    cgTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTa                                       <  2:751061/139‑1 (MQ=255)                                                                                                       cACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              <  2:1159032/113‑1 (MQ=255)                                                                                                       cACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              >  1:1159032/1‑113 (MQ=255)                                                                                                                 cTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCg                          >  2:169610/1‑139 (MQ=255)                                                                                                                    tACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGcc                                                                                                                              >  2:371704/1‑36 (MQ=255)                                                                                                                    tACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGcc                                                                                                                              <  1:371704/36‑1 (MQ=255)                                                                                                                    tACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                   >  1:916540/1‑63 (MQ=255)                                                                                                                    tACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                   <  2:916540/63‑1 (MQ=255)                                                                                                                    tACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAc                       >  2:633522/1‑139 (MQ=255)                                                                                                                    tACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAc                       >  2:398075/1‑139 (MQ=255)                                                                                                                     aCTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                      <  2:62567/139‑1 (MQ=255)                                                                                                                       tGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                    <  1:818210/139‑1 (MQ=255)                                                                                                                       tGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                    <  1:633522/139‑1 (MQ=255)                                                                                                                        gCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGaa                                                                                                 <  1:1089392/61‑1 (MQ=255)                                                                                                                        gCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGaa                                                                                                 >  2:1089392/1‑61 (MQ=255)                                                                                                                        gCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGaa                                                                                                 >  1:476915/1‑61 (MQ=255)                                                                                                                        gCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGaa                                                                                                 <  2:476915/61‑1 (MQ=255)                                                                                                                         cA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCgg                                                    >  1:435315/1‑105 (MQ=255)                                                                                                                         cA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCgg                                                    <  2:435315/105‑1 (MQ=255)                                                                                                                         cA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCa                  >  2:1119992/1‑139 (MQ=255)                                                                                                                         cA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCa                  <  2:1191139/139‑1 (MQ=255)                                                                                                                            ccTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGGGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCAtt                >  2:120457/1‑139 (MQ=255)                                                                                                                              tGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTAc              <  2:112457/139‑1 (MQ=255)                                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              <  1:777339/90‑1 (MQ=255)                                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                              >  2:777339/1‑90 (MQ=255)                                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                                <  2:989907/104‑1 (MQ=255)                                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                                >  1:989907/1‑104 (MQ=255)                                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTAt                                      >  1:1161844/1‑114 (MQ=255)                                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTAt                                      <  2:1161844/114‑1 (MQ=255)                                                                                                                                  ccATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggt          >  2:610104/1‑139 (MQ=255)                                                                                                                                  ccATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggt          >  2:966147/1‑139 (MQ=255)                                                                                                                                   cATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggtt         <  1:398075/139‑1 (MQ=255)                                                                                                                                     tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                            >  1:225205/1‑54 (MQ=255)                                                                                                                                     tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                            <  2:225205/54‑1 (MQ=255)                                                                                                                                     tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                                >  2:796334/1‑98 (MQ=255)                                                                                                                                     tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                                <  1:796334/98‑1 (MQ=255)                                                                                                                                      gTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggttggt      >  2:679936/1‑139 (MQ=255)                                                                                                                                       tAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTggttggtt     >  1:879576/1‑139 (MQ=255)                                                                                                                                          gCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                    >  2:741879/1‑121 (MQ=255)                                                                                                                                          gCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                    <  1:741879/121‑1 (MQ=255)                                                                                                                                          gCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAcc  >  2:1196349/1‑139 (MQ=255)                                                                                                                                          gCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAcc  >  2:435741/1‑139 (MQ=255)                                                                                                                                            |                                                                                                                                         ACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACC  >  NC_000913/566185‑566459

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