Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,615 C→T intergenic (‑93/‑118) ettA ← / → slt energy‑dependent translational throttle protein EttA/soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,6150CT95.0% 177.9 / ‑1.6 61intergenic (‑93/‑118)ettA/sltenergy‑dependent translational throttle protein EttA/soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (31/26);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAG  >  NC_000913/4630365‑4630858
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
gCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGCCCCATGACACC‑AATTTTTGCCCCAGTGACGAAACTCAGAGAGATGTTTTTCAACATATGACGTTTCGGCGGAACAACTTTGCCGACCCGATGCATGGTCTAAACGATTTGAGCCTCGTTGGCCTTCGCCTCAATTTTTATCGTGATGCTGAGTTTTCAAAGGCGAAGTGTCTCCTTTTTTCCCGCCTAATCCCAGTTGGTCGAT‑ACCCt                                                                                                                                                                                                                                                     >  2:294110/1‑251 (MQ=255)
        gcaccgtggggtagaattACCCGAGCCATTCAGACCCATGA‑ACCTAATTTTTGCCCCAGGGAAGAATCTCACAGATATGTTTTTCAAAATATGTCGGTTCGGCGGAACAATTTTTCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGCCTTCGCCTCTATGTTTATATTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCTTAATCCCAGTTGGTCGATCACCTTCGCagta                                                                                                                                                                                                                                              <  2:604898/235‑1 (MQ=255)
          gcagcaGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagtaaa                                                                                                                                                                                                                                            <  1:329570/251‑1 (MQ=255)
          gcagcaGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagtaa                                                                                                                                                                                                                                             >  1:589543/1‑250 (MQ=255)
                        ttaccccagcagttCAGACCCAGGACACC‑CATTTTTGCCCCAGGGAAGCAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACCACTTTGCCGACACGATGCATGGTATCAACGAAATGAGCCACGTTGGACTTCGCCTCTTAGTTTATCGTGATGAGGAGTTTTCAAAGGCGCAGTGTAGCCTTTTTCCCCGCCTGATCCCAGTTGGTCGATCAACTTCGCAGTAAATGTaaaaaatggac                                                                                                                                                                                                                              >  2:283220/13‑246 (MQ=255)
                          aCCCGCGCCGTTCAGACCCAGGAAACC‑AATTTTTGCCCCAGGGAAGAAACTTCGAGAGATGTTTTTAAAAATATGACTGTTCGGCGGAACAACCTGGCCGACCCGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGAGGAGTTTGCAAAGGCGAATTGTAGCATTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTcc                                                                                                                                                                                                                             <  2:192713/250‑1 (MQ=255)
                            ccgggccgttCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAAAAACTCAGAGAGATGTTTTGCAAAATATGACGTTTCGGCGGAACAACTTTGCGGACACGATGCATGGTATAAACGAATTGAGCCTCGTTGGATTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGTCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTa                                                                                                                                                                                                                          <  2:584174/243‑1 (MQ=255)
                                 ccgttCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTg                                                                                                                                                                                                                    <  1:60363/248‑1 (MQ=255)
                                  cgttCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTGCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTgg                                                                                                                                                                                                                   <  1:213010/249‑1 (MQ=255)
                                      aaGACCCAGGACACC‑AATTTTTGCCCCAGGGAAAAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTGCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTa                                                                                                                                                                                                               <  2:640634/250‑1 (MQ=255)
                                       aGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTa                                                                                                                                                                                                               >  1:426374/1‑250 (MQ=255)
                                           ccAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACgg                                                                                                                                                                                                          <  1:303735/251‑1 (MQ=255)
                                           ccAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACgg                                                                                                                                                                                                          <  1:257040/251‑1 (MQ=255)
                                                 cacC‑AATTTTTGCCCCAGGGAAGAAACTAAGAGAGATGTTTTTCACAATATGACGTTTCGGCGGAACAACTTTGCCGACACGTTGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTATCCGTAAC‑GTGGCGTAAACGGCAATg                                                                                                                                                                                                     >  2:465817/1‑250 (MQ=255)
                                                   cc‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAGAAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGACGAGTTTTCAAAGGCGAAGCATAGCCTGTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGAAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACt                                                                                                                                                                                                  >  2:533939/1‑251 (MQ=255)
                                                              ccccAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCAt                                                                                                                                                                                         >  1:257683/1‑250 (MQ=255)
                                                                      aagaaACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTAtt                                                                                                                                                                                <  1:465817/251‑1 (MQ=255)
                                                                         aaaCTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑Gt                                                                                                                                                                                                                     >  1:661780/1‑211 (MQ=255)
                                                                              cAGAGAGATGTTTTTAAAAAAATGACTTTTTGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GGGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCa                                                                                                                                                                        <  2:60789/251‑1 (MQ=255)
                                                                                 agagaTTTTTTTCAAAATATTACGTTTCGGCGGAACACCTTTGCCTACACGAGGCTTGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCACAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGa                                                                                                                                                                     <  1:20182/251‑1 (MQ=255)
                                                                                    gaTGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGc                                                                                                                                                                  >  2:171294/1‑251 (MQ=255)
                                                                                    aaTGTTTTTAAAAATAGGACGTTGCGGCGGAACATCTTTGCCGATACGATGCATGGTATAAACGAATTGAGCCCTTTTGTACTTCGCCTCTATGTGTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTATCATAAATCTATTACGCGGCATGACGc                                                                                                                                                                  <  2:457484/250‑1 (MQ=255)
                                                                                         ttttCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCAt                                                                                                                                                             >  1:321099/1‑251 (MQ=255)
                                                                                             cAAAATATGACGTTTCGACGGAACAACTTTGCCGACACGATGCATGGAATAAACGACTTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCAATAATCTATTACGTGGCATGACgaggaattga                                                                                                                                                          >  2:569744/1‑241 (MQ=255)
                                                                                                       cGTTTCGACGGAACAGCTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGAATTCGCCTCTATGATTATCGTGATGATGAGATTTCACAGGCGAAGTGTAGCCTTTTTCACCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAACGTGTCCGTAAC‑GTGGCGTAAACGGCAATGCATGGTTAGCATAAATCTTTTACGCGGCATGACGCTGCATTGATGTATTTacac                                                                                                                                               >  2:270830/1‑251 (MQ=255)
                                                                                                               cggAACAACTTTGCCGACACGATGCATGGGATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGAAGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGAAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGAAGACGCATTGATGTCTTTACACTTAGAgg                                                                                                                                        >  2:464466/1‑250 (MQ=255)
                                                                                                                ggAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGAGTAGCCTTTTTCCCCGCCTAATCACAGTTGGCCGATCACCTTCGCCGTAAAAATAAAAAAGTGTCCG‑AACAGTGGCGTAAACGACAATAACTGGATAGCATAAATCTATTACGCGGCATGACGCTGCCTTGCAGTATTTACACTTcgagtat                                                                                                                                      >  2:510282/1‑248 (MQ=255)
                                                                                                                ggAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGAt                                                                                                                                      >  1:260371/1‑251 (MQ=255)
                                                                                                                 gAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTg                                                                                                                                                                                                 <  2:11742/191‑1 (MQ=255)
                                                                                                                 gAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTg                                                                                                                                                                                                 >  1:11742/1‑191 (MQ=255)
                                                                                                                  aacaacTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATgc                                                                                                                                    <  1:387470/251‑1 (MQ=255)
                                                                                                                   acatcTTTGCCGACACGTTGCATTGTATAATCGATTTTAGCCCTGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCTATCCCCTTCTCAGTAAAAGTAAAAAAGTGTCCGTa                                                                                                                                                                                                                          <  2:580997/161‑1 (MQ=255)
                                                                                                                   acaacTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTa                                                                                                                                                                                                                          >  1:580997/1‑165 (MQ=255)
                                                                                                                     aacTTTGCCGACACGTTGCATGGTATAAACGAATTGAGCCATGTTGGACTTCGCCTCTATGTTTAACGTGTTGATTAGTTTTCAAAGTCGAAGTGTAGCTTTTTTCCCCGACTAATCCCAGTTGGTCGATCACCGTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCt                                                                                                                                 <  2:207128/251‑1 (MQ=255)
                                                                                                                               acacGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGaaaa                                                                                                                        >  1:88572/1‑250 (MQ=255)
                                                                                                                                    aTGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTATCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGCGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGACTTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCa                                                                                                                   >  1:318266/1‑250 (MQ=255)
                                                                                                                                    aTGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCAACGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCa                                                                                                                   >  1:477726/1‑250 (MQ=255)
                                                                                                                                               aaaCGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGGGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCt                                                                                                       <  2:257683/251‑1 (MQ=255)
                                                                                                                                               aaaCGAATTGAGCCACGTTGGACTTCGCCGCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTAcc                                                                                                        <  1:542660/250‑1 (MQ=255)
                                                                                                                                                  cgcaTTGAGCCGCGATGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAACGGAGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGAAAAAAAGTGTCCGTAAC‑TTGGCGAAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGCAGCTGCATTGATGTATTTACACTTAGCGAATGCGCTTGTGGATAAAGCCAAACAAGTTACCTggc                                                                                                    >  2:560023/4‑251 (MQ=255)
                                                                                                                                                             aCGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCt                                                                                         >  1:122978/1‑251 (MQ=255)
                                                                                                                                                                   gACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGgt                                                                                   <  1:388700/251‑1 (MQ=255)
                                                                                                                                                                         cccTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGtgtctgt                                                                             <  1:626187/250‑1 (MQ=255)
                                                                                                                                                                             ctATGTTTATCGTGAAGAGGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGGCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTgtctgtct                                                                           <  2:416388/249‑1 (MQ=255)
                                                                                                                                                                               atttttATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGtctgtctg                                                                          <  2:260371/245‑1 (MQ=255)
                                                                                                                                                                               aTGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTctgctg                                                                       >  1:488307/1‑251 (MQ=255)
                                                                                                                                                                                             gaTGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGt                                                         >  1:152116/2‑251 (MQ=255)
                                                                                                                                                                                                   ttttCAAAGGCGCAGTGTAGCCTTTTTCCCCGCCTAATCCCAATTGTTCGATCACCTTCGCAGTAAAAGTAAAACAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTATGCTGACGGTCAGCAGCGTGGcgcg                                                   >  2:532176/1‑251 (MQ=255)
                                                                                                                                                                                                      tCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGc                                                >  1:248744/1‑251 (MQ=255)
                                                                                                                                                                                                         aaGGCGAAGTGTAGCATTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGCCGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAACAAGCCAAACAAGTTACCTGGCGGCTGTTGGATGCCGGTGTCTGTCTGCTGACGGTCAGCAACATGGCGCGATCCg                                              >  2:660107/1‑250 (MQ=255)
                                                                                                                                                                                                              gcaGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCAACTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTAAGCAGCGTAGCGCGGGCCGACTCAc                                        >  2:29728/3‑251 (MQ=255)
                                                                                                                                                                                                                 gtgtAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTCCCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTgg                                     >  2:461966/1‑251 (MQ=255)
                                                                                                                                                                                                                 gtgtAGCCTCTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTgg                                     >  1:450062/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGtt                   >  1:295135/1‑250 (MQ=255)
                                                                                                                                                                                                                                      aaTCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCAc                                        <  1:651138/227‑1 (MQ=255)
                                                                                                                                                                                                                                      aaTCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCAc                                        >  2:651138/1‑227 (MQ=255)
                                                                                                                                                                                                                                      aaTCCCAGTGGGGCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTATC‑GTGTCGTAAACGGCAATTACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCAATGATGTATTTACACTTAGAGGATTCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTGGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCTCTGGATGATCAGCGTAGCCGTTat                 <  2:605671/250‑2 (MQ=255)
                                                                                                                                                                                                                                                   cGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCa   <  1:171294/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    gATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGc    >  1:181659/1‑249 (MQ=255)
                                                                                                                                                                                                                                                    gATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAg  <  1:360874/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    gATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAg  <  1:283220/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAAC‑GTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAG  >  NC_000913/4630365‑4630858

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