Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,598,497 T→C 5.7% N89D (AAC→GAC) ‡ dctB ← C4‑dicarboxylate transport sensor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,598,4970TC 5.7% 256.8 / 4.3 87N89D (AAC→GAC) ‡dctBC4‑dicarboxylate transport sensor protein
Reads supporting (aligned to +/- strand):  ref base T (39/43);  new base C (3/2);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TAGGGGCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTGCTGGCTGGCACGCTTGGCGTC  >  NC_002947/1598349‑1598644
                                                                                                                                                    |                                                                                                                                                   
tAGGGGCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGtt                                                                                                                                                     >  4:316182/1‑149 (MQ=255)
  ggggCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTc                                                                                                                                                    <  1:68831/148‑1 (MQ=255)
  ggggCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTc                                                                                                                                                    <  7:351032/148‑1 (MQ=255)
  ggggCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCt                                                                                                                                                   <  4:381749/149‑1 (MQ=255)
   gggCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTg                                                                                                                                                  <  7:347221/149‑1 (MQ=255)
    ggCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGa                                                                                                                                                 <  1:51873/149‑1 (MQ=255)
      cGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAAc                                                                                                                                               <  4:344306/149‑1 (MQ=255)
       ggAAACCGTAGTTGGCGCCGACGTAGCTGCTGCGCTGCAGCCAGTTGCTGTCGGCTATGGACAGCCCGGTGCGGTGGAGCAGCTCGAGGGTGGATGAATTGGCGTCGCCATTGATGCGTTCGAGCTTGAGGTTCAGAACGCCCTGAAcc                                                                                                                                              <  5:404510/149‑1 (MQ=255)
          aaCCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGAACGCCCTGAACCtca                                                                                                                                           <  4:196509/149‑3 (MQ=255)
           aCCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCtttt                                                                                                                                          >  8:332687/1‑149 (MQ=255)
             cGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTtcg                                                                                                                                        >  8:84940/1‑149 (MQ=255)
                  ttGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGAc                                                                                                                                   >  2:165693/1‑149 (MQ=255)
                  ttGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGAc                                                                                                                                   >  4:147618/1‑149 (MQ=255)
                    ggAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCCTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCg                                                                                                                                 >  7:71508/1‑149 (MQ=255)
                          gACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCcgc                                                                                                                           <  1:200367/149‑1 (MQ=255)
                            cGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCgcgc                                                                                                                         <  1:263882/149‑1 (MQ=255)
                               aGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAg                                                                                                                      >  7:144377/1‑149 (MQ=255)
                                    ctgGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTg                                                                                                                  <  8:144377/148‑1 (MQ=255)
                                    ctgGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGc                                                                                                                 >  2:209578/1‑149 (MQ=255)
                                      gggCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGa                                                                                                               >  2:264509/1‑149 (MQ=255)
                                      gggCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGAGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTtc                                                                                                                                         >  4:333336/1‑123 (MQ=255)
                                      gggCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGAGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTtc                                                                                                                                         <  3:333336/123‑1 (MQ=255)
                                       ggCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGTTGCGAt                                                                                                              <  3:147618/149‑1 (MQ=255)
                                        gCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCTCAAGGCTGCGATc                                                                                                             <  1:165693/149‑1 (MQ=255)
                                        gCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATc                                                                                                             <  3:23106/149‑1 (MQ=255)
                                        gCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATc                                                                                                             <  7:332687/149‑1 (MQ=255)
                                        gCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATc                                                                                                             <  5:99180/149‑1 (MQ=255)
                                          aGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAg                                                                                                           <  2:291347/149‑1 (MQ=255)
                                            ccgccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCt                                                                                                         >  1:64519/1‑149 (MQ=255)
                                             cgccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTc                                                                                                        >  2:2415/1‑149 (MQ=255)
                                              gccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTcc                                                                                                       >  6:212557/1‑149 (MQ=255)
                                              gccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCCTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTcc                                                                                                       <  8:71508/149‑1 (MQ=255)
                                                cAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCgg                                                                                                     >  1:21472/1‑149 (MQ=255)
                                                     gCTGGCGGCAATGGCCAGCCCGGTGTGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTc                                                                                                  >  4:384209/1‑147 (MQ=255)
                                                       tGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAgc                                                                                              >  1:264673/1‑149 (MQ=255)
                                                         gcggcAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTg                                                                                                                  >  8:71925/1‑127 (MQ=255)
                                                         gcggcAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTg                                                                                                                  <  7:71925/127‑1 (MQ=255)
                                                         gcggcAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTc                                                                                                  <  6:341937/143‑1 (MQ=255)
                                                         gcggcAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTc                                                                                                  >  5:341937/1‑143 (MQ=255)
                                                         gcggcAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAgc                                                                                              <  8:17/147‑1 (MQ=255)
                                                         gcggcAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAgc                                                                                              <  2:51719/147‑1 (MQ=255)
                                                              aaTGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGtt                                                                                                                                                     >  3:67846/1‑87 (MQ=255)
                                                              aaTGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGtt                                                                                                                                                     <  4:67846/87‑1 (MQ=255)
                                                              aaTGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCggg                                                                                                    <  8:146204/136‑1 (MQ=255)
                                                              aaTGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCggg                                                                                                    >  7:146204/1‑136 (MQ=255)
                                                              aaTGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCAATGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCa                                                                                       >  1:111507/1‑149 (MQ=255)
                                                                  gCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCAcggc                                                                                   >  6:82655/1‑149 (MQ=255)
                                                                    cAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCAcggcgg                                                                                 <  5:212557/149‑1 (MQ=255)
                                                                    cAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCAcggcgg                                                                                 <  1:2415/149‑1 (MQ=255)
                                                                    cAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCAcggcgg                                                                                 >  1:287862/1‑149 (MQ=255)
                                                                    cAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCAcggcgg                                                                                 <  2:264673/149‑1 (MQ=255)
                                                                    cAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAACTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCAGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTc                                                                                                  >  4:157546/1‑132 (MQ=255)
                                                                    cAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAACTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCAGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTc                                                                                                  <  3:157546/132‑1 (MQ=255)
                                                                        ccGGTGGGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAg                                                                             <  5:207259/149‑1 (MQ=255)
                                                                         cGGTGCGGTCGAGCAGCTCGGGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAgg                                                                            <  1:17005/149‑1 (MQ=255)
                                                                               ggTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACgg                                                                      <  2:64519/149‑1 (MQ=255)
                                                                                  cGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAg                                                                   <  3:316182/149‑1 (MQ=255)
                                                                                   gAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGc                                                                  >  8:83217/1‑149 (MQ=255)
                                                                                   gAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGc                                                                  >  4:344535/1‑149 (MQ=255)
                                                                                    agcagcTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGa                                                                              >  4:78053/1‑136 (MQ=255)
                                                                                    agcagcTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGa                                                                              <  3:78053/136‑1 (MQ=255)
                                                                                     gcagcTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGcga                                                                <  8:53930/149‑1 (MQ=255)
                                                                                          tCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGcgatcgat                                                           >  1:235060/1‑149 (MQ=255)
                                                                                                  ggaggaATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGgtgt                                                   >  8:178357/1‑149 (MQ=255)
                                                                                                   gaggaATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGtgtg                                                  <  3:384209/149‑1 (MQ=255)
                                                                                                     ggaATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGtgtg                                                  >  5:247165/1‑147 (MQ=255)
                                                                                                      gaATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACAGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAg                                               >  5:221319/1‑149 (MQ=255)
                                                                                                        aTTGGCGGAGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCCCGGCGGGGAGGGCACGGTCGCGATCGATCAGTGTGTGCCGTg                                             >  8:159721/1‑149 (MQ=255)
                                                                                                            gcggcgCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGtt                                         >  8:190328/1‑149 (MQ=255)
                                                                                                                cgcCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCg                                     >  6:421693/1‑149 (MQ=255)
                                                                                                                cgcCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCg                                     <  7:190328/149‑1 (MQ=255)
                                                                                                                        aTGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGTTCAGGGTGTGCAGTGCGTTGGCGTACAGGcc                             >  7:156554/1‑149 (MQ=255)
                                                                                                                         tGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGcc                                                                                           <  6:118780/86‑1 (MQ=255)
                                                                                                                            gTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGt                         >  5:387718/1‑149 (MQ=255)
                                                                                                                             ttCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGcgatcgat                                                           <  2:38174/114‑1 (MQ=255)
                                                                                                                             ttCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGtt                        <  5:82655/149‑1 (MQ=255)
                                                                                                                               cGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTg                       <  2:21472/148‑1 (MQ=255)
                                                                                                                                gAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgc                      <  7:83217/148‑1 (MQ=255)
                                                                                                                                 acgaaaaGGTTCAGAACGCCCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgctg                    >  7:338285/7‑149 (MQ=255)
                                                                                                                                   cTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgctggc                  <  8:156554/149‑1 (MQ=255)
                                                                                                                                     tGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgctggctg                >  5:155280/1‑149 (MQ=255)
                                                                                                                                     aaaGGTTCAGAACGCCCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgctggctg                >  6:31165/3‑149 (MQ=255)
                                                                                                                                     aaaGGTTCAGAACGCCCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgctggctg                >  5:390975/3‑149 (MQ=255)
                                                                                                                                      gAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgctggct                 <  1:209578/147‑1 (MQ=255)
                                                                                                                                      gAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTgctggct                 <  1:264509/147‑1 (MQ=255)
                                                                                                                                      gAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTGctggctgg               >  4:339924/1‑149 (MQ=255)
                                                                                                                                         gTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCATTGGCGTACAGGCCCAGTTGCTGGCTGGCAc            >  3:45801/1‑149 (MQ=255)
                                                                                                                                           tCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTGCTGGCTGGCACGc          >  6:341655/1‑149 (MQ=255)
                                                                                                                                             aGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTGCTGGCTGGCACGCtt        <  3:344535/149‑1 (MQ=255)
                                                                                                                                                gagTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTGCTGGCTGGCACGCTTGGc     >  3:346707/3‑149 (MQ=255)
                                                                                                                                                   ttCTGAACCTTTTAGTCGACAGGGCCGCGCAAGGCTGCGATCAGCGCCGGCACCAGCGCCAGCACGGCGGGGAGGGCACGGTGGCGCTCGAACAGGGTGTGGAGTGCGTTGCCGGACAGGCCGAGGTGCTGGCTGGCACGCTTGGcggc  >  4:89240/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAGGGGCGGAAACCGTAGTTGGAGCCGACGTAGCTGCTGGGCAGCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGGGTCCAGCGCCAGCACGGCGGGGAGGGCACGGTAGCGATCGATCAGGGTGTGCAGTGCGTTGGCGTACAGGCCCAGTTGCTGGCTGGCACGCTTGGCGTC  >  NC_002947/1598349‑1598644

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