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

TTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAC  >  NC_000913/566199‑566458                                                                                                                                      |                                                                                                                               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                                                                                                                  >  1:1476012/1‑86 (MQ=255)                                                               gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                                                                                  <  2:1476012/86‑1 (MQ=255)                                                               gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                             >  2:777962/1‑91 (MQ=255)                                                               gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                             <  1:777962/91‑1 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                             >  1:1443356/1‑139 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                             >  1:1564854/1‑139 (MQ=255)                                                               gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                             <  1:1701266/139‑1 (MQ=255)                                                                 aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCt                                                           <  2:1564854/139‑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)                                                                                ttttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                              <  1:1200950/137‑1 (MQ=255)                                                                                ttttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                              >  2:1200950/1‑137 (MQ=255)                                                                                    cgcgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGc                                                                                                              >  1:297604/1‑69 (MQ=255)                                                                                    cgcgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGc                                                                                                              <  2:297604/69‑1 (MQ=255)                                                                                      cgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  <  1:1786369/63‑1 (MQ=255)                                                                                      cgTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  >  2:1786369/1‑63 (MQ=255)                                                                                        tCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATc                                    <  2:1443356/139‑1 (MQ=255)                                                                                           cccACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  >  1:763694/1‑58 (MQ=255)                                                                                           cccACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                                  <  2:763694/58‑1 (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                                                                                                                            >  2:1202636/1‑37 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCt                                                                                                                            <  1:1202636/37‑1 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                               >  2:741015/1‑114 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                               <  1:741015/114‑1 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATG‑AATCCCGGGAGGTACGAGTATTGCCg                                                                       <  2:337439/89‑1 (MQ=255)                                                                                                      tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATG‑AATCCCGGGAGGTACGAGTATTGCCg                                                                       >  1:337439/1‑89 (MQ=255)                                                                                                         tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTa                                                        <  1:1170593/102‑1 (MQ=255)                                                                                                         tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTa                                                        >  2:1170593/1‑102 (MQ=255)                                                                                                         tGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgcg                   <  1:634695/139‑1 (MQ=255)                                                                                                          gCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                         <  2:1555608/68‑1 (MQ=255)                                                                                                          gCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                         >  1:1555608/1‑68 (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         >  1:198467/1‑136 (MQ=255)                                                                                                                      tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggt         <  2:198467/136‑1 (MQ=255)                                                                                                                         aaGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAA‑TTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAc  <  1:419823/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                               TTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAC  >  NC_000913/566199‑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: ‑