Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,326	0	T	C	28.9%	115.9 / 75.1	92	pseudogene (51/279 nt)	exoD	pseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base C (14/12);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-02
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAC  >  NC_000913/566192‑566458                                                                                                                                       |                                                                                                                                     aTGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTa                                                                                                                                  <  1:1205003/139‑1 (MQ=255)   gTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCCTCACGCTTCATACGCGGGTCATAGTTCGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtg                                                                                                                                <  2:1085336/139‑2 (MQ=255)        tttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTc                                                                                                                           >  2:480126/1‑139 (MQ=255)         ttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTATTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCg                                                                                                                          >  1:1761099/1‑139 (MQ=255)         ttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCg                                                                                                                          >  1:1284497/1‑139 (MQ=255)            aTTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGaaa                                                                                                                       >  2:1286650/1‑139 (MQ=255)             tcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAac                                                                                                                      >  2:398921/3‑139 (MQ=255)                cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                                                                                   <  1:1286650/139‑1 (MQ=255)                cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccac                                                                                                                   <  1:453578/139‑1 (MQ=255)                cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccac                                                                                                                   <  1:827450/139‑1 (MQ=255)                               tAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCa                                                                                                    >  1:1217486/1‑139 (MQ=255)                                agtGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACTAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAGACCACCGAGCCGGAACTTCAt                                                                                                   <  1:900528/137‑1 (MQ=255)                                   gCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGaa                                                                                                <  2:1435158/139‑1 (MQ=255)                                       gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATcc                                                                                            >  2:764941/1‑139 (MQ=255)                                       gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATcc                                                                                            >  2:1701266/1‑139 (MQ=255)                                        ccTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATccc                                                                                           >  2:484048/1‑139 (MQ=255)                                         cTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCg                                                                                          <  1:764941/139‑1 (MQ=255)                                         cTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCg                                                                                          <  1:484048/139‑1 (MQ=255)                                         cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCg                                                                                          <  1:808401/139‑1 (MQ=255)                                          ttCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                         <  1:1656965/139‑1 (MQ=255)                                            cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACTCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa                                                                                       >  1:34866/1‑139 (MQ=255)                                            cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                                                                       >  2:566973/1‑139 (MQ=255)                                                 ttCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGA‑TTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGt                                                                                    <  2:1153859/136‑1 (MQ=255)                                                 ttCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGA‑TTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGt                                                                                    >  1:1153859/1‑136 (MQ=255)                                                   catacgcgGGTCATAG‑TGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacg                                                                               >  2:1138376/5‑135 (MQ=255)                                                   cATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                           <  2:1095215/128‑1 (MQ=255)                                                   cATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                                                           >  1:1095215/1‑128 (MQ=255)                                                   cATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaac                                                                                <  2:1778605/139‑4 (MQ=255)                                                         cgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCAGAACTTCATGAAATCCCGGGAGGTaaacgggcat                                                                          >  1:85635/1‑130 (MQ=255)                                                               aTAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGaa                                                                    >  2:1296211/1‑139 (MQ=255)                                                                 aGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGc                                                                  <  2:439400/139‑1 (MQ=255)                                                                      gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                                                                                  <  2:1476012/86‑1 (MQ=255)                                                                      gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                                                                                  >  1:1476012/1‑86 (MQ=255)                                                                      gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                             >  2:777962/1‑91 (MQ=255)                                                                      gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                             <  1:777962/91‑1 (MQ=255)                                                                      gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                             <  1:1701266/139‑1 (MQ=255)                                                                      gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                             >  1:1443356/1‑139 (MQ=255)                                                                      gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                             >  1:1564854/1‑139 (MQ=255)                                                                        aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCt                                                           <  2:1564854/139‑1 (MQ=255)                                                                           gTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTa                                                                                                                                  >  2:1567843/1‑65 (MQ=255)                                                                           gTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTa                                                                                                                                  <  1:1567843/65‑1 (MQ=255)                                                                                     tCTTTTTGCGTCACCCACCTGCTGTACTGCACCTGGGCCATGTAAGCAGATTTTATTGCCTCGAAACCACCTAGCCGGAATTTCATGAAATCCCGGGGGGTACGGGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                              <  1:438980/139‑1 (MQ=255)                                                                                     tCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                              <  2:1125450/139‑1 (MQ=255)                                                                                     tCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAAGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                              <  1:1068925/139‑1 (MQ=255)                                                                                      cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAAGCACCGAGCAGGAATTTCATGAAATCACGGGAGGTACGAGTATTGCCGGAAGCGTGGGCTGTATCCGGATGCag                                             >  1:1234441/1‑139 (MQ=255)                                                                                      cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGACCATGTAAGCCGATTTTAt                                                                                                                                 >  2:394815/1‑55 (MQ=255)                                                                                      cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGACCATGTAAGCCGATTTTAt                                                                                                                                 <  1:394815/55‑1 (MQ=255)                                                                                       ttttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                              <  1:1200950/137‑1 (MQ=255)                                                                                       ttttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                              >  2:1200950/1‑137 (MQ=255)                                                                                           cgcgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGc                                                                                                              <  2:297604/69‑1 (MQ=255)                                                                                           cgcgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGc                                                                                                              >  1:297604/1‑69 (MQ=255)                                                                                             cgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  >  2:1786369/1‑63 (MQ=255)                                                                                             cgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  <  1:1786369/63‑1 (MQ=255)                                                                                               tCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATc                                    <  2:1443356/139‑1 (MQ=255)                                                                                                  cccACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  <  2:763694/58‑1 (MQ=255)                                                                                                  cccACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  >  1:763694/1‑58 (MQ=255)                                                                                                   ccACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                                                                                   <  2:875366/56‑1 (MQ=255)                                                                                                   ccACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                                                                                   >  1:875366/1‑56 (MQ=255)                                                                                                   ccACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTg                                >  2:634695/1‑139 (MQ=255)                                                                                                           tGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGa                        <  2:34866/139‑1 (MQ=255)                                                                                                             tACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCt                                                                                                                            <  1:1202636/37‑1 (MQ=255)                                                                                                             tACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCt                                                                                                                            >  2:1202636/1‑37 (MQ=255)                                                                                                             tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                               >  2:741015/1‑114 (MQ=255)                                                                                                             tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                               <  1:741015/114‑1 (MQ=255)                                                                                                             tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATG‑AATCCCGGGAGGTACGAGTATTGCCg                                                                       >  1:337439/1‑89 (MQ=255)                                                                                                             tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATG‑AATCCCGGGAGGTACGAGTATTGCCg                                                                       <  2:337439/89‑1 (MQ=255)                                                                                                                tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTa                                                        >  2:1170593/1‑102 (MQ=255)                                                                                                                tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTa                                                        <  1:1170593/102‑1 (MQ=255)                                                                                                                tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                   <  1:634695/139‑1 (MQ=255)                                                                                                                 gCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                         >  1:1555608/1‑68 (MQ=255)                                                                                                                 gCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                         <  2:1555608/68‑1 (MQ=255)                                                                                                                    ccTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                     >  1:1168973/1‑133 (MQ=255)                                                                                                                    ccTGGCCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                     <  2:1168973/133‑1 (MQ=255)                                                                                                                     cTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTa              >  1:1327548/1‑139 (MQ=255)                                                                                                                     cTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAA‑TTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTAc             >  2:419823/1‑139 (MQ=255)                                                                                                                       gggCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaac                                                                                                        >  2:1457269/1‑46 (MQ=255)                                                                                                                       gggCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaac                                                                                                        <  1:1457269/47‑2 (MQ=255)                                                                                                                       gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                        >  2:585247/1‑63 (MQ=255)                                                                                                                       gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                        <  1:585247/63‑1 (MQ=255)                                                                                                                       gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACt            >  2:1185529/1‑139 (MQ=255)                                                                                                                         gCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtgg          <  2:783555/139‑1 (MQ=255)                                                                                                                            aTGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCggg                                                                                        <  1:715189/58‑1 (MQ=255)                                                                                                                            aTGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCggg                                                                                        >  2:715189/1‑58 (MQ=255)                                                                                                                            aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAAt                                                                                                        <  1:579494/42‑1 (MQ=255)                                                                                                                            aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAAt                                                                                                        >  2:579494/1‑42 (MQ=255)                                                                                                                             tGTAAGCCGATTTTATTGCCTCGAAACCGCCGAGCCGGAATTTCATGAAATCCCggg                                                                                        <  1:1675867/57‑1 (MQ=255)                                                                                                                             tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                        >  2:1675867/1‑57 (MQ=255)                                                                                                                             tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                         >  2:400005/1‑56 (MQ=255)                                                                                                                             tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                         <  1:400005/56‑1 (MQ=255)                                                                                                                             tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggt         <  2:198467/136‑1 (MQ=255)                                                                                                                             tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggt         >  1:198467/1‑136 (MQ=255)                                                                                                                                aaGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAA‑TTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAc  <  1:419823/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     ATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAC  >  NC_000913/566192‑566458

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