Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,277,497 G→A 14.4% V55M (GTG→ATG)  yhaV → toxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,277,4970GA14.4% 173.8 / 3.8 73V55M (GTG→ATG) yhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (3/58);  new base A (9/2);  total (12/61)
Fisher's exact test for biased strand distribution p-value = 1.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

AAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATA  >  NC_000913/3277275‑3277716
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
aaagaacttggtgtggggaaTGGACGTAAACTTTGAGTTGGATATTGGCGCCGACTATTCATGGATTTCCCACAGAGGTTAAAGGGTTGGGCGTTATATTCTCATCCCTTGTTGCAGGAAACCAACGACTCTTTATTTGCCAAAGTCGATACATTAAAGGGAAAATATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTg                                                                                                                                                                                                                            <  1:741621/210‑1 (MQ=255)
  agaaaCTTTGCGCTGGCAGGGACGTAAACAGTGATGATGGGATGGGACACAGCGAATAAGGGATTTTCCACAAGGGGTTAATGGTTGGGGGTTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATAATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgag                                                                                                                                                                                                                          <  2:412583/229‑1 (MQ=255)
    aatcTTGTCGCTGGCAAGGATGTCAACATTGATGATGAGGTTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGATATAATCTCATCCTTGTTTACAGGAAACCTACGACGTTTTAGTTGCCGAAGCAGAAACATTAAAGGAAAAAGATCCTGAAAATTAACAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgagga                                                                                                                                                                                                                        <  1:149928/226‑1 (MQ=255)
            cGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTCCGACGATTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCa                                                                                                                                                                                                                <  1:53544/229‑1 (MQ=255)
                      tgcGTAAACATGGATGAGGAGATGGACGACGACGAATAATGGCTTCTGCACAAGGGGTTAATGGTTGGGCACTTGACGCTATTCCCGGTTTTAAGGAAACCTACGCACATTAAGTAAGCGAAGTCGATACATTAAAGGAAAAAGACCATGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATcc                                                                                                                                                                                                      <  1:342149/227‑1 (MQ=255)
                                        gagaTTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAGAACATGAGGGGCAAATCACGGTCAATACATCATCAACGGCATTCCg                                                                                                                                                                                  >  1:624599/1‑229 (MQ=255)
                                             tGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCAGCCACAAAGTTATTGGCGGTAGACAATAAAGGGATTTAGGGGCATAACACGGGCAAACcacacataccgtcattccctcctg                                                                                                                                                                               >  1:290456/1‑207 (MQ=255)
                                                gcaccacaaaactttgtttttCCCAAAATGGTTAAGGGTTGTGCTCTAGATCTACATCCCTTTTTTATAAAACCCGACGACATTTTATTTCCAAAGGTCGAAACTTAAAAGGCAAAGGACCCTAAAAATTCTCTGAAAAATGCCGCCCCAAAGTTATTGGCGGTAGCCCATAAAGTGATTTAGGAGCATATCGCGGTCAATCCATCATCACCGGCATTCCGTCATGGcg                                                                                                                                                                          <  2:386779/212‑2 (MQ=255)
                                                ctaggaagaATAATGGATTTTCCACAAAGGGTTAATGGTTGGGTGATAAATGCTCACCCTTGTTTGCAGGAAACCTACGCCGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCa                                                                                                                                                                          <  2:160865/222‑1 (MQ=255)
                                                 tgagtctaataatGGGTTTCCCCAAAGGGGTTAATGGTGGGCGCGTAAGCGCTCATCCCTGTTTCCAAGAACCCGACGACGATTTAGTTGAAGAAGTCGAGAAATAAAGGGAAAAAGAGCCTGAAAAGTACAAGAGAAAAGCCGCCACAAAGGTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGAATTCCGTCATGGCa                                                                                                                                                                          <  1:350867/221‑1 (MQ=255)
                                                 cggcgaccaattagggttttCCCCAAAAGGTTGAATGGTTGCGCGTTATCGGGTCATCCCTCTTTCCAGGAAACTAACGGCGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCa                                                                                                                                                                          <  1:366521/213‑1 (MQ=255)
                                                      ccaataagGGATTTCCCAAAAGGGGTTAATGGTTGGGCGCTTAATGCTAACCCTTGTTTTCAGGAAACCTACGACGTTTAAGTACCCAAAGCCGAGAATTTAAAGGAAAAAGACCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGt                                                                                                                                                                   <  1:202446/227‑1 (MQ=255)
                                                       aaataatGGATTTCCCACAAAGGGTTAAGGGTTGGGCGCTATATGCTCATCCGTTTTTTCAGGAAACTTACGACGTTTTAGTTGCCGAAGCCGAGAAATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGtt                                                                                                                                                                  <  1:521007/228‑1 (MQ=255)
                                                             ttgATTTTCCAAAAAGGTTAAAGGGTTGGGCTCTATATGCTCATCCCTGTTTTCAGGAAACCTGCGACGTTTTAGTTGCCAAAGTCGAAAAATTAAAGGGAAAATACCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGct                                                                                                                                                             <  1:647480/226‑1 (MQ=255)
                                                                aTTTCCCAAAAGGGGTTAAGGGTTGGGCGCTATATGCTCATCCCTGTTTTAAGGAAACCTCCGACTCTTTAGTTGCCGAAGTCGAGACATAAAAGGAAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTgg                                                                                                                                                         <  1:476182/229‑1 (MQ=255)
                                                                  tttCCCAAAGGGGTTAATGGTGGGGGTATTGATGCACACCCTTTTTTTCAGAAAACCTACGCTGTTTTAGTTGCAGAAGACAGAACATAAAAGGAAAAAGACCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                       <  2:523544/229‑1 (MQ=255)
                                                                  tttCCAAAAAGGTTTAAGGGTTGGGTGTTTTATCTCTACCATTTGTTCAAGGAAACCTCCACGTATTTACGTCCAGAGGTGGAAACATTAAAGGAAAAAGACACTGAAAATTATCAGAAAAAGGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCTTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                       <  2:713590/229‑1 (MQ=255)
                                                                  ttcccACAAAGGGTTGATGGTTGGGCATTATATCCATACACTGTTTTGCAGGAAACCTACGACGCTTAATTTCCCGAAGTCGAGACATTAAAGGGAAAAGTTCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                       <  1:466420/226‑1 (MQ=255)
                                                                  ttcaccaaaggggTTAAGGGTGGGGCGTTAGAGGCACCTCACGTTTTACAGGAAAGCGAGGAGTATTTGTTTGCCAAAGTGGAGACATTAAAGGGAAAAGTTCCTAAAAATCATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCTCGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCTTTAGGCTCTGGGa                                                                                                                                                       <  2:704017/219‑1 (MQ=255)
                                                                  gtcaccaaaaGGGTTAATTGTGGGGCCTTTAGTCCCCAGCCCTTTTATAAGAAAACCAGCGACTGTTTGGTTGCCGAAGGCGACAAATAAAAGGAAAAAGATCGTAAAAATTATCAGAAAAAGTCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                       <  2:506308/222‑1 (MQ=255)
                                                                   tccctcaAAGGTGTAGTGGTTGGGTGTTACATTCTCAGCCTTGTTTCAAGAAATCCGTCGACCTTTTAGTTGCGCAAGTCGAGACATAAAAGGGAAAGGCTCCAGAAAATTATCAGAAAAAGTCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                       <  2:75508/223‑1 (MQ=255)
                                                                         aaaGGGTAAATGGTTGGGCGTTATATCGTAAGCCTTTTTTCAAGGAAACCAACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGAAAAATATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATaa                                                                                                                                                <  1:301614/229‑1 (MQ=255)
                                                                             ggTTAACTGGTTGGGTGTATTATGTTCATCCCTTTTTTCAGTAAACCTACGACGTTTAAGTTGCCGAAGTCGAGAATTTAAGGGAAATAGGTCCTGAAAATTAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATTGCAAGTCGTTAGGCTCTGGGAAAAATAAAGa                                                                                                                                             <  1:35232/229‑1 (MQ=255)
                                                                                  aaTGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAGTATTGAGGAGCATAATACGGTCAAACCATCATCACCGGCCTTCCCGCCTGGGCAGATGTTTGGGCTCTGGGCAAATAAAAACTggg                                                                                                                                        >  1:524614/1‑228 (MQ=255)
                                                                                   atgggTTGGGCGCTATATCCTCAACCGTTTTTCAAGGAAACCAACGCCGCTTTGGTGGCCGAATCCGAGACATTAAAGGAAAAAGATCCGAAAAATTACCAGACAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGt                                                                                                                                        <  2:624598/226‑1 (MQ=255)
                                                                                    gtGGTTGGGCACTATCAGCTACTCCCTGGTTGAAAGAACCATACGCAGCTTGAATTGCCGAAGGCGAGACAGTAAAGGGAAAAGGTCCTAAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTc                                                                                                                                      <  1:34831/228‑1 (MQ=255)
                                                                                      ggTTGGGCTCTATATGCTCATCCCTGTTTTCAGAAAACCTACGAGGTTTTAGTTGCAGAAGTCGGGACATTAAAGGAAAAAGACCTTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTAATGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCAc                                                                                                                                    <  1:728126/229‑1 (MQ=255)
                                                                                           aaggacccgatGCTTTTCCCTGTTTAGAGGAAACCAACCACGATTGATGTGGCAAAGTCGAGAAATTAAAGGAAAAAGCCCCAGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGt                                                                                                                                <  2:253052/220‑1 (MQ=255)
                                                                                            gccgttataTGTACACCCCTGTTTTAAGGAAACTTACGACATATTAGCTGCAGAAAGGGAGAACTAAAAGGGAAAAAGTTCAAAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTa                                                                                                                               <  2:752300/224‑1 (MQ=255)
                                                                                                       gccaaaccttGTTTTAAGGAAACTTACGACGCTTAAGTTCCAGAATCCGAGACTTTAAAGGGAAAAGATCCTGAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                     <  1:219269/220‑1 (MQ=255)
                                                                                                       gCTCACCCGTGTTTTAAGGAACCCTTCGCTGATTAAGTTGCCGAAGGAGAGCCATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGTAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                     <  1:279643/229‑1 (MQ=255)
                                                                                                        gacATCACTGTTTTGAGGAAACCTACGACCTTTAGGTTGCCAAGGTCGAGACATTAAGGGGAAAATACCCTGAAAATTTTCAGAAAAAGTCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                     <  1:174348/226‑1 (MQ=255)
                                                                                                            tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAAGAATTCAGAGGCTTAAACGGTGAATCCATCATCACCCGCAATACGTCAAGGCAAAGCGCTTGGCTCCGGGTAAAAAAAAAAATCGTCACGGGGAAAAATTTGGGATGGt                                                                                                                >  2:219269/1‑229 (MQ=255)
                                                                                                            tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGGTGGCCATAAAATGATTTAGGAGCATATCACGGTCAATCCATAATCACCGGCAATCCGTACCGGGAAGTCGGTAAGGCATGGGGAAAATTAAAAACGGTTACGGGGAAAAATTgggtgtggt                                                                                                                >  1:108607/1‑222 (MQ=255)
                                                                                                              ccTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAAGATTGAGGAGCAAATCACGGTCAAACCATCATCACCGGCATTCCCTCATGGCAAGGCGTTAGGCTCTGGGGAAAATAAAAACTGGTCACGGGGCAAATTTTGTGCTTGt                                                                                                                >  1:529410/1‑227 (MQ=255)
                                                                                                                tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAATGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTTCGTCCTGGCAAGTCGTTAGGCCCTGGGGAAAATAAAAAATGGGCACGGGGAAAAATTTGGGCTgtgtct                                                                                                             >  1:253052/1‑223 (MQ=255)
                                                                                                                tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAAGGTTGAGGAGCATATAACAGGCAATCACTCATCACCCGCCTTCCCTCATTGGAAGACGTTTGGGGCTGGGGAAAAAAAAAACTAGTACCCGGGGAAAATTTGTGATGGGCGtt                                                                                                            >  2:569612/1‑229 (MQ=255)
                                                                                                                tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCCTAAAGAGATTTAGGGGCATATCACGCGCAATCCAGCATCCCCCCCAATCCGTCATTGCAAGACGTTAGGGCCTGGGTAAAAAAAAAACAGGGCCCGGGGaaaaatagtggttggtgggt                                                                                                            >  2:532081/1‑213 (MQ=255)
                                                                                                                  ttttCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGGCATAAAATGATTGAGGAGCATAACACGCTCAATACATCATAACCGGCATTCCGTCATAGGAAACCGTTAGGCTCTGGGGAAAAAAAAAACTGGGCACCGGGAACAATTTGTGCTgggtgctat                                                                                                          >  2:711865/1‑222 (MQ=255)
                                                                                                                  ttttAAGGAAACCTACGACGTTTTAGTTGCGGAAGTCGGAACATTAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                                                                            <  2:246773/135‑1 (MQ=255)
                                                                                                                     tCAGGAAACCTACGACGCTTTAGTCGCCGAAGTCGAGACATTAAAGGGAAAATATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGt                                                                                                       <  1:507552/229‑1 (MQ=255)
                                                                                                                       taaGGAACCCTAGTACGTTTAAGTTGCCAAATGCAGGAAATAAAGGGGAAAAGACCCTAAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGt                                                                                                       <  1:227779/227‑1 (MQ=255)
                                                                                                                       taaGGAAAACTACCGCGTTTTAGGTGCCGAAGGCGAAAAATTAAGGGGAAAAGGTTCTGAAAAATAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGt                                                                                                       <  2:290456/227‑1 (MQ=255)
                                                                                                                                       cgttttAGTTGCCGAAGTCGAGAAATTAAAGGGAAAAGACCCGGAAAATTATCAGAGAAAAGCCGGCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTta                                                                                         <  1:77360/225‑1 (MQ=255)
                                                                                                                                         tgtttatgtggcGAAGTCGAAAAGTAAAAGGGAAAAGACCCTGAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCTGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTta                                                                                         <  1:110922/216‑1 (MQ=255)
                                                                                                                                                      aagaggagaAATAAAAGGAAAAAGATCCGGAAAATAATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGaa                                                                         <  2:740378/224‑1 (MQ=255)
                                                                                                                                                         tCGAGAAATTAAAGGAAAAAGATCCTGAAATTTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGt                                                                      <  1:468815/229‑1 (MQ=255)
                                                                                                                                                           agaaCATTAAGGGAAAAAGATCGTGAAAATTAACAGGAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                                                                    <  1:711866/226‑1 (MQ=255)
                                                                                                                                                            agaCATTAAAGTGAAAAGACGCTGAAAATTAAAAGAGAAAACCCGCCCCAAAGTTATTGGCGGTAGTCCATAAAGTGAATGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                                                                    <  1:723252/228‑1 (MQ=255)
                                                                                                                                                                 taaaaGGGAAAATATCCAGAAAATGAAAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGCCCATAAAGGGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCAtt                                                               <  2:529410/226‑1 (MQ=255)
                                                                                                                                                                 taaaaGGAAAAAGACCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCAtt                                                               <  1:636275/226‑1 (MQ=255)
                                                                                                                                                                  taaagaggaaaaagccgaaaaaaTTAGCAGAGAAAAACCGCCACAAAGTTATTGGGGGGAAGCCATAAAGGGATTGAGGAGCATATCACGGGCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCGCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTc                                                              <  2:576610/210‑1 (MQ=255)
                                                                                                                                                                  tAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTTGTGCATCAAAGGATTTAGGGGGCAATCACGGACAAATCATCCATCCCGGCGTTCCGTCTTTGCAAGTCCGTAGGCTTGGGGAAAAAACAAAAATGGGCACGGGGAAAAATTTGTGGACGTCGTTGACGTCCATTATTTTTTTTTAGTTGAAAAAGGAAAAgaaccaagta                                                              >  2:495101/1‑220 (MQ=255)
                                                                                                                                                                          aaaGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGACCCTAAAATAGTTGGGGAGCAAATCACGGTCAATCCATCCTCACCGGCATTCCTTCATGGCGAGTAGTTAGGCTCTGGGCAACATCAAGACTGGGGACGGGGAAACATTTGTGCTGGGCGTTATCGTTTCTTCTTTCGTTTTAATGAAAAAAAGAAAAGTATCATTCTGggatggg                                                     >  1:421645/1‑228 (MQ=255)
                                                                                                                                                                                  cTGAAAATTAGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                     <  1:65512/159‑1 (MQ=255)
                                                                                                                                                                                            tgagagAAAAGCCGCCACAAATGTATTGGCGGTCGTCCAGAAAGTGATTGAGGGTGATATCACGGTCAATCGATCATCACCGGCATCCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACAct                                    <  1:359966/227‑1 (MQ=255)
                                                                                                                                                                                             acgagaAAAGGCCGCCACAATGGTATTCGCGGTACTCCATAAAGTGATGGAGGAGCATATCACTGTCAATACCTCATAACCGGCCTTCGGTCAGGAAAAGTTGGTAGGCTCTGGCAAAAATAAAGGCTGGGCTGGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCGTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAacaaa                                    <  2:579833/227‑3 (MQ=255)
                                                                                                                                                                                               gaaaaaaGCCGCCACAAAGGTATTGGCGGTATTCCTAAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTg                                 <  2:571664/226‑1 (MQ=255)
                                                                                                                                                                                                   aaaggccgcccacacgTTATTTGCGGTAGTCCAGAAAGTGATTGAGGCGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTGAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTgcgc                              <  1:724928/214‑1 (MQ=255)
                                                                                                                                                                                                                     tttGGTGGTTGTCCATAAACTTATTTGGGCGTATCTCATGGTCAATCCATCATCACCGGCAGTCCGTCAGGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATGTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGACCGATGAAAACACTCTGCGCACCTCCGGTAAAAAAACa            <  2:490594/228‑1 (MQ=255)
                                                                                                                                                                                                                     tttGGCTGTAGTCCATAAAAGTATTGGCGAGCATACCACGGTCAATCCATCAACACCGGCATTCCGTCAGGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCATGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTGCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa            <  2:108607/228‑1 (MQ=255)
                                                                                                                                                                                                                     cttGGCGGTGGTACAAAAAGGGATGGAGGAGAATATCACGGGCAAACCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGTTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa            <  2:709899/228‑1 (MQ=255)
                                                                                                                                                                                                                      ttGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa            <  1:445417/228‑1 (MQ=255)
                                                                                                                                                                                                                      ttGGCGGTAGGCCATAAAGGGATTGAGGGGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAg           <  1:451466/229‑1 (MQ=255)
                                                                                                                                                                                                                      atGGCGGTAGTCCATAAAGGTATTGAGGAGCATATCACGGTCCATCCATCAGCACCGGCATGCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCCCCTACGGTAAAAAAACAg           <  1:725052/228‑1 (MQ=255)
                                                                                                                                                                                                                       tGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCAGCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa          <  1:651521/229‑1 (MQ=255)
                                                                                                                                                                                                                       tGGCGGTAGTCCAAAAAATTAGTGTGGGTGATAACACGGTCAATCCAACAGCCCCGGCATTCCGTGATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa          <  2:62319/229‑1 (MQ=255)
                                                                                                                                                                                                                       gggCGGTAGGCCATGAAGTGGGTGAGGAGCATTTCACGGTAAATCCACCATCCGCGGCATTCCGTCATGGCAAGTCGTTTGGCTCTGGGAAAAATAAAGACTGTGCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa          <  2:116363/228‑1 (MQ=255)
                                                                                                                                                                                                                          gcggttcccgtAAAGTTGTTGAGGAGCATTTCACGGTCAAGCCCTCAACGACGGCTCTACTACATGGCCACTCTCTAGTCTCTGGTAAAAATAAAGACTGTTCGCGTGTAAAATTTGGTGCTGGTCGTTAGCGTCTCTTCTTGCGTGATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGACCGATTAAAACACTCTGCCCACCTACGGTAAAAAAACAGATGc       <  2:525557/219‑1 (MQ=255)
                                                                                                                                                                                                                          cGGTAGTCCATAAAGTGATTGATGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGc       <  1:676909/229‑1 (MQ=255)
                                                                                                                                                                                                                          cGGTAGTCCATAAAGGGATTGAGGAGCATATTACGGTAAAGCCATCAACACCGGCATCCCGGCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGc       <  2:44349/229‑1 (MQ=255)
                                                                                                                                                                                                                            gttgTCCATAAAGTTATTGGGGAGCATAATGCGGTCACTCCATCAACACCGGCATTTCGTCCTGGCAAGGTGTTAGGCTCTGGGAAAAATAAAGACTTGTCCCGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTGTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGCAATCACGCTGCGCACCTACGGTAAAAAAACAGATGCCt     <  2:512527/226‑1 (MQ=255)
                                                                                                                                                                                                                               gTCCATAAAGTGATTTAGGAGCATATAACGGTTCATCCATCAAACCCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTTCTGGTCGGTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCtata  <  2:47639/229‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
AAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATA  >  NC_000913/3277275‑3277716

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