Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,326	0	T	C	36.7%	53.9 / 81.4	79	pseudogene (51/279 nt)	exoD	pseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base C (16/13);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTG  >  NC_000913/566203‑566453                                                                                                                            |                                                                                                                                atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGaaa                                                                                                                  >  1:380008/4‑139 (MQ=255)   tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAacc                                                                                                                >  2:857485/1‑139 (MQ=255)   cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAacc                                                                                                                >  2:599226/2‑139 (MQ=255)   cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAacc                                                                                                                >  2:307657/2‑139 (MQ=255)   cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAacc                                                                                                                >  1:330568/2‑139 (MQ=255)     cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccac                                                                                                              >  1:141135/1‑139 (MQ=255)           aTCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                        >  1:1066605/1‑139 (MQ=255)            tCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGACATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCg                                                                                                       <  2:713231/139‑1 (MQ=255)             cACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                                                      <  2:141135/139‑1 (MQ=255)               cGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaa                                                                                                    <  2:330568/139‑1 (MQ=255)                 aCATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaact                                                                                                  <  2:1066605/139‑3 (MQ=255)                              cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCg                                                                                     <  1:599226/139‑1 (MQ=255)                                 cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa                                                                                  >  2:493258/1‑139 (MQ=255)                                    gCTTGATGCGCAGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGt                                                                               <  2:96741/139‑1 (MQ=255)                                      ttCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaa                                                                             <  2:714634/139‑2 (MQ=255)                                        cATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGa                                                                           <  2:1117152/139‑1 (MQ=255)                                        cATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                                                      <  2:478917/112‑1 (MQ=255)                                        cATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                                                      >  1:478917/1‑112 (MQ=255)                                          tactcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                        <  1:843527/104‑1 (MQ=255)                                          tacgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTCAGCTGACTTTATGGCCTCGAAACCACCGAGcc                                                                                                        >  2:843527/3‑108 (MQ=255)                                              cgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcat                                                                     >  1:1083345/1‑130 (MQ=255)                                              cgGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTg                                                                     >  2:439346/1‑139 (MQ=255)                                                           gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtggc                                                                                                                         <  2:1298124/74‑4 (MQ=255)                                                           gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGc                                                                                                                         >  1:1298124/1‑74 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTa                                                                                                                             <  1:897471/70‑1 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTa                                                                                                                             >  2:897471/1‑70 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGc                                                                                                                         >  2:155610/1‑74 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGc                                                                                                                         <  1:155610/74‑1 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTGGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  2:941630/1‑139 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACg                                                                            >  2:705360/1‑119 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACg                                                                            <  1:705360/119‑1 (MQ=255)                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTA‑‑‑CACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTg                                                     >  1:46561/1‑139 (MQ=255)                                                               aGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTAGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  1:966424/1‑135 (MQ=255)                                                               aGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTAGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        <  2:966424/135‑1 (MQ=255)                                                                 tACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTg                                                     >  1:79742/1‑136 (MQ=255)                                                                 tACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTg                                                     <  2:79742/136‑1 (MQ=255)                                                                         aTCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                          <  1:439346/139‑1 (MQ=255)                                                                          tCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                         <  1:1266571/139‑1 (MQ=255)                                                                           tttttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGaa                                                                                                                   >  1:87629/2‑64 (MQ=255)                                                                           tttttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGaa                                                                                                                   <  2:87629/63‑1 (MQ=255)                                                                           cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                   >  2:166881/1‑96 (MQ=255)                                                                           cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                   <  1:166881/96‑1 (MQ=255)                                                                             tttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCagag                                      >  2:1081656/1‑139 (MQ=255)                                                                                    tCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        <  1:1153197/114‑1 (MQ=255)                                                                                    tCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  2:1153197/1‑114 (MQ=255)                                                                                      aCCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGcc                                                                                                        <  1:1014595/64‑1 (MQ=255)                                                                                      aCCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGcc                                                                                                        >  2:1014595/1‑64 (MQ=255)                                                                                      aCCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        <  1:388188/112‑1 (MQ=255)                                                                                      aCCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  2:388188/1‑112 (MQ=255)                                                                                               cTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaa                                                                             >  1:436709/1‑81 (MQ=255)                                                                                               cTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        <  1:376480/103‑1 (MQ=255)                                                                                               cTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  2:376480/1‑103 (MQ=255)                                                                                               cTGTACTGCACCTGGGCAATGTAAGCTAACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaa                                                                             <  2:436709/82‑2 (MQ=255)                                                                                                  tACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCg                                                                                     >  2:109526/1‑71 (MQ=255)                                                                                                  tACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCg                                                                                     <  1:109526/71‑1 (MQ=255)                                                                                                  tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGaa                                                                                           <  2:673738/65‑1 (MQ=255)                                                                                                  tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGaa                                                                                           >  1:673738/1‑65 (MQ=255)                                                                                                  tACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  2:817120/1‑100 (MQ=255)                                                                                                  tACTGCACCTGGGCCATGTAAGCCGATTTTACTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        <  1:817120/100‑1 (MQ=255)                                                                                                    cTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcgc               <  2:686930/139‑1 (MQ=255)                                                                                                      gCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                >  1:578592/1‑136 (MQ=255)                                                                                                      gCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                <  2:578592/136‑1 (MQ=255)                                                                                                       cACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaact                                                                                                  >  2:264925/1‑51 (MQ=255)                                                                                                       cACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaact                                                                                                  <  1:264925/53‑3 (MQ=255)                                                                                                            gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                             <  1:89026/53‑1 (MQ=255)                                                                                                            gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                             >  2:89026/1‑53 (MQ=255)                                                                                                            gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGaa                                                                                           <  1:8330/55‑1 (MQ=255)                                                                                                            gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGaa                                                                                           >  2:8330/1‑55 (MQ=255)                                                                                                            gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                         >  2:995804/1‑105 (MQ=255)                                                                                                            gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                         <  1:995804/105‑1 (MQ=255)                                                                                                            gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACt       >  2:1294627/1‑139 (MQ=255)                                                                                                                 aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggttg  >  1:942819/1‑139 (MQ=255)                                                                                                                  tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCg                                                                                                       >  1:1091729/1‑37 (MQ=255)                                                                                                                  tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCg                                                                                                       <  2:1091729/37‑1 (MQ=255)                                                                                                                  tGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATcc                                                                                       >  1:562937/1‑53 (MQ=255)                                                                                                                  tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  1:1002506/1‑84 (MQ=255)                                                                                                                  tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        <  2:1002506/84‑1 (MQ=255)                                                                                                                  tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        <  2:1104778/84‑1 (MQ=255)                                                                                                                  tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                        >  1:1104778/1‑84 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCAt                                                                                              >  2:1102984/1‑36 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCAt                                                                                              <  1:1102984/36‑1 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTa                                                                              >  2:1226401/1‑52 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTa                                                                              <  1:1226401/52‑1 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTa                                                   >  2:1311900/1‑79 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTa                                                   <  1:1311900/79‑1 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                <  1:632008/114‑1 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                >  2:632008/1‑114 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACgcgc             >  2:717279/1‑117 (MQ=255)                                                                                                                            ttttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACgcgc             <  1:717279/117‑1 (MQ=255)                                                                                                                            |                                                                                                                                ATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTG  >  NC_000913/566203‑566453

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