Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,647,677 C→T intergenic (+26/+28) gor → / ← dinQ glutathione reductase (NADPH)/membrane toxin DinQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,6770CT100.0% 219.0 / NA 69intergenic (+26/+28)gor/dinQglutathione reductase (NADPH)/membrane toxin DinQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/34);  total (35/34)

ATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACC  >  NC_000913/3647430‑3647911
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
gcgttgaCCGCCTTCCCCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCAGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGttct                                                                                                                                                                                                                                         <  2:10009/245‑1 (MQ=255)
    ataCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAg                                                                                                                                                                                                                                      <  1:38481/250‑1 (MQ=255)
         gccgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACAGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTa                                                                                                                                                                                                                                <  1:185708/251‑1 (MQ=255)
            gTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAAtt                                                                                                                                                                                                                             >  1:119371/1‑251 (MQ=255)
               cccacTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACg                                                                                                                                                                                                                          <  1:140607/250‑1 (MQ=255)
                 cacTCCCCGCCAGCCGTGCCGCATGAAGCTTGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATGGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGtt                                                                                                                                                                                                                        <  2:293493/251‑1 (MQ=255)
                    tCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAACGACTTCGACAATACCGTCGCCATTCACCCAACGGCGCCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCAAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTcc                                                                                                                                                                                                                     >  1:404512/1‑251 (MQ=255)
                     cACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATGCCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTcc                                                                                                                                                                                                                     <  1:244255/250‑1 (MQ=255)
                         gCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCAGTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTAt                                                                                                                                                                                                                <  1:207135/251‑1 (MQ=255)
                         gCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGGCGGTAGTCACGGCATTGGCGTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCATCGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTAt                                                                                                                                                                                                                <  1:526176/251‑1 (MQ=255)
                              ccGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGtt                                                                                                                                                                                                           >  1:426241/1‑251 (MQ=255)
                                  gCCGCATGAAACTGGTGTGCGTTGGAACGGAAGAGAAGATTGTCGGTATTAACGGCATTGGCTTTGGGATGGACGAAATGTTGCAGGGCTTCGCGGGGGCGCTGAAGATAGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATGCACCCAACGGCGGCAGAAGATGTCGTGAACATGAGTGAAATGTTAAAGGGCTAAGAGGAGTGTGTTCTTAAACCTTAATTACGTGTCCGCTATCAGTTc                                                                                                                                                                                                          >  2:417648/1‑248 (MQ=255)
                                       aTGAAGCTGGGGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTAGGGACGAAATTTTGCAGGGCTTCGCGGTGGCGCTTAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCTGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTg                                                                                                                                                                                                  <  2:558581/251‑1 (MQ=255)
                                             cTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTAc                                                                                                                                                                                                                           <  2:93572/220‑1 (MQ=255)
                                             cTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTAc                                                                                                                                                                                                                           >  1:93572/1‑220 (MQ=255)
                                                         ggagcGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCtt                                                                                                                                                                                                                                        <  2:564328/191‑1 (MQ=255)
                                                         ggATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCtt                                                                                                                                                                                                                                        >  1:564328/1‑195 (MQ=255)
                                                            tCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTc                                                                                                                                                                             >  1:62798/1‑251 (MQ=255)
                                                             cGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTcc                                                                                                                                                                            >  1:447839/1‑251 (MQ=255)
                                                               gaagagaagaTTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGc                                                                                                                                                                                       <  1:94234/238‑1 (MQ=255)
                                                               gaagagaagaTTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGGAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGc                                                                                                                                                                                       >  2:94234/1‑238 (MQ=255)
                                                                aagagaagaTTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGAAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGTGCTAGGAGTCGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAACAGAAAGCGGATCAGTTCcagc                                                                                                                                                                         >  2:40474/1‑251 (MQ=255)
                                                                 agcgaagaTTGACGGTAATCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGCAGATGGTGGCAACCACAAAAGACTTCGACAATCCCGTCGCCATTCACCAAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAGATGTTAAAAGGCTAAGCGAAGTGTGTTCTTAGACCTTACTTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAGGCGGATCAGTTCcacc                                                                                                                                                                         >  2:257644/4‑248 (MQ=255)
                                                                   agaagaTTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCagcagc                                                                                                                                                                      >  1:195856/1‑251 (MQ=255)
                                                                   agaagaTTGTAGGTATTCACGGCATTGGCTTTGGTATGGACGACATGATGCAGGGCTTAGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTAACCCAACGGCGGCAGAAAAGATCGTGACAATGCGTTAACTGTTAAAGGGCTAAGAGTAGTGTGTTCATAGCCCTTACTTACGTTTACGTTATAAGTTCAGCAGCTGACGCAGacatgggataaggtaccgc                                                                                                                                                                         >  2:349395/1‑230 (MQ=255)
                                                                       gaTTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCaa                                                                                                                                                                  <  2:119371/251‑1 (MQ=255)
                                                                          tGTCGGTCTTCACGGTATTGGCTTTGGTATGGACGAAATGGTGCAGGGCCTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCCCCCAAAGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAAtt                                                                                                                                                                <  1:417648/250‑1 (MQ=255)
                                                                               gTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGcc                                                                                                                                                          >  1:548760/1‑251 (MQ=255)
                                                                                    cACGGAATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGGAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCt                                                                                                                                                       >  2:123153/1‑249 (MQ=255)
                                                                                      cGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGa                                                                                                                                                    >  1:261739/1‑250 (MQ=255)
                                                                                       ggCATAGGCTTTGCTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAACGACTTCGACAATACCGTCGCCATTCACCCACCGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTCTCAGTTCCGAATCTGAAGCAGAAAGCGGCTCAGTTCCAACAGCGCCATTCACGCCCCTAGAAc                                                                                                                                                  >  2:271811/1‑251 (MQ=255)
                                                                                         cATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACga                                                                                                                                                <  1:499255/251‑1 (MQ=255)
                                                                                         cATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACAGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCACCTAGaa                                                                                                                                                   >  1:107241/1‑248 (MQ=255)
                                                                                            tGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACgatg                                                                                                                                              >  1:474787/1‑250 (MQ=255)
                                                                                                 ttGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGc                                                                                                                                         <  1:429018/250‑1 (MQ=255)
                                                                                                   ggTAAGGACGAAATGTTGCCGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTGCGTGACAATGCTTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGACCCTAGAACGATGATTGCttt                                                                                                                                      <  2:312167/251‑1 (MQ=255)
                                                                                                         gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaat                                                                                                                                <  1:429873/251‑1 (MQ=255)
                                                                                                             aaaTGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagca                                                                                                                                                                        <  1:210761/207‑1 (MQ=255)
                                                                                                             aaaTGTTGCAGGGCTTCGCGGGGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGAAAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGAGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGACGCTGAAGCAGAAAGCGGATCAGTTCcagca                                                                                                                                                                        >  2:210761/1‑207 (MQ=255)
                                                                                                                  ttGCAGAGCTTCGCGGTGGCGCTGAAGATGGGGGCCACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGAAAATACGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCACTTAATTACGTTTCCGATATCAGTTCATAAGCTGTAGCAGAAAGCGGATCAATTCCAGCAGCGCAATTCACGCCCCTAGCACGATGCTTGCTTTatacatccccagttt                                                                                                                       >  2:459803/1‑238 (MQ=255)
                                                                                                                   tGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagc                                                                                                                                                                         <  1:504621/200‑1 (MQ=255)
                                                                                                                         gCTTCGCGGTGGCGCTTAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCATCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCAt                                                                                                                <  1:42696/251‑1 (MQ=255)
                                                                                                                         gCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCAt                                                                                                                <  1:449337/251‑1 (MQ=255)
                                                                                                                             cgcgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCg                                                                                                             >  1:550110/1‑250 (MQ=255)
                                                                                                                                         gAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTgaga                                                                                                  <  1:419472/249‑1 (MQ=255)
                                                                                                                                              tGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATc                                                                                           >  1:497815/1‑251 (MQ=255)
                                                                                                                                              tGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGAAGAAGAGATCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGAAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGACCATAGAACAATGATTGCTTTATTAATCACCCGTTTTCTCCATGCGATGGAGCGAGCAGGCATc                                                                                           >  2:131978/1‑251 (MQ=255)
                                                                                                                                                       ccAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTc                                                                                  <  1:84353/251‑1 (MQ=255)
                                                                                                                                                          aaaaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATc                                                                               <  1:101157/251‑1 (MQ=255)
                                                                                                                                                                      aCAATTCCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCa                                                                   <  1:275689/251‑1 (MQ=255)
                                                                                                                                                                      aCAATACAGTCGCCATTCACCCAACGGCGGCAGCAGAGTTCGTGACAATGCGTTAACTGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCCTGCGATGGAGGGAGAATGCATCCGCTTCCTCATCCCCTGCCTGTcc                                                                   >  2:31039/1‑250 (MQ=255)
                                                                                                                                                                        aaTACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACg                                                                 >  1:203803/1‑251 (MQ=255)
                                                                                                                                                                           aCCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGgc                                                               <  2:195856/250‑1 (MQ=255)
                                                                                                                                                                                 gCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGc                                                                                              >  1:539799/1‑213 (MQ=255)
                                                                                                                                                                                 gCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGtc                                                        >  1:523929/1‑251 (MQ=255)
                                                                                                                                                                                            aaCGGCGGCAGAAGAGTTCGCGACAATGCGTTAAATGTTAAAGGGCTAAGATTAGTGTGTTCTTAGCCCTTAAGTTCGTTTCCTCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCACTAGAACGATGATTGCTTTACCAATCTCCCGTTTTCCCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCTCTGCCTGTCACGGCGCATGTCTCATTGTTAGa                                             <  2:472535/251‑1 (MQ=255)
                                                                                                                                                                                                gcggcAGAAGAGTGCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGGAGTGTGTTCTTAGCCCTTAATGTCGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAATCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTCTATCAATCACCCGGTTTCTCCATGCGATGGAGTGAGAATGCTTCCGCCTACTCATCCACTGCCTGTCTCGGCGCATGTCTCATTGTTAGATAAg                                         <  1:257644/251‑1 (MQ=255)
                                                                                                                                                                                                    cAGAAGAGTTCGTGACAAAGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAAct                                     >  1:150725/1‑251 (MQ=255)
                                                                                                                                                                                                          agTTCGTGACAGTGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCa                               <  1:479590/251‑1 (MQ=255)
                                                                                                                                                                                                          agTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCa                               <  2:497815/251‑1 (MQ=255)
                                                                                                                                                                                                              cGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTATCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCACTTAACGCCCCTAAAACGCTGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGCGTGTGAATGCTTCCGCTTACTCATCCACTGCCTGTCACGGCGAATGTCTCATTGTTAGATAAGAACTCTCTCACTcc                           >  2:39983/1‑251 (MQ=255)
                                                                                                                                                                                                                  aCAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGcc                       >  1:523113/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaTGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCag                     <  1:208284/251‑1 (MQ=255)
                                                                                                                                                                                                                       gCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCagag                   <  1:213010/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTTAAATGTTAAAGGGCGAAGAGCAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTGCAGAAGCTGCAGCAGAAAGCGGATCAGTTCGAGCAGAGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCACCATGCGATGGAGAGAGAATGCATCCGCTAACTCATCCAATGCCTGTCACAGCGCATGTCTCATTGTTCGATAAGAACTCTCTCACCCCAGCCAGAGccc                >  2:185148/1‑249 (MQ=255)
                                                                                                                                                                                                                          ttaaattgtAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGACCGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATTTCTCATTGTTAGATAATAACTCTCTCACTCCGGCCAGAGCAt                <  2:536636/242‑1 (MQ=255)
                                                                                                                                                                                                                              aTGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGtt           <  1:40474/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTGTTAAAGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGtt           >  1:140155/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aGGGCTAAGAGTAGTGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGAATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGcc    >  1:291256/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ggCTAAGAGTAGAGTGTTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCGCTCCGGCCAGAGCATCAGTTGACGGcacc  >  1:563881/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
ATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACC  >  NC_000913/3647430‑3647911

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