Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,332	0	T	G	23.5%	181.3 / 73.7	85	pseudogene (45/279 nt)	exoD	pseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base G (12/8);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGT  >  NC_000913/566199‑566466                                                                                                                                      |                                                                                                                                       tttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTc                                                                                                                                   >  1:502898/1‑139 (MQ=255)  ttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCg                                                                                                                                  >  2:431954/1‑139 (MQ=255)       tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAacc                                                                                                                             >  2:62447/1‑139 (MQ=255)          gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                          >  2:815143/1‑139 (MQ=255)          gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTCGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGTCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                          <  1:674363/139‑1 (MQ=255)          gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                                                                                          >  1:461226/1‑139 (MQ=255)          gCTCAATCACGAAATAATGCAGGCCTTCACGCTTCATACGCGCGACATAGATGGCAAAGTACCCGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGCCTTTATGGCCTCGAAcccacc                                                                                                                          >  1:179433/1‑139 (MQ=255)               aTCTCGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                                     <  1:404651/139‑1 (MQ=255)                    gACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAAt                                                                                                                <  1:431954/139‑1 (MQ=255)                      cATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAAttt                                                                                                              <  1:815143/139‑1 (MQ=255)                        tAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCa                                                                                                            <  2:461226/139‑1 (MQ=255)                         aaTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCAt                                                                                                           <  2:502898/139‑1 (MQ=255)                          aTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                                                                                          <  1:14100/139‑1 (MQ=255)                           tGCAGTCCTTCACGCTTCATGCGCGGGTCATAGTTGGCATAGTACCAGGCCTCTTTTCGTGTCACGCACATGATGTACTGCACCTGGGCCATGTAAGCCGATTTCATTGCCTCGAAACCACCGAGCCGGAATTTCAGGa                                                                                                         >  1:502097/1‑139 (MQ=255)                           tGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGa                                                                                                         >  2:768527/1‑139 (MQ=255)                           tGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGa                                                                                                         >  1:65243/1‑139 (MQ=255)                           tGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGa                                                                                                         >  1:17975/1‑139 (MQ=255)                                   ttCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                                 <  2:606535/139‑1 (MQ=255)                                    tCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                                <  1:367178/139‑1 (MQ=255)                                     cACGCTTCATGCGCGGGTCGTAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa                                                                                               >  2:738785/1‑139 (MQ=255)                                     cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa                                                                                               >  2:915354/1‑139 (MQ=255)                                     cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                                                                               >  2:1005066/1‑139 (MQ=255)                                     cACGCTTCATACGCGGGGCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                                                                               >  2:694402/1‑139 (MQ=255)                                               gcgcGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTa                                                                                     <  1:408813/139‑1 (MQ=255)                                               acgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                                                                     <  2:729241/138‑7 (MQ=255)                                                  cgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcat                                                                                  >  1:235044/1‑130 (MQ=255)                                                   gggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGc                                                                                 <  1:993722/139‑1 (MQ=255)                                                    ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                                     <  1:960813/102‑1 (MQ=255)                                                    ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                                     >  2:960813/1‑102 (MQ=255)                                                    ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                                   <  2:1063889/120‑1 (MQ=255)                                                    ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                                   >  1:1063889/1‑120 (MQ=255)                                                    ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                                                                                >  2:332952/1‑139 (MQ=255)                                                      tCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcatttca                                                                              >  2:547038/1‑126 (MQ=255)                                                           gTTGGCAAAGTACCAGACATCTTTTCGCGCCACCCACATGCTATACTTCACCTGGTCCATGTAATCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCg                                                                         <  1:1069117/139‑1 (MQ=255)                                                            ttGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGt                                                                        <  1:62447/139‑1 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                                <  2:521481/112‑1 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                                >  1:521481/1‑112 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                     <  1:659417/139‑1 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                     <  1:915354/139‑1 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                     <  2:496538/139‑1 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                     <  2:1082014/139‑1 (MQ=255)                                                               gCAAAGGACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                     <  2:930882/139‑1 (MQ=255)                                                                    gTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGAACGAGTATTGCCGGAAGCGTGGCCTGTa                                                                >  2:665494/1‑139 (MQ=255)                                                                        cAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCg                                                            >  2:531111/1‑139 (MQ=255)                                                                         aGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCgg                                                           <  2:426485/139‑1 (MQ=255)                                                                         aGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCgg                                                           >  2:363495/1‑139 (MQ=255)                                                                          ggCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGa                                                          >  2:187009/1‑139 (MQ=255)                                                                              tCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                                      >  2:221171/1‑139 (MQ=255)                                                                               cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                                      >  2:865165/1‑138 (MQ=255)                                                                               cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                                      <  1:865165/138‑1 (MQ=255)                                                                                                     gTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGaa                               <  2:919932/139‑1 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt                                                                                                                                    <  2:16763/37‑1 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt                                                                                                                                    >  1:16763/1‑37 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                                      <  1:492370/115‑1 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                                      >  2:492370/1‑115 (MQ=255)                                                                                                       aCTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAATCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                             <  1:665494/139‑1 (MQ=255)                                                                                                       aCTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                             <  2:992139/139‑1 (MQ=255)                                                                                                         tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                           <  1:221171/139‑1 (MQ=255)                                                                                                         tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                           <  1:738785/139‑1 (MQ=255)                                                                                                           cACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCa                         >  2:1075158/1‑139 (MQ=255)                                                                                                             ccTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                     >  1:88515/1‑93 (MQ=255)                                                                                                             ccTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                     <  2:88515/93‑1 (MQ=255)                                                                                                                gggCCATGTAAGCCGATTTTATTGCCTCGAAACCATCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGa                                                                             <  2:935869/82‑1 (MQ=255)                                                                                                                gggCCATGTAAGCCGATTTTATTGCCTCGAAACCATCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGa                                                                             >  1:935869/1‑82 (MQ=255)                                                                                                                gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAATTCGAACGCGCATTACt                    >  2:752743/1‑139 (MQ=255)                                                                                                                gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACt                    >  2:434347/1‑139 (MQ=255)                                                                                                                gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACt                    >  1:244711/1‑139 (MQ=255)                                                                                                                gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACt                    >  2:249311/1‑139 (MQ=255)                                                                                                                    aaTGTAAGCCTATTTTATTGCCTCGAAACAACCGAGCCGGAATTTCATGCAATCCCAGGAGGTACGGGTATTGCTGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCAAACGCGCATTACtggtt                <  1:249311/138‑1 (MQ=255)                                                                                                                     aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                          >  1:389149/1‑48 (MQ=255)                                                                                                                     aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                          <  2:389149/48‑1 (MQ=255)                                                                                                                     aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggttg               >  2:195910/1‑139 (MQ=255)                                                                                                                      tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCggg                                                                                                <  1:901056/57‑1 (MQ=255)                                                                                                                      tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCggg                                                                                                >  2:901056/1‑57 (MQ=255)                                                                                                                      tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                                                                     <  1:759378/68‑7 (MQ=255)                                                                                                                      tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                                                                     >  2:759378/1‑62 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGACACCACCAAGCCGGAATTTCATGAAATCCCGGGAGGTCCGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGGCTTATCCGTGGAAATAGAACGCGCATTACtggtcgg              >  1:971023/1‑139 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGCGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggttgg              >  1:764868/1‑139 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCg                                                                         <  2:950897/80‑1 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCg                                                                         >  1:950897/1‑80 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTg                                                                  <  2:171811/87‑1 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTg                                                                  >  1:171811/1‑87 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggtt                >  1:538807/1‑137 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggtt                <  2:538807/137‑1 (MQ=255)                                                                                                                            ccGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGCGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCa        >  1:555332/1‑139 (MQ=255)                                                                                                                             cGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTTTATCCGGGTGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACaaa       >  2:39792/1‑139 (MQ=255)                                                                                                                                  ttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGt  <  1:531111/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       TTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGT  >  NC_000913/566199‑566466

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