Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,277 G→A L171L (CTC→CTT)  ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,2770GA100.0% 175.1 / NA 54L171L (CTC→CTT) ytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/39);  total (15/39)

GGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCAT  >  NC_000913/4624031‑4624512
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
gggCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑GGCAGCGTGCGGGTCAGCACTACGCCGCTTGCCAGTGAAagcc                                                                                                                                                                                                                                           >  1:432368/1‑250 (MQ=255)
   catctGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGAGGGGACGTTTACCCTGT‑AGCAGCGTGCGGGGCAGACCTACGCCGATTTCCAGTGAAagcagc                                                                                                                                                                                                                                         >  1:173649/5‑250 (MQ=255)
           aaaaaaaACACTCCCCCTCCGCCAGGTTAAACCACGGGTCTTAAGATGCAGAGATTAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTGGCCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCATCAACATCAgg                                                                                                                                                                                                                                <  2:233386/250‑1 (MQ=255)
            aaaaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGc                                                                                                                                                                                                                               <  1:651444/251‑1 (MQ=255)
             aaaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGGCAGCACTACGCCGATTGGCAGTGAAAGCAACAACATCAGGCg                                                                                                                                                                                                                              >  1:483043/1‑251 (MQ=255)
             aaaaTACACGCCTCCCCCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCg                                                                                                                                                                                                                              <  1:574751/251‑1 (MQ=255)
              aaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGt                                                                                                                                                                                                                             <  1:489380/251‑1 (MQ=255)
               aaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTa                                                                                                                                                                                                                            <  1:444450/251‑1 (MQ=255)
                          ccACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtg                                                                                                                                                                                                                 <  1:547248/251‑1 (MQ=255)
                             ccGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtggtg                                                                                                                                                                                                              <  1:94511/251‑1 (MQ=255)
                               gCCAGGTTACACCACGGGTCGTAAGAGGCCGAGATGAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTGGTAGGAATGGGTCACACTTTGTTCTCGCTTTCTTCCTCTTCCGGCACCGTTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGtt                                                                                                                                                                                                            <  2:464672/251‑1 (MQ=255)
                                  aGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATc                                                                                                                                                                                                        <  1:107205/251‑1 (MQ=255)
                                            aCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTgg                                                                                                                                                                                              <  1:595700/251‑1 (MQ=255)
                                             cggtttgTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTggg                                                                                                                                                                                             <  2:487277/245‑1 (MQ=255)
                                              gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGc                                                                                                                                                                                            <  1:617043/251‑1 (MQ=255)
                                                  cGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCg                                                                                                                                                                                        <  2:312729/251‑1 (MQ=255)
                                                   gTAAGAGCTAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTGTCTTGTAGTAATGGGTCACTCTTTTTGCGCGCTTTCTTCCTCTTCCGGCACCGGTTTTCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCgg                                                                                                                                                                                       <  2:14083/251‑1 (MQ=255)
                                                    tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAAGGAGGGCGATTGCTGCCCGCGGGGGCGGTTTCCCTTTCACAACCGTGGGGGTCGCCAC‑ACGCCGGTTGCCAGGGAAAGCCGCAACAACAGGGGGAAAATGTTTATGGTGTT‑ATCCCCCTGTTGGGCTTCggg                                                                                                                                                                                      >  2:680848/1‑250 (MQ=255)
                                                     aaGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTg                                                                                                                                                                                     <  1:199025/251‑1 (MQ=255)
                                                                aTGAGCAGGCTTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATg                                                                                                                                                                          <  2:243935/251‑1 (MQ=255)
                                                                        gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAACAGCAACAACAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGGGCGTATgggggtggg                                                                                                                                                                  >  1:637140/1‑245 (MQ=255)
                                                                               gtgGACATAAAGAGTTCCTTTCTTGTAGTAATGGGGCACTCGTTTTTCTCGCTTTCTTCCTCTTCCGGCTCCGGGTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAg                                                                                                                                                           <  2:251036/251‑1 (MQ=255)
                                                                                    cATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagca                                                                                                                                                                                                                                           <  1:317838/167‑1 (MQ=255)
                                                                                     aTAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAg                                                                                                                                                     <  2:46955/251‑1 (MQ=255)
                                                                                                  ttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCgg                                                                                                                                         <  1:545272/250‑1 (MQ=255)
                                                                                                       gtagtaATTGGTCACTCTTTTTTCTCGCTTTCTTACCCTTCCGGCACCGGTTTGCTGGCGGTGAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtt                                                                                                                                   <  2:148634/251‑1 (MQ=255)
                                                                                                                ggTCACTCTGTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccggc                                                                                                                          <  1:680848/251‑7 (MQ=255)
                                                                                                                   cgctctTTTTTCGCGCTTTCTTCCTCTTCCGGCTCCGGTTTGCTGGCGGTTAACAGGAAGGGCGAGTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAAt                                                                                                                      <  2:222499/249‑1 (MQ=255)
                                                                                                                      tctTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGACGATTGCCAGTGAAAGCCGCAAAATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCCGGTAGCCGAGCGGTCCGGTTTTCCCGG‑CAATTgg                                                                                                                   >  1:405510/1‑251 (MQ=255)
                                                                                                                       ctTTTTGCTCGCTTTCTTCCTCTTCCGGCACCGGGTTGCGGGCGGGTAACAGTAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAAtta                                                                                                                    <  2:19881/249‑2 (MQ=255)
                                                                                                                                  cTTTCTTCCTCTTCCGGCACCGATTTGCTGGCGGTTAAAAGGGAGGGCGATTACTTCCAGCGGGGACGTTTACCCTGT‑AGCAGCGTACGGGTCAGCCCTACGCCGCTTGCCAGAGCAAGCAGCACCATCAGGCGGAAAATGTTGGGGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGGGCGTGTTGGGGTCGCGCGTCAAGCGCCGATAACCGAGCGGTCC‑TTTTTTCCCGTCAATTGGTTCGGCaaacct                                                                                                       >  1:172941/1‑247 (MQ=255)
                                                                                                                                           cctgtcGGCACCGGTTTGCTGGCGTTTAACAGGAAGGGCGATTGCTGCCAGCGTGTACGCTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCTATTTCCAGTGAAAGCAGCAACATCAGGCGTAAAAGGTTAGTGGTGTT‑ATCCACCTGTGGGGCTTCGTTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCTGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTaaa                                                                                              <  1:440958/246‑1 (MQ=255)
                                                                                                                                               ccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAg                                                                                          <  1:180873/251‑1 (MQ=255)
                                                                                                                                                cGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagca                                                                                                                                                                                                                                           <  1:322652/107‑1 (MQ=255)
                                                                                                                                                    cccGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑AGCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAgccgcc                                                                                     <  1:177685/250‑1 (MQ=255)
                                                                                                                                                                 gggTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGG‑AGCAGCGTTCGGGGCAGCACTACGCCGATTGCCAGTGAAAGCAGTAACAGCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCg                                                                        <  2:405510/250‑1 (MQ=255)
                                                                                                                                                                          ggAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGcc                                                               <  1:483998/251‑1 (MQ=255)
                                                                                                                                                                             aGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGc                                                            >  1:150537/1‑251 (MQ=255)
                                                                                                                                                                                      gctgcCAGCGGGTACGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTcgcg                                                    <  1:575827/250‑1 (MQ=255)
                                                                                                                                                                                                ggtacGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGCTTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGc                                         >  1:600191/5‑251 (MQ=255)
                                                                                                                                                                                                  tacGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCtt                                       >  1:114486/3‑251 (MQ=255)
                                                                                                                                                                                                     gTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                            <  2:227223/164‑5 (MQ=255)
                                                                                                                                                                                                     gTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                            >  1:227223/1‑160 (MQ=255)
                                                                                                                                                                                                       ttACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGGGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGcc                                  >  1:195227/1‑251 (MQ=255)
                                                                                                                                                                                                              gT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAAc                             <  1:477607/249‑1 (MQ=255)
                                                                                                                                                                                                                 agcggcGTGCGGGTCAGCACTTAGCCACTTGCAAGGGAAAGCAGCCACAGCAGGGGTAAAATTGTAG‑GGTGTTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGTTGTCGAGCGTCAGGCGCGGAGAGCCGAGCGGTCCGTTTTTCCCGG‑CAGTTGGTGCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCCGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                         <  2:173649/247‑4 (MQ=255)
                                                                                                                                                                                                                 agcagcTTGCGGGTCAGCACTACTCCGATTGCCAGTGAAAGCATCAACATCAGGCGTAAAATGTTTGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                         <  2:637140/251‑4 (MQ=255)
                                                                                                                                                                                                                     gcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                     <  1:139829/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                     >  1:42438/1‑250 (MQ=255)
                                                                                                                                                                                                                      cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                     >  1:348039/1‑250 (MQ=255)
                                                                                                                                                                                                                            gtgcatcactacggcGATTGCCATTGAAAGCAGCATCATCATGCGTAAAATGTTAGTGGTGTT‑ATCCGCCTGTTGTGCTACTGTGGCGGTCGTATTGGTGTCGAGCGGCAGGCGCAGATAGCCGTGCGGTCCTTTTTTCCCGG‑CAATTGTTTCGACAATCTGCTGGTTAAAATAGCCGCCCGTTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCAGCTCGACGCTTTCGCCAGAACGGGCGATAAGATCGc              <  2:342863/237‑1 (MQ=255)
                                                                                                                                                                                                                                aGCATTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCTTAAAATGTTAGTGGTGTT‑ATCCCCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCTAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCTCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTg          <  1:272932/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ctacgccggttgctcgTGAAAGCAGCAACAGCAGGCGTAAAATGTTAGTGGGGGT‑AGCATCCTGTTGGGCTTCGTTGGCTGGCGTATGGGTGTTGCGCGTCAGGCGCGGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAACCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGtt        <  2:375496/234‑1 (MQ=255)
                                                                                                                                                                                                                                        gccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCAt  >  2:567033/2‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
GGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTT‑ATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCAT  >  NC_000913/4624031‑4624512

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