Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,173	0	C	G	46.8%	14.1 / 141.2	79	pseudogene (204/279 nt)	exoD	pseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base G (20/17);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACC  >  NC_000913/566044‑566308                                                                                                                                  |                                                                                                                                        cgcAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCa                                                                                                                                >  2:114891/1‑139 (MQ=255)    aaTGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCAGCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCaaaa                                                                                                                             <  1:587376/139‑1 (MQ=255)          ggCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                       <  2:710603/139‑1 (MQ=255)              tGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGt                                                                                                                   <  1:114891/139‑1 (MQ=255)                ccGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                       <  2:1196070/133‑1 (MQ=255)                ccGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                       >  1:1196070/1‑133 (MQ=255)                         cgcctACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCAtcc                                                                                                        <  1:137786/135‑3 (MQ=255)                             tACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCt                                                                                                    <  2:979837/139‑1 (MQ=255)                             tACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCt                                                                                                    >  2:663073/1‑139 (MQ=255)                               cccGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                       <  2:166016/118‑1 (MQ=255)                               cccGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGcc                                                                                                                       >  1:166016/1‑118 (MQ=255)                                  ggATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATC‑CGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                              <  1:955503/139‑1 (MQ=255)                                   gATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                              <  1:332788/139‑1 (MQ=255)                                        tatCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                         >  1:21547/1‑139 (MQ=255)                                        tatCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                         <  1:408604/139‑1 (MQ=255)                                         atCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGcc                                                                                                                       >  1:1009312/1‑108 (MQ=255)                                         atCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGcc                                                                                                                       <  2:1009312/108‑1 (MQ=255)                                          tCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAAc                                                                                                                            <  2:668119/102‑1 (MQ=255)                                          tCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAAc                                                                                                                            >  1:668119/1‑102 (MQ=255)                                          tCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAAc                                                                                                                            <  2:227517/102‑1 (MQ=255)                                          tCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAAc                                                                                                                            >  1:227517/1‑102 (MQ=255)                                             tGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATg                                                                                    >  1:442201/1‑139 (MQ=255)                                             tGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATg                                                                                    >  1:760679/1‑139 (MQ=255)                                              gAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGc                                                                                   >  2:585341/1‑139 (MQ=255)                                               aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTc                                                                                                   >  1:946126/1‑122 (MQ=255)                                               aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTc                                                                                                   <  2:946126/122‑1 (MQ=255)                                               aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCa                                                                                  >  2:768458/1‑139 (MQ=255)                                               aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGc                                                                                                     >  2:1163962/1‑120 (MQ=255)                                               aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGc                                                                                                     <  1:1163962/120‑1 (MQ=255)                                                 gATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAgg                                                                                >  1:924006/1‑139 (MQ=255)                                                 gATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAgg                                                                                >  2:1334331/1‑139 (MQ=255)                                                   tGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                         >  2:770546/1‑128 (MQ=255)                                                   tGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                         <  1:770546/128‑1 (MQ=255)                                                    gCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCt                                                                             >  2:918114/1‑139 (MQ=255)                                                      gTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGGTGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGAAGGCCTTc                                                                           >  2:812485/1‑139 (MQ=255)                                                      gTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTc                                                                           >  2:1137215/1‑139 (MQ=255)                                                           cGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCt                                                                      >  2:302396/1‑139 (MQ=255)                                                            gCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGt                                                                                                                                  <  2:448224/78‑1 (MQ=255)                                                            gCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGt                                                                                                                                  >  1:448224/1‑78 (MQ=255)                                                                 tGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgc                                                                >  1:322762/1‑139 (MQ=255)                                                                  gCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                               <  2:760679/139‑1 (MQ=255)                                                                       ccAAATACAAAACCAATTTAAGCCCGTGCCCCGCCCATTTTTTCGCTGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGt                                                          >  1:336448/1‑139 (MQ=255)                                                                        cAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGAt                                                                                                                                       <  2:302451/61‑1 (MQ=255)                                                                        cAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGAt                                                                                                                                       >  1:302451/1‑61 (MQ=255)                                                                                aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTgg                                                 >  2:63786/1‑139 (MQ=255)                                                                                aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACttt                                                                                                              <  2:450307/78‑1 (MQ=255)                                                                                aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACttt                                                                                                              >  1:450307/1‑78 (MQ=255)                                                                                aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTgg                                                 >  2:923002/1‑139 (MQ=255)                                                                                aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTgg                                                 >  1:52273/1‑139 (MQ=255)                                                                                  aCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCa                                               <  1:1059703/139‑1 (MQ=255)                                                                                      aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                              >  1:1055194/1‑88 (MQ=255)                                                                                      aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                              <  2:1055194/88‑1 (MQ=255)                                                                                       tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCt                                                                                                    >  2:964853/1‑81 (MQ=255)                                                                                       tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCt                                                                                                    <  1:964853/81‑1 (MQ=255)                                                                                       tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTa                                          >  2:428705/1‑139 (MQ=255)                                                                                            gCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                              <  2:1322154/82‑1 (MQ=255)                                                                                            gCCAGTCCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                              >  1:1322154/1‑82 (MQ=255)                                                                                             ccAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCAtc                                                                                                         <  1:333094/70‑2 (MQ=255)                                                                                             ccAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCAtc                                                                                                         >  2:333094/1‑69 (MQ=255)                                                                                              cAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                         <  2:1187352/85‑1 (MQ=255)                                                                                              cAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                         >  1:1187352/1‑85 (MQ=255)                                                                                                 tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCt                                >  2:572078/1‑139 (MQ=255)                                                                                                 tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCt                                >  1:1251574/1‑139 (MQ=255)                                                                                                 tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCt                                >  1:571079/1‑139 (MQ=255)                                                                                                         ccATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCAt                                                                                                                     <  2:452754/46‑1 (MQ=255)                                                                                                         ccATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCAt                                                                                                                     >  1:452754/1‑46 (MQ=255)                                                                                                           aTTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCa                      <  1:1137215/139‑1 (MQ=255)                                                                                                             tttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCAtacg                                                               <  1:1206827/96‑4 (MQ=255)                                                                                                             tttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                               >  2:1206827/1‑96 (MQ=255)                                                                                                              ttttCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                   <  2:322762/139‑1 (MQ=255)                                                                                                                   gATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATg              <  2:336448/139‑1 (MQ=255)                                                                                                                    aTGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGc             <  1:663073/139‑1 (MQ=255)                                                                                                                     tGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCt            <  1:614298/139‑1 (MQ=255)                                                                                                                       aaCTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCggg                                                           <  1:503499/90‑1 (MQ=255)                                                                                                                       aaCTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCggg                                                           >  2:503499/1‑90 (MQ=255)                                                                                                                         cTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc        <  1:441724/139‑1 (MQ=255)                                                                                                                         cTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc        <  1:1319049/139‑1 (MQ=255)                                                                                                                              gCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAc   <  1:923002/139‑1 (MQ=255)                                                                                                                               cACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc  <  2:1251574/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        CGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACC  >  NC_000913/566044‑566308

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: ‑