Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,162 C→T H255Y (CAC→TAC)  spoT → bifunctional (p)ppGpp synthase/hydrolase SpoT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,1620CT100.0% 128.4 / NA 38H255Y (CAC→TAC) spoTbifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/17);  total (21/17)

GGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGG  >  NC_000913/3822919‑3823399
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
ggTATCCACCACATTAAAACCGAACTCGAATAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAAGACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATGTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGAt                                                                                                                                                                                                                                        <  1:166047/251‑1 (MQ=255)
      ccccacATTAAAACCGAACGCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGa                                                                                                                                                                                                                                  <  2:246248/249‑1 (MQ=255)
      caccacATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGAGTTTTTACTCGATCATGGa                                                                                                                                                                                                                                  <  2:246255/251‑1 (MQ=255)
                           gAAGAGCTGGGTTTTGAGGCGCTGTAGCCCAACCGTTAGCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCGTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGGCGCGAGAAGCAGCTTCATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTGTGTACTCGATCATGGACATCTACGCTTTCCGCGTGAt                                                                                                                                                                                                             <  2:188185/251‑1 (MQ=255)
                            aaGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATc                                                                                                                                                                                                            >  2:193123/1‑251 (MQ=255)
                                gCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCa                                                                                                                                                                                                        <  2:343237/251‑1 (MQ=255)
                                    ggTTTTGAGGCGCTGGATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAAGCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCAGCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGTGTTTTTACTCGATCATGGACATCTACGCGTTCCGCGTTATCGTCAATGa                                                                                                                                                                                                    <  2:212318/251‑1 (MQ=255)
                                     gTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTCATCCAAGAAGTCGCGAAAGCCGCGCGCGCCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAACGGCGCTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGCGCAGCGTTTTTACTCGATAATGCACATCTACGCTTTCCGCGTGATCGTCAATGAt                                                                                                                                                                                                   >  1:230867/1‑251 (MQ=255)
                                                cTGTCTCCCAACCGTTATCGCGTAATCCAAGAAGTGGTGTAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGAGTTTGCAGTAAGCGGGAATACCGTGCCGCGCCAGTGGGCGCGAGAAGCATCTTTATTCGATTTACTCCAAAGTGGGGCTCAAAGAGCAGCGTTTTTACTCGCTCATGTACATCTAGGCTTTCCGCGTGCTCGTCCATAATTTTGAAACCTg                                                                                                                                                                                        >  2:197732/1‑251 (MQ=255)
                                                      cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATc                                                                                                                                                                                   <  1:193123/250‑1 (MQ=255)
                                                        cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcgcg                                                                                                                                                                                <  1:109487/251‑1 (MQ=255)
                                                                       aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCATAAGATGCTGTCTGAAATCGAAGGGCGTGTGCAGGAAGCGGGAATACCGTGCCTCGTCAGTGGTCGCGAGAGGCATCTTTATTGGATTTACTGCAAAGTGGTGCTCAAAGAGCGGCGTTTTTACTCGATCATGGACATCTACGCTGTCCGCATGATCGTGAATGATTCTGACATGTGTTATCTCGTGCTGGGCCAGATg                                                                                                                                                                  <  2:216937/250‑1 (MQ=255)
                                                                             agaagTGGTGAACGCCGCGAGCGGCAACCGTAAAGAGCTGCACCAGAAGATTATTTCTGAAATCGAAGGGCGTTTGCAGGACGAAGGACTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGcc                                                                                                                                                           >  2:213829/1‑251 (MQ=255)
                                                                              gaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCt                                                                                                                                                          <  1:201673/251‑1 (MQ=255)
                                                                                agTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATTCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGACGCATCTTTATTCGATTTACTGCGAAATGGCGATCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGGTTCTGACACCTGTTATCGCGTGCTGGGACGGATGCACAGCCTGt                                                                                                                                                        >  1:31732/1‑251 (MQ=255)
                                                                                  tggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAc                                                                                                                                                      >  1:103471/1‑251 (MQ=255)
                                                                                           ccGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGt                                                                                                                                             >  1:73930/1‑251 (MQ=255)
                                                                                            cgcgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTc                                                                                                                                            <  1:20243/251‑1 (MQ=255)
                                                                                              cgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGAGGCACAGCCTGTACAAGCCGCGTCCg                                                                                                                                          <  2:99606/251‑1 (MQ=255)
                                                                                                  cgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCcgc                                                                                                                                               >  1:44948/1‑242 (MQ=255)
                                                                                                  cgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCAGGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCcgc                                                                                                                                               <  2:44948/242‑1 (MQ=255)
                                                                                                          gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGCAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaa                                                                                                                               >  1:126221/1‑250 (MQ=255)
                                                                                                          gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                                                                              <  2:38743/251‑1 (MQ=255)
                                                                                                          gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaa                                                                                                                               >  2:179279/1‑250 (MQ=255)
                                                                                                                  atgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                                                                      >  1:129038/1‑251 (MQ=255)
                                                                                                                                tctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGACTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTCTTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGTTCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTCCAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                                         <  2:189603/250‑1 (MQ=255)
                                                                                                                                                                  ggAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatgat                                                                      <  2:59248/251‑1 (MQ=255)
                                                                                                                                                                   gAATACCGTGCCCCGTCAGTGGTCGCGAGCAGCATCTTTACTCGATTTACTGCAAAATGGTGCTCCAAGCGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGATCGGCTATCAGTCTTTGCACACCTCGATGATc                                                                     >  2:256569/1‑251 (MQ=255)
                                                                                                                                                                                      tGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCgg                                                  >  1:278306/1‑251 (MQ=255)
                                                                                                                                                                                        gTCGCGAGAAGAATCTTTATTCGATTTACTGCAAAATGGTCCTCAAAGAGCAGCGTCTTTACTCGCTCATGGACATCTACGCTTTCCGAGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCAAGATGCACAGCCTGTACAAGCCGCGTCCGGGACGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGtt                                                >  1:366390/1‑251 (MQ=255)
                                                                                                                                                                                            cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACACTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAgg                                            >  2:182479/1‑251 (MQ=255)
                                                                                                                                                                                               gaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTcc                                         >  2:107775/1‑251 (MQ=255)
                                                                                                                                                                                               gaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGAAGCGTTTTTACTCGATCATGGACATATACGCTTTCAGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTACAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGAGTGCCGGTTGAGGTcc                                         >  1:334959/1‑251 (MQ=255)
                                                                                                                                                                                                                    cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgg                    >  1:337418/1‑251 (MQ=255)
                                                                                                                                                                                                                            tGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCAAGATCCGTACCGAAGATATGGACCAGATg            >  1:358814/1‑251 (MQ=255)
                                                                                                                                                                                                                             ggTGCTCACAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCCTTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATgg           >  2:64426/1‑251 (MQ=255)
                                                                                                                                                                                                                                      agagCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATgg  <  2:292714/251‑1 (MQ=255)
                                                                                                                                                                                                                                      agagCAGAGTTTTTGCTCGATCTAGGAAATCTACGCTTTCAGCGAGATCGTCAATGATTCTGACACCTGTTATCGCGTGCAGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATAACGCCATTCCAAAAGCGAACGGCTATAAGCCTTTGCACACATCGATGATCGGCCCGCACGGTGTGCAGGGTGAGGTACAGATCCGTACCGAAGATATGGATCAGAtagagtcgcagg  >  2:230537/1‑240 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
GGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGG  >  NC_000913/3822919‑3823399

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