Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,566,493 G→C 5.8% A131G (GCC→GGC)  eutA ← reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,566,4930GC 5.8% 301.3 / 9.0 86A131G (GCC→GGC) eutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base G (40/41);  new base C (2/3);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

CGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTC  >  NC_000913/2566349‑2566641
                                                                                                                                                |                                                                                                                                                    
cGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGAt                                                                                                                                                  <  2:409788/149‑1 (MQ=255)
cGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGAt                                                                                                                                                  <  4:423512/149‑1 (MQ=255)
    caccGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGATCGCCGGTCACg                                                                                                                                              <  1:716119/149‑1 (MQ=255)
      ccGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGa                                                                                                                                            >  1:298189/1‑149 (MQ=255)
       cGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATATCGATATTCAGTACCCGACACAGCCGTTGCCCAGGAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGACCACGGAt                                                                                                                                           >  2:443003/1‑149 (MQ=255)
           gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGATCGCCGGTCACGGagccgg                                                                                                                                       <  2:883885/149‑6 (MQ=255)
           gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGATCGCCGGTCACGGagccgg                                                                                                                                       >  1:883885/1‑144 (MQ=255)
           gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGa                                                                                                                                       <  4:1002110/149‑1 (MQ=255)
           gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGa                                                                                                                                       <  4:304932/149‑1 (MQ=255)
           gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGa                                                                                                                                       <  2:1056072/149‑1 (MQ=255)
           gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGa                                                                                                                                       <  4:52712/149‑1 (MQ=255)
           gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGAGCCCCGGCTCCGTGACCGGCGATCACGGATTCGa                                                                                                                                       <  4:518739/149‑1 (MQ=255)
            ttGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAg                                                                                                                                      <  1:430237/149‑1 (MQ=255)
                 gcaggcagTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCg                                                                                                                                 >  4:580893/1‑149 (MQ=255)
                        gTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGgc                                                                                                                          >  2:280901/1‑149 (MQ=255)
                            cGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTg                                                                                                                      >  2:115402/1‑149 (MQ=255)
                                  tttttCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                                >  3:285115/1‑149 (MQ=255)
                                   ttttCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGa                                                                                                               <  4:285115/149‑1 (MQ=255)
                                            aTCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGc                                                                                                      >  1:12156/1‑149 (MQ=255)
                                            aTCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGc                                                                                                      <  1:293189/149‑1 (MQ=255)
                                                 aCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGc                                                                                                 >  2:839179/1‑149 (MQ=255)
                                                 aCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGc                                                                                                 <  2:298189/149‑1 (MQ=255)
                                                 aCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGc                                                                                                 >  1:594451/1‑149 (MQ=255)
                                                  cAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCg                                                                                                >  4:669451/1‑149 (MQ=255)
                                                     ggCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACt                                                                                             >  4:813801/1‑149 (MQ=255)
                                                      gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATc                                                                                                        >  3:285901/1‑137 (MQ=255)
                                                      gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATc                                                                                                        <  4:285901/137‑1 (MQ=255)
                                                       cGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTga                                                                                           <  2:557842/149‑1 (MQ=255)
                                                         tAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgaga                                                                                         >  3:152439/1‑149 (MQ=255)
                                                           gTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgagaga                                                                                       >  3:851205/1‑149 (MQ=255)
                                                                   tGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATc                                                                               <  2:123139/149‑1 (MQ=255)
                                                                       aCCGCCGATGTCGATATTCAGTACCCGACACAGCTGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccg                                                                           >  4:888953/1‑149 (MQ=255)
                                                                       aCCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccg                                                                           >  4:339654/1‑149 (MQ=255)
                                                                        ccgccgATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGccc                                                                                                                             <  2:634925/98‑1 (MQ=255)
                                                                        ccgccgATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGccc                                                                                                                             >  1:634925/1‑98 (MQ=255)
                                                                          gccgATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccgccg                                                                        <  2:202201/149‑1 (MQ=255)
                                                                              aTGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTga                                                                                           >  3:656871/1‑126 (MQ=255)
                                                                              aTGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTga                                                                                           <  4:656871/126‑1 (MQ=255)
                                                                              aTGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCg                                                                    >  4:331523/1‑149 (MQ=255)
                                                                                gTCGATATTGAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAg                                                                  >  3:754589/1‑149 (MQ=255)
                                                                                 tCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGc                                                                 >  1:367429/1‑149 (MQ=255)
                                                                                   gATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCAt                                                               >  4:438840/1‑149 (MQ=255)
                                                                                    atatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCAtt                                                              >  2:217135/1‑149 (MQ=255)
                                                                                       ttCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCg                                                           >  4:656537/1‑149 (MQ=255)
                                                                                        tCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCgg                                                          <  4:128372/149‑1 (MQ=255)
                                                                                        tCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCgg                                                          <  3:813801/149‑1 (MQ=255)
                                                                                        tCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCgg                                                          <  1:280901/149‑1 (MQ=255)
                                                                                            tACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGttt                                                      <  2:591416/149‑1 (MQ=255)
                                                                                            tACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGttt                                                      >  2:710385/1‑149 (MQ=255)
                                                                                              cccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTc                                                    <  4:114896/149‑1 (MQ=255)
                                                                                              cccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTc                                                    >  3:187856/1‑149 (MQ=255)
                                                                                              cccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTc                                                    <  2:1102558/149‑1 (MQ=255)
                                                                                               ccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcg                                                   >  3:924060/1‑149 (MQ=255)
                                                                                                  acacAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCgcgc                                                >  4:893356/1‑149 (MQ=255)
                                                                                                    acaGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGcc                                                                                  >  3:382957/1‑113 (MQ=255)
                                                                                                    acaGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGcc                                                                                  <  4:382957/113‑1 (MQ=255)
                                                                                                       gCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGa                                                                                                               <  4:106219/81‑1 (MQ=255)
                                                                                                       gCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGa                                                                                                               >  3:106219/1‑81 (MQ=255)
                                                                                                         cGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCAcc                                         >  2:992602/1‑149 (MQ=255)
                                                                                                          gttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGc                                                                 >  3:40843/1‑124 (MQ=255)
                                                                                                          gttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGc                                                                 <  4:40843/124‑1 (MQ=255)
                                                                                                          gttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCg                                        <  4:324283/149‑1 (MQ=255)
                                                                                                            tcttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGt                                      <  2:319235/147‑1 (MQ=255)
                                                                                                               tCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgat                                   <  1:217135/149‑1 (MQ=255)
                                                                                                                cAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatg                                  <  2:22167/149‑1 (MQ=255)
                                                                                                                 aGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatga                                 <  2:950640/149‑1 (MQ=255)
                                                                                                                 aGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatga                                 >  1:736582/1‑149 (MQ=255)
                                                                                                                  gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGc                                                                 >  3:1010656/1‑116 (MQ=255)
                                                                                                                  gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGc                                                                 <  4:1010656/116‑1 (MQ=255)
                                                                                                                  gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgat                                >  2:1076452/1‑149 (MQ=255)
                                                                                                                     aGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTgatgatgat                             <  3:557896/149‑1 (MQ=255)
                                                                                                                       ggTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATgg                           <  1:1076452/149‑1 (MQ=255)
                                                                                                                          ttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCAc                        <  2:514568/149‑1 (MQ=255)
                                                                                                                          ttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCAc                        <  3:68993/149‑1 (MQ=255)
                                                                                                                                 ccGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATc                 >  1:215371/1‑149 (MQ=255)
                                                                                                                                    gCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAAc              <  4:187856/149‑1 (MQ=255)
                                                                                                                                    gCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAAc              <  1:189103/149‑1 (MQ=255)
                                                                                                                                    gCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAAc              <  1:710385/149‑1 (MQ=255)
                                                                                                                                       ccGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTc                                                    >  2:144227/1‑108 (MQ=255)
                                                                                                                                       ccGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTc                                                    <  1:144227/108‑1 (MQ=255)
                                                                                                                                         gTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGt                                      <  2:352765/120‑1 (MQ=255)
                                                                                                                                         gTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGt                                      >  1:352765/1‑120 (MQ=255)
                                                                                                                                         gTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCttt         >  3:6328/1‑149 (MQ=255)
                                                                                                                                         gTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCttt         >  1:19978/1‑149 (MQ=255)
                                                                                                                                          tGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCGCCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTc        >  3:838678/1‑149 (MQ=255)
                                                                                                                                           gACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTcc       <  1:839179/149‑1 (MQ=255)
                                                                                                                                               ggCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGCGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCt   >  2:884248/1‑149 (MQ=255)
                                                                                                                                               ggCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCAGCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCt   >  3:710165/1‑149 (MQ=255)
                                                                                                                                               ggCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCt   <  2:113153/149‑1 (MQ=255)
                                                                                                                                               ggCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCt   >  2:26899/1‑149 (MQ=255)
                                                                                                                                                gCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGACGATGATGGCACCAGAATCAACGCTTTCCGGCTc  <  4:590918/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTC  >  NC_000913/2566349‑2566641

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