Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,173	0	C	G	57.0%	51.5 / 144.0	86	pseudogene (204/279 nt)	exoD	pseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (18/19);  new base G (23/26);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTA  >  NC_000913/566053‑566301                                                                                                                          |                                                                                                                                 ggCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                <  2:229871/139‑1 (MQ=255)       ccGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                <  1:1060894/133‑1 (MQ=255)       ccGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                >  2:1060894/1‑133 (MQ=255)       ccGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTAc                                                                                                          >  1:342225/1‑139 (MQ=255)                  ggTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTcc                                                                                                <  1:633115/139‑1 (MQ=255)                  cctACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCg                                                                                               <  2:332324/137‑1 (MQ=255)                      cccGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCAtt                                                                                                  <  1:235887/132‑1 (MQ=255)                      cccGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCAtt                                                                                                  >  2:235887/1‑132 (MQ=255)                      cccGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCa                                                                                           >  2:928238/1‑139 (MQ=255)                             attatCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGa                                                                                    >  2:365485/1‑139 (MQ=255)                              ttatCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                       <  2:956137/135‑1 (MQ=255)                              ttatCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                       >  1:956137/1‑135 (MQ=255)                               tatCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACa                                                                                  <  2:548259/139‑1 (MQ=255)                                 tCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCg                                                                                               >  2:769827/1‑124 (MQ=255)                                 tCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCg                                                                                               <  1:769827/124‑1 (MQ=255)                                 tCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAAt                                                                                         <  2:527696/130‑1 (MQ=255)                                 tCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAAt                                                                                         >  1:527696/1‑130 (MQ=255)                                 tCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATa                                                                                >  1:72085/1‑139 (MQ=255)                                    tGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATg                                                                             >  2:75669/1‑139 (MQ=255)                                    tGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATg                                                                             >  1:706151/1‑139 (MQ=255)                                    tGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGc                                                                                              <  2:627391/122‑1 (MQ=255)                                    tGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGc                                                                                              >  1:627391/1‑122 (MQ=255)                                     gAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCAtc                                                                                                  >  1:202240/1‑116 (MQ=255)                                     gAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCAtc                                                                                                  <  2:202240/117‑2 (MQ=255)                                     gAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGc                                                                            >  2:976162/1‑139 (MQ=255)                                     gAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACAACATAATGc                                                                            >  1:440397/1‑139 (MQ=255)                                      aGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGAc                                                                                   <  1:25245/131‑1 (MQ=255)                                      aGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGAc                                                                                   >  2:25245/1‑131 (MQ=255)                                      aGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCa                                                                           >  2:331387/1‑139 (MQ=255)                                          tGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGcc                                                                       >  1:239482/1‑139 (MQ=255)                                           gCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGt                                                                                                            <  2:762719/101‑1 (MQ=255)                                           gCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACa                                                                                  >  2:945516/1‑127 (MQ=255)                                           gCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACa                                                                                  <  1:945516/127‑1 (MQ=255)                                           gCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCt                                                                      >  2:628138/1‑139 (MQ=255)                                           gCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGAGGAACTCCGGCACGATCTCGTCAAAACTCGCCATGt                                                                                                            >  1:762719/1‑101 (MQ=255)                                              tCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCa                                                                   <  2:814205/139‑1 (MQ=255)                                              tCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCa                                                                   <  1:365485/139‑1 (MQ=255)                                                  cGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATg                                                                                                             <  2:170195/93‑1 (MQ=255)                                                  cGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATg                                                                                                             >  1:170195/1‑93 (MQ=255)                                                        tGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATg                                                                                                             >  2:201262/1‑87 (MQ=255)                                                        tGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATg                                                                                                             <  1:201262/87‑1 (MQ=255)                                                        tGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                  <  1:628315/114‑1 (MQ=255)                                                        tGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                  >  2:628315/1‑114 (MQ=255)                                                        tGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgc                                                         >  2:932558/1‑139 (MQ=255)                                                         gCT‑CCCCAAATACAAAACCACTTCCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAATCTCGCCATGTACTTTTCATTCCGCGCAATCACGACAGAATGCAGGCCTTCACGCTTCATgcg                                                        <  1:1014887/139‑1 (MQ=255)                                                         gCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCAtacg                                                        <  2:570318/139‑4 (MQ=255)                                                         gCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                        <  2:521061/139‑1 (MQ=255)                                                          cT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacgc                                                       <  1:975588/139‑1 (MQ=255)                                                                cAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATg                                                          <  2:984770/131‑1 (MQ=255)                                                                cAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATg                                                          >  1:984770/1‑131 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTgg                                          >  2:799064/1‑139 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTgg                                          >  2:812902/1‑139 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                       >  1:877843/1‑94 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                       <  2:877843/94‑1 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                        >  2:623088/1‑125 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                        <  1:623088/125‑1 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCAt                                                <  2:364088/133‑1 (MQ=255)                                                                        aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCAt                                                >  1:364088/1‑133 (MQ=255)                                                                          aCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCa                                                                                                               <  2:420179/68‑1 (MQ=255)                                                                          aCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCa                                                                                                               >  1:420179/1‑68 (MQ=255)                                                                          aCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCa                                        <  1:976162/139‑1 (MQ=255)                                                                              aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCa                                                                                                    <  2:716354/75‑1 (MQ=255)                                                                              aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCa                                                                                                    >  1:716354/1‑75 (MQ=255)                                                                               tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCa                                                                                                               <  2:1116884/63‑1 (MQ=255)                                                                               tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCa                                                                                                               >  1:1116884/1‑63 (MQ=255)                                                                                         taccTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCt                         >  2:1041016/3‑139 (MQ=255)                                                                                         tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTAGCAAAGTACCAGGCATCt                         >  2:1004796/1‑139 (MQ=255)                                                                                             tCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                       <  2:821118/73‑1 (MQ=255)                                                                                             tCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                       >  1:821118/1‑73 (MQ=255)                                                                                              cGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTcg                    <  1:991447/139‑1 (MQ=255)                                                                                               gTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacgc                                                       >  1:866167/1‑103 (MQ=255)                                                                                               gTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacgc                                                       <  2:866167/103‑1 (MQ=255)                                                                                                tCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGcc                                                                                                                >  1:184623/1‑45 (MQ=255)                                                                                                tCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGcc                                                                                                                <  2:184623/45‑1 (MQ=255)                                                                                                 ccATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                 <  1:785672/139‑1 (MQ=255)                                                                                                 ccATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACAACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                 <  2:440397/139‑1 (MQ=255)                                                                                                    ttttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCAc                                                                  <  2:98272/87‑1 (MQ=255)                                                                                                    ttttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCAc                                                                  >  1:98272/1‑87 (MQ=255)                                                                                                    ttttttCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAc              <  1:1041016/139‑1 (MQ=255)                                                                                                     tttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATg                                                                                                             <  1:45865/43‑1 (MQ=255)                                                                                                     tttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATg                                                                                                             >  2:45865/1‑43 (MQ=255)                                                                                                      ttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc            <  2:488923/139‑1 (MQ=255)                                                                                                        ttCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcac          <  1:928238/139‑1 (MQ=255)                                                                                                         tCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcaca         <  2:551269/139‑1 (MQ=255)                                                                                                           gATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATg       <  1:374404/139‑1 (MQ=255)                                                                                                                aCTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTa  <  2:239482/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                 GGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCT‑CCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTA  >  NC_000913/566053‑566301

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