Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,953,617	G→A	8.3%	A796V (GCG→GTG)	recB ←	exodeoxyribonuclease V subunit RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,953,617	0	G	A	8.3%	150.7 / 4.2	48	A796V (GCG→GTG)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base A (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TTAACCTGCCAGGGTTGGTTATCACCAGTTTGTGCCGTTTGCCAGGCAATATCATCATCGCATAACGCTTCAATACAGGTGCGAAGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACAATCTGCACCAGATGTTTATCACTTTCGAGACGCATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGC  >  NC_000913/2953362‑2953868                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                             tcccctgaccAGGGTTGGTTATCACCAGTTTGGGCCGTTTGCCCGGCAATCACCTCACCGCATAACGCTTCAAGACCGGTGCGAAGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTc                                                                                                                                                                                                                                              <  2:963322/264‑1 (MQ=255)                     ttCACCAGTTTTTGCCGTTTTCCAGGCAATATCATCATCGCATAACGCTTCAATACAGGTGCGAAGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTg                                                                                                                                                                                                                          <  1:249310/271‑1 (MQ=255)                             tttGTGCCGTTTGCCAGGCAATATCATCATCGCATAACGCTTCAATACAGGGGCGAAGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCgcagca                                                                                                                                                                                                                  <  1:671258/272‑1 (MQ=255)                             tttGTGCCGTTTGCCACGCAATATCATCATCGCATAACGCTTCAATACAGGTGCGTAGCCCTGCCTCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCgcagca                                                                                                                                                                                                                  <  1:816373/272‑1 (MQ=255)                             tgtcgccgttttGCCCGGCAGTATCACCATGCCATACCGCTGCAATACAGGCGCGAAGCCCTGCCGCATCTTGCGGGTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCgcagca                                                                                                                                                                                                                  <  2:261051/263‑1 (MQ=255)                              ttGTGCCTTTTGCCAGGCAATATCCTCATCGCATAACGCTTCAATACAGGTGCGAAGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCgcagca                                                                                                                                                                                                                  <  1:786624/271‑1 (MQ=255)                                                                aaCGCTTCAATACAGGTGCGAAGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGATGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCAttca                                                                                                                                                                                                              >  1:122516/1‑239 (MQ=255)                                                                                   gAAGCCCTGCCGCAGCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCg                                                                                                                                                                                                                                                              <  2:87177/174‑1 (MQ=255)                                                                                     aGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTg                                                                                                                                                          >  1:868640/1‑272 (MQ=255)                                                                                     aGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCACCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTaa                                                                                                                                                                                                              >  1:403434/1‑220 (MQ=255)                                                                                      gCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCACCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAg                                                                                                                                                                                                             <  2:403433/220‑1 (MQ=255)                                                                                          tGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAgg                                                                                                                                                                                                                                       >  2:729809/1‑190 (MQ=255)                                                                                           gCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGt                                                                                                                                                                                                                                      <  1:729810/190‑1 (MQ=255)                                                                                                                         aaaCGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGGGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCGTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCaa                                                                                                                      <  1:780797/272‑1 (MQ=255)                                                                                                                                                gACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCgcagc                                                                                                                                                                                                                   >  2:975966/1‑156 (MQ=255)                                                                                                                                                 aCGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCgcagca                                                                                                                                                                                                                  <  1:975967/156‑1 (MQ=255)                                                                                                                                                      ggTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGGGATAAACGGCCACCAAACAAATGGGAATTCCAAGCCTTTCAGATTtggg                                                                                          >  2:243834/1‑269 (MQ=255)                                                                                                                                                           cACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGc                                                                                                                                                                                                                                >  2:464715/1‑132 (MQ=255)                                                                                                                                                            aCCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCt                                                                                                                                                                                                                               <  1:464716/132‑1 (MQ=255)                                                                                                                                                               ttttttAGCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACa                                                                                <  1:1126585/271‑1 (MQ=255)                                                                                                                                                                                ggCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACAATCTGCACCAGATGTTTa                                                              <  1:243834/272‑1 (MQ=255)                                                                                                                                                                                                    ccGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCATACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGaa                                                                                                                                                                                                     >  1:132925/1‑118 (MQ=255)                                                                                                                                                                                                     cGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCATACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAAc                                                                                                                                                                                                    <  2:132925/118‑1 (MQ=255)                                                                                                                                                                                                          cTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATcc                                                                                                                                                                                                         >  2:7050/1‑108 (MQ=255)                                                                                                                                                                                                          cTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCaa                                                                                                                      >  2:143500/1‑191 (MQ=255)                                                                                                                                                                                                          cTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCa                                                                                                                       >  1:604956/1‑190 (MQ=255)                                                                                                                                                                                                          cTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACAATCTGCACCAGATGTTTATCACTTTCGAGACGCATttgt                                         >  1:270349/1‑267 (MQ=255)                                                                                                                                                                                                           tGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCa                                                                                                                                                                                                        <  1:7050/108‑1 (MQ=255)                                                                                                                                                                                                            gCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACAATCTGCACCAGATGTTTATCACTTTCGAGACGCATTTGTTGGCTgg                                  <  1:539279/272‑1 (MQ=255)                                                                                                                                                                                                              aaTGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCt                                                                                                                                                                                                                           >  1:945038/1‑86 (MQ=255)                                                                                                                                                                                                               aTGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTg                                                                                                                                                                                                                          <  2:945037/86‑1 (MQ=255)                                                                                                                                                                                                                     aaCCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCaa                                                                                                                      >  2:681028/1‑180 (MQ=255)                                                                                                                                                                                                                      aCCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAAt                                                                                                                     <  1:681029/180‑1 (MQ=255)                                                                                                                                                                                                                             gtgtCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCcgcg                                                                                                                                                                                                                                          >  1:385042/1‑56 (MQ=255)                                                                                                                                                                                                                              tgtCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGa                                                                                                                                                                                                                                         <  2:385041/56‑1 (MQ=255)                                                                                                                                                                                                                                            aaaGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCgg                                                                                                                                                   >  1:417076/1‑128 (MQ=255)                                                                                                                                                                                                                                            aaaGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATa                                                                                                                                       >  1:1103030/1‑140 (MQ=255)                                                                                                                                                                                                                                            aaaGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAg                                                                                                                           >  1:835200/1‑152 (MQ=255)                                                                                                                                                                                                                                             aacaaaaCGCGGATGCTCCGGCCAGCCGTTCGGCTTCCGCGGGGTGGCCGCTTTCTGGCGAAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACAATCTGCACCAGATGTTTATCACTTTCGAGACGCATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGc  <  2:941766/268‑1 (MQ=255)                                                                                                                                                                                                                                             aaGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGa                                                                                                                                                  <  2:417075/128‑1 (MQ=255)                                                                                                                                                                                                                                             aaGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATaa                                                                                                                                      <  2:1103029/140‑1 (MQ=255)                                                                                                                                                                                                                                             aaGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGa                                                                                                                          <  2:835199/152‑1 (MQ=255)                                                                                                                                                                                                                                                 aaaCGCAGAT‑CTTCCACCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAg                                                                                                                           >  1:868600/1‑147 (MQ=255)                                                                                                                                                                                                                                                  aaCGCAGAT‑CTTCCACCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGa                                                                                                                          <  2:868599/147‑1 (MQ=255)                                                                                                                                                                                                                                                      cAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCa                                                                                                                                >  1:567689/1‑137 (MQ=255)                                                                                                                                                                                                                                                      cAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCa                                                                                                                                >  1:893428/1‑137 (MQ=255)                                                                                                                                                                                                                                                      cAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACAATCTGCACCAGATGTTTATCACTTTCGAGACGCATTTGTTGGCTGGAGGCATTACTTTCTGGCTCGAGGATATGttt   >  1:941767/1‑260 (MQ=255)                                                                                                                                                                                                                                                       aGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAg                                                                                                                               <  2:567688/137‑1 (MQ=255)                                                                                                                                                                                                                                                       aGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAg                                                                                                                               <  2:893427/137‑1 (MQ=255)                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                             TTAACCTGCCAGGGTTGGTTATCACCAGTTTGTGCCGTTTGCCAGGCAATATCATCATCGCATAACGCTTCAATACAGGTGCGAAGCCCTGCCGCATCTTGCGGTTCCCCTTTTTGCAGCAAACGCCCGAGCGCACTTTGGTGGACGTCGGTGTCACCTTTTTTATCGCCACGACGGCGCACCAGCGGTGCAACGCCGAGACTGCAATGCCAAACCGAACGTGTCAGCGCCACGTAAAGCAAACGCAGAT‑CTTCCGCCAGACGTTCGGCCTCCGCGAGGTCGACGCTTTCTGGCGCAGCATTAAGATCCAGAACTGCCTCAAACGAGTGGCGATCGTGATAAAACGCCTGCTCCTGGACGCGGAAATTGGTGATAAACGGCAGCCAGACCAATGGATATTCCAGCCCTTTC‑GATTTGTGGATCGTGACAATCTGCACCAGATGTTTATCACTTTCGAGACGCATTTGTTGGCTGGAGGCATTACTGTCTGGCTCGAGGATATGTTGC  >  NC_000913/2953362‑2953868

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