Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,825,939 T→A 9.2% intergenic (‑314/+16) spy ← / ← astE periplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,9390TA 9.2% 120.6 / 2.5 43intergenic (‑314/+16)spy/astEperiplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base A (2/2);  total (19/24)
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.27e-01

AGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGG  >  NC_000913/1825645‑1826238
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                           
aGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTAATAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGc                                                                                                                                                                                                                                                                                                       >  1:642899/1‑301 (MQ=255)
 tttatttatTACTGCTTTCTTAAGTAGCGTGCCTTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCg                                                                                                                                                                                                                                                                                                      <  1:319650/300‑1 (MQ=255)
   atttttttACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCggg                                                                                                                                                                                                                                                                                                    <  1:321824/299‑1 (MQ=255)
    atttattACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGt                                                                                                                                                                                                                                                                                                   >  2:634256/1‑301 (MQ=255)
       tattaCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGAAGAAATCATTCTGCgggagaa                                                                                                                                                                                                                                                                                                >  2:322228/1‑297 (MQ=255)
              gCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCAGCAGAATGATTTATGCGGGTAAGTATTAGc                                                                                                                                                                                                                                                                                         >  1:211748/1‑301 (MQ=255)
                             gTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCa                                                                                                                                                                                                                                                                          <  2:35415/301‑1 (MQ=255)
                                              agttcagacgtgtgctcttccgatcttAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAAt                                                                                                                                                                                                                                                         <  2:264726/276‑1 (MQ=255)
                                                aaCTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATcc                                                                                                                                                                                                                                                       >  1:254851/1‑301 (MQ=255)
                                                     tACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACg                                                                                                                                                                                                                                                  <  2:114953/301‑1 (MQ=255)
                                                      aCTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTGTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGc                                                                                                                                                                                                                                                 <  1:565911/301‑1 (MQ=255)
                                                                  gttaaTTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCAAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGAtt                                                                                                                                                                                                                                     <  1:363524/298‑1 (MQ=255)
                                                                      aTTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAggg                                                                                                                                                                                                                                 <  1:127925/301‑1 (MQ=255)
                                                                        tAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAAtagatcggaagagcgtcgtgtagggaa                                                                                                                                                                                                                               >  1:264458/1‑275 (MQ=255)
                                                                             agagCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAaca                                                                                                                                                                                                                          >  2:279641/1‑301 (MQ=255)
                                                                               agCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                                                                                                                                                                        <  2:211991/301‑1 (MQ=255)
                                                                                gCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACAt                                                                                                                                                                                                                       <  2:41994/301‑1 (MQ=255)
                                                                                                                  ctacacgacgctcttccgatctaaaGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGAGTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctc                                                                                                                                                                                     <  1:164433/279‑1 (MQ=255)
                                                                                                                                      aaaaaGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATg                                                                                                                                                                 >  1:496232/1‑301 (MQ=255)
                                                                                                                                        aaaGAAAATTTGTTGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGtt                                                                                                                                                               <  1:279338/301‑1 (MQ=255)
                                                                                                                                        aaaGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctcagaacggaagagcacacgtcta                                                                                                                                                               >  2:164585/1‑279 (MQ=255)
                                                                                                                                                tttGATGGGGAGTGATGAGAGAATATTACAACACGTTGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtctc                                                                                                                                                       <  1:194786/301‑1 (MQ=255)
                                                                                                                                                       gagaGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATTTTCCTTTCTCAAACGGc                                                                                                                                                >  1:332751/1‑301 (MQ=255)
                                                                                                                                                       gagaGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGc                                                                                                                                                >  2:98981/1‑301 (MQ=255)
                                                                                                                                                                           acaacaCGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACt                                                                                                                            >  2:394269/1‑301 (MQ=255)
                                                                                                                                                                                acGATGATTTTGCATAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCa                                                                                                                       <  2:643776/301‑1 (MQ=255)
                                                                                                                                                                                                       aagaaCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTCTTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgag                                                                                                >  1:424882/1‑301 (MQ=255)
                                                                                                                                                                                                              ttaCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcgag                                                                                         <  2:85068/300‑1 (MQ=255)
                                                                                                                                                                                                              ataCCGGATGATTGGTGATAAATAATGCAAATAACCAGGATTAATTTTTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcgag                                                                                         <  1:98857/301‑1 (MQ=255)
                                                                                                                                                                                                               taCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcgaga                                                                                        >  1:84944/1‑300 (MQ=255)
                                                                                                                                                                                                                           ggTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAACCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGACACTTGCCATATACATTTCGAAGAACGGCGAGTGGCGAGTAATATGCGAAac                                                                            >  1:608186/1‑301 (MQ=255)
                                                                                                                                                                                                                                aaaattaaGCAAATACCTAGGATTTTTCTTTTTTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTTCGGGTAAGTATTAGCTTATTTTTTCTATCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCAccc                                                                       <  1:633379/300‑1 (MQ=255)
                                                                                                                                                                                                                                                 cAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGcggcggg                                                      >  2:59329/1‑300 (MQ=255)
                                                                                                                                                                                                                                                 cAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGACCGGCGgt                                                      >  1:180019/1‑301 (MQ=255)
                                                                                                                                                                                                                                                   ggATTATTCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGgtgt                                                    <  1:59261/301‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaTCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCa                                               >  2:563826/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                     ttATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACAct                                  <  2:255119/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gacgtgtgctcttccgatctcccGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACActc                                 <  2:93139/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCAcc                            <  2:425445/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                              aGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGa                         <  2:333156/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      taTACCCGCATATTGATTTTTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAgcagcg                 <  2:496879/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          cccGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACActcagatcggaagagcgtcgtgt             >  1:93015/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                                     aTTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCAGTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCGCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTgg  <  1:393755/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                           
AGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGG  >  NC_000913/1825645‑1826238

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