Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,176,625 G→A S70F (TCC→TTC)  cpdA ← cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,6250GA100.0% 233.5 / NA 70S70F (TCC→TTC) cpdAcAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/40);  total (30/40)

GGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACC  >  NC_000913/3176378‑3176867
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ggCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAAtt                                                                                                                                                                                                                                                 >  2:96304/1‑251 (MQ=255)
 cccAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCAGCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCATTCATCGAAATGCTGATAGGCCGCAGAGAATTg                                                                                                                                                                                                                                                <  2:450164/250‑1 (MQ=255)
   cAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCGATCCAGCAACAGGATTTGCCATTGCTCACCAATATACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGTTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGAt                                                                                                                                                                                                                                              <  1:289054/251‑1 (MQ=255)
     gTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTAACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTAACCAATAAACACGCGCTTCGCCGGGGAGATACACGCATCCTGTAACGAGCTGTACATCGCGGGCTGGAAATCGTGGTTGCACGGCAGACAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGACCGCAGAGAATTGATcc                                                                                                                                                                                                                                            >  2:149887/1‑251 (MQ=255)
      tttACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATcc                                                                                                                                                                                                                                            <  1:161411/250‑1 (MQ=255)
               cAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAt                                                                                                                                                                                                                                  <  2:130637/251‑1 (MQ=255)
               cAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAt                                                                                                                                                                                                                                  <  1:264543/251‑1 (MQ=255)
                       cAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGt                                                                                                                                                                                                                           <  2:334206/250‑1 (MQ=255)
                       cAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCAGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGTCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTc                                                                                                                                                                                                                          >  1:289834/1‑251 (MQ=255)
                        aaGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTcg                                                                                                                                                                                                                         <  2:178813/251‑1 (MQ=255)
                        aaGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATGTGCCATTGCTCACCAATAAACACGCGCGTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACGTGCGATGCCTTCAGCGAAAGGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAACTCACCTGTcg                                                                                                                                                                                                                         <  2:121922/251‑1 (MQ=255)
                         aGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTcgc                                                                                                                                                                                                                        >  2:361208/1‑251 (MQ=255)
                             gAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGa                                                                                                                                                                                                                      <  1:170633/249‑1 (MQ=255)
                             gAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCAGTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGTCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACa                                                                                                                                                                                                                    <  1:66543/251‑1 (MQ=255)
                             gAAACGCGCTCAGCTCACCGTGCGGCCCGCCAAACTCCTGGCTCTCCAGCAACATGATTTGCCATTGCTCCCCAATAAACACGCGCTTCGCCGGGGAGATACCCGCTTCCTGTAACGCGCTGTACAACGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAACACACCTGTCGCGAc                                                                                                                                                                                                                     <  2:329163/250‑1 (MQ=255)
                                aCTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAGTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATc                                                                                                                                                                                                                 <  1:332563/251‑1 (MQ=255)
                                       cagctccCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCCCCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAAGTGATCCTGCGCTAAACCACCTTTCTCGAGAATCAGGTCGa                                                                                                                                                                                                          <  2:57256/251‑1 (MQ=255)
                                       cagctcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTc                                                                                                                                                                                                            <  1:356263/249‑1 (MQ=255)
                                          ctcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAAtt                                                                                                                                                                                                       >  1:157797/1‑251 (MQ=255)
                                          ctcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAAtt                                                                                                                                                                                                       >  1:211041/1‑251 (MQ=255)
                                                gTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCt                                                                                                                                                                                                 >  2:208631/1‑251 (MQ=255)
                                                    ggCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGgtg                                                                                                                                                                                             <  2:157797/251‑1 (MQ=255)
                                                         gCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGACGGGGAGATAACCGCCTCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGCTAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACACTCAGGTCGAATTCGTGCTGGTGTGGcc                                                                                                                                                                                        >  2:2925/1‑251 (MQ=255)
                                                                 ccTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATc                                                                                                                                                                                   <  2:275825/248‑1 (MQ=255)
                                                                 ccTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATc                                                                                                                                                                                   >  1:275825/1‑248 (MQ=255)
                                                                      cTATCCAGCAACAGGCTTTGCCATTGCTCACCAATAACCACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGCCGCAGGGCGCACGAACACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTcc                                                                                                                                                                             >  2:266310/1‑249 (MQ=255)
                                                                      cTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTcc                                                                                                                                                                             >  1:263724/1‑249 (MQ=255)
                                                                                   ggATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAGCTTGCGATGCCTTCAGCGAAATGCTTAGAGGCCGCAGAGAATTGa                                                                                                                                                                                                                                               <  2:66449/170‑1 (MQ=255)
                                                                                   ggATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATAGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGa                                                                                                                                                                                                                                               >  1:66449/1‑170 (MQ=255)
                                                                                         gCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACtctc                                                                                                                                                         <  2:359613/250‑1 (MQ=255)
                                                                                         gCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTcc                                                                                                                                                        >  2:234560/1‑251 (MQ=255)
                                                                                           cATTGATCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAAAGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATCGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGAAACTCTcccc                                                                                                                                                      >  2:458269/1‑250 (MQ=255)
                                                                                            aTTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAg                                                                                                                                                     <  1:266310/251‑1 (MQ=255)
                                                                                              tGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAgg                                                                                                                                                    <  2:338382/250‑1 (MQ=255)
                                                                                                          aaCACGAGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGGGTTTACCCCTa                                                                                                                                        >  2:50595/1‑250 (MQ=255)
                                                                                                            cacGCGCTTCGAAGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTCACTCTCCCAGGTGTTTCCCCCTAAc                                                                                                                                      >  2:385329/1‑250 (MQ=255)
                                                                                                               gcgcTTCGCCGGTGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACCATCCGGTCGAATTCGTGCTGGTGTGGCCGaa                                                                                                                                                                                     >  2:326877/1‑200 (MQ=255)
                                                                                                               gcgcTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGTCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCGCCCAGGTTTTTACCCCTAACAggg                                                                                                                                  <  2:86076/251‑1 (MQ=255)
                                                                                                               gcgcTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGaa                                                                                                                                                                                     <  1:326877/200‑1 (MQ=255)
                                                                                                                         ggggAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGACAGACGCAGGGCGCACGAAAACTTGAGATGCCTTCAGCGCAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACCATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCGAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTACCAggcactcgtgctt                                                                                                                        >  2:96494/1‑240 (MQ=255)
                                                                                                                         ggggAGATACCCGCATCCTGTAACGCGCAGTACAGCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCGTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAACTCACCTTTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGAGTTTACCCCTAACAGGGCTTCGTGCtt                                                                                                                        <  2:106522/251‑1 (MQ=255)
                                                                                                                           ggAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCtttt                                                                                                                      <  1:85735/251‑1 (MQ=255)
                                                                                                                                   ccGCATCCGGTAACGCGCTGTCCAGCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGAGGCAGGGAGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGGCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCGCCAGCACCGCCTGGTAACTCGCCCAGGTGTTTACCCCTAACAGGGCGTCGTGCTGTTGTGCAAAc                                                                                                              <  2:321083/251‑1 (MQ=255)
                                                                                                                                         ccTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                                         >  2:361335/1‑250 (MQ=255)
                                                                                                                                             gaaCGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAgtgt                                                                                                    <  1:396509/250‑1 (MQ=255)
                                                                                                                                                      gTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAAttt                                                                                           <  1:463548/251‑1 (MQ=255)
                                                                                                                                                      gTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAAt                                                                                             <  1:271395/249‑1 (MQ=255)
                                                                                                                                                        aCATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTg                                                                                          <  1:356216/250‑1 (MQ=255)
                                                                                                                                                        aCATCGCGGGATCGAAATCGGGATCGCCCGGCAGCCAGACGCAGAGCCCACGAAAACTTGCGATGCATTCAGCGAAATGCTGATAGGCCGCCGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGGGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTCAATCTCCCAGGTGTTTACCCCTAACAGGGCTACGTGCTTTTCAGCAAACAGATGAGAGTCGGTAATTTGt                                                                                         >  2:51800/1‑251 (MQ=255)
                                                                                                                                                         cagcgcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGGTGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAATGCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTa                                                                                        <  1:10961/248‑1 (MQ=255)
                                                                                                                                                                ggCTGGAAATCGAGGGTGCCCGGCAGCCAGACGCAGGTCGCACGAAAACTTGCGATGCCTTCAGCTAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTTTAAAATCCt                                                                                 <  2:162868/251‑1 (MQ=255)
                                                                                                                                                                    ggAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGTTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGGAAGTTGTAAAATCCTGa                                                                               <  1:146915/249‑1 (MQ=255)
                                                                                                                                                                        aTCGTGGTTGCCCGGCAGCCAGATGCAGGGCGCACGAAAACTTGCGATTCCTTCAGCGAAATTCTGATAGGCCGCAGAGAATGGATCCTGCGCTAAATCTCCGGTCGCGACAATCAGGTCGTATCCGTGCTGGTTTGGCCGAATCGCTTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCGTTTGTGCAAACAGGTGAGTGTCGTTAATTTGTAAAACCCTGTCTCTgg                                                                         <  2:394486/251‑1 (MQ=255)
                                                                                                                                                                                 gCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAcca                                                                 <  1:44077/250‑1 (MQ=255)
                                                                                                                                                                                                      cgcACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCAGGAGGAAGGGGTAACAGGc                                           >  2:211951/1‑251 (MQ=255)
                                                                                                                                                                                                             aaaCTTGCGATGAATTCAGCGAAATGCTGATAGGCCGCAGCGAATTGATCCTGAGCTAACTCACCTGTCGCGACAATCAGGTCGAATTAGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTACCAGGGGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGAATCTGGCCTCAACAGCCAGAGGAAGGGTTAACAAGCTTTCCa                                     >  2:327281/1‑250 (MQ=255)
                                                                                                                                                                                                                   gCGATGCCTTCAGCGAAATTCTGATAGGCCGCAGAGAATTGATCCTTCGCTAAATCACCGGGCGCGACAAGCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACGCTGGCCTCACCAGCCAGAGTAAGGGTTAACAGGCTTTCCAAATGgtg                              <  2:225007/251‑1 (MQ=255)
                                                                                                                                                                                                                   gCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgtg                              <  1:175484/251‑1 (MQ=255)
                                                                                                                                                                                                                       tGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTACCAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCt                          >  1:323014/1‑251 (MQ=255)
                                                                                                                                                                                                                            tCAGCGAAATGGTGATAGGCACCAGAGAATTGCCCCTGCGCTAAATCCCCTGTCGCGACACTCCGGTCGAATTCGTGCTGGCGTGGCCGAACCGCCTCCCGCACAGCCTGGCAACTATCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGCGCAAAAAGGCGAGTGTCGGAAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTTCAAATGGTGTCCTTAGGt                     >  2:451417/1‑251 (MQ=255)
                                                                                                                                                                                                                             cAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGtt                    >  1:375402/1‑251 (MQ=255)
                                                                                                                                                                                                                             cAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGt                     <  1:282430/250‑1 (MQ=255)
                                                                                                                                                                                                                               gCGAAATGCTGCTAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGCATTCGTGCTGGTGTGGCCGAATCGCCTCAAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTc                  >  1:310130/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGAAATGCTGATAGGCCGCAGAGAATTGATCAAGCGCTATATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGAGGCAGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACCCTGGCCTCACCAGCCAGAGGAAGGGTTATCAGGCTTTCCAAATGGTGTCCTTAGGTTTc                  >  2:242043/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaa                                                                                      <  1:41568/178‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGACGACTCGCCTCCAGCACCGCCTGGTCACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTCATTTGTaaa                                                                                      >  2:41568/1‑178 (MQ=255)
                                                                                                                                                                                                                                    aTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAGCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacgac             <  1:257966/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaa                                                                                      >  2:464469/1‑174 (MQ=255)
                                                                                                                                                                                                                                        tGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaa                                                                                      <  1:464469/174‑1 (MQ=255)
                                                                                                                                                                                                                                               ccGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATAAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGCGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAAcc  >  1:176897/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACC  >  NC_000913/3176378‑3176867

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