Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,360,310	0	A	T	11.4%	97.5 / 4.5	35	intergenic (‑279/+86)	cadB/cadC	putative lysine/cadaverine transporter/cadBA operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base T (1/3);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01
Rejected: Not supported by required number of reads on both strands.

TTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATC  >  NC_000913/4360018‑4360602                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                     ttGGCAGAACTCATGCTCTTTTCCTAATTTCATTGTTGGATTTTGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGc                                                                                                                                                                                                                                                                                              <  1:4037/301‑1 (MQ=255)   ggCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCg                                                                                                                                                                                                                                                                                            <  1:293874/301‑1 (MQ=255)    gCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATACGGCAAATAAATTGCCGc                                                                                                                                                                                                                                                                                           >  1:422890/1‑301 (MQ=255)         aCTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAATTTATTTGCCGCataat                                                                                                                                                                                                                                                                                      <  2:305551/301‑5 (MQ=255)            aaTGCTCTTCTCCTAATGTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTACTGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAATTTATTTGCCGCataatttt                                                                                                                                                                                                                                                                                   <  1:600094/300‑8 (MQ=255)             aTGCTCTTCTCCTAATTTCATTTGTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGAACAAACGTGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCGGTAAGTGAGAACGTGAGGTTTTTTATTAACACATCAGGAGCGCAAGTTGATATCATGAAAAGATAAACATTGAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAATTTATTTGCCGCataattgtt                                                                                                                                                                                                                                                                                  <  1:492183/301‑9 (MQ=255)                 tcttctCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATAcc                                                                                                                                                                                                                                                                              <  2:168276/301‑1 (MQ=255)                    tctcCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAAc                                                                                                                                                                                                                                                                           <  2:185217/301‑1 (MQ=255)                      tcCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAg                                                                                                                                                                                                                                                                         <  2:453296/301‑1 (MQ=255)                         tAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGaa                                                                                                                                                                                                                                                                      >  1:579876/1‑301 (MQ=255)                           aTTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGAGGCAAATAAATTGCCGCAACATATTATACCAACAGAAACa                                                                                                                                                                                                                                                                    >  1:415984/1‑301 (MQ=255)                              tCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATAc                                                                                                                                                                                                                                                                 <  2:423453/301‑1 (MQ=255)                                     tGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACt                                                                                                                                                                                                                                                          <  2:464643/301‑1 (MQ=255)                                       aaTTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTca                                                                                                                                                                                                                                                        <  2:75068/301‑1 (MQ=255)                                                     gTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATGTCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTc                                                                                                                                                                                                                                          <  2:605849/301‑1 (MQ=255)                                                       tatgatgTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCg                                                                                                                                                                                                                                        <  2:416547/300‑1 (MQ=255)                                                                taACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAAt                                                                                                                                                                                                                               <  2:580600/301‑1 (MQ=255)                                                                 aaCTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGGTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATc                                                                                                                                                                                                                              <  1:281093/301‑1 (MQ=255)                                                                  aCTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACACAACAACAAATATTTCCAACCATacaaca                                                                                                                                                                                                                             >  1:204189/1‑296 (MQ=255)                                                                          cTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAg                                                                                                                                                                                                                     <  2:427330/301‑1 (MQ=255)                                                                                      aaCTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACCTACAAAAACTCAACAACAAATATTTCCAAGCATAAATCAACCCGGAGATACATAATCtt                                                                                                                                                                                                         >  1:46748/1‑300 (MQ=255)                                                                                                                                                cacaAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATAca                                                                                                                                               >  2:529240/1‑300 (MQ=255)                                                                                                                                                 acaAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAAGAACAGAAGTCTGGAATATACCa                                                                                                                                              >  2:461530/1‑301 (MQ=255)                                                                                                                                                            tGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTAAATCCa                                                                                                                                   >  2:597906/1‑301 (MQ=255)                                                                                                                                                                              tttttttATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGcccc                                                                                                                 <  2:204432/300‑1 (MQ=255)                                                                                                                                                                                tttttATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCACCATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAACTCAACCCGGAGTTACTTATTCTGAAGCAAGAACTTTGTCGAGATAAGGTACAACATAACGAACAGACGTCTGGAATATACCATTATCACACCAGAAAAGGAAGTATACCCCTg                                                                                                               >  1:231205/1‑301 (MQ=255)                                                                                                                                                                                            caTCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTaa                                                                                                   >  2:391416/1‑301 (MQ=255)                                                                                                                                                                                                                     aaaaaGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAgctgc                                                                          <  1:528567/300‑1 (MQ=255)                                                                                                                                                                                                                                                 ggATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATaa                                              >  2:558617/1‑301 (MQ=255)                                                                                                                                                                                                                                                       cGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCttt                                        >  2:176174/1‑300 (MQ=255)                                                                                                                                                                                                                                                            taTTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAg                                   <  2:47783/300‑1 (MQ=255)                                                                                                                                                                                                                                                                ctctGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAac                               <  1:597145/301‑1 (MQ=255)                                                                                                                                                                                                                                                                  ctGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGAATCAGCAGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCACGCAacac                             >  1:411052/1‑301 (MQ=255)                                                                                                                                                                                                                                                                          aTTTATGTTATAAAAATTATGCGGCAAATAAATTGCCGCAAAATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTa                     <  1:557903/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                             tGCGGCAATTTATTTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAACCCTTGCCAAGCAACCCATAATTTAGCCAGGACATTCCAAGATc  >  1:317800/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                     TTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTGAAATAACACAAAAGTTTCTGTAAGTGAGAACTTGAGGTTTTTTATTAACACATCAGGATCGCAAGTTGATATCATGAAAAGATAAACATTTAATGTTTACAATGGATTGCGTGACATTCTCTGGTTAAATTTATGTAATAAAAATTATGCGGCAAATAAATTGCCGCAACATATTATACCAACAGGAACATACAAAAACTCAACAACAAATATTTCCGAGCATAAATCAACCCGGAGTTACTTATTCTGAAGCAAGAAATTTGTCGAGATAAGGTACAACATAAGGAACAGAAGTCTGGAATATACCATTTTCAATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATC  >  NC_000913/4360018‑4360602

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