Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 487,833 C→T 5.6% G157S (GGC→AGC) ‡ pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947487,8330CT 5.6% 274.4 / 3.4 89G157S (GGC→AGC) ‡pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (45/39);  new base T (3/2);  total (48/41)
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.33e-01

CACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGCA  >  NC_002947/487687‑487980
                                                                                                                                                  |                                                                                                                                                   
cacGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCcgc                                                                                                                                                   <  7:373230/149‑1 (MQ=255)
cacGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCcgc                                                                                                                                                   <  6:166809/149‑1 (MQ=255)
cacGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCcgc                                                                                                                                                   <  8:292992/149‑1 (MQ=255)
 acGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCcgcg                                                                                                                                                  <  7:240805/149‑1 (MQ=255)
 acGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCcgcg                                                                                                                                                  >  8:286836/1‑149 (MQ=255)
      gATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACCCGCGGCCTGCGCGTgg                                                                                                                                             <  6:224173/149‑1 (MQ=255)
          cgcgGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCa                                                                                                                                          <  3:253240/148‑1 (MQ=255)
           gcgGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCATGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAg                                                                                                                                         >  2:307575/1‑148 (MQ=255)
           gcgGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAg                                                                                                                                         >  7:246905/1‑148 (MQ=255)
           gcgGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAg                                                                                                                                         >  3:192336/1‑148 (MQ=255)
                   gCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTAc                                                                                                                                >  7:97428/1‑149 (MQ=255)
                    cAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTAc                                                                                                                                >  6:59057/1‑148 (MQ=255)
                     aaTGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCAAGCGCAGGCCCCGGGTGGCCAGCATCATGACCa                                                                                                                              >  8:309338/1‑149 (MQ=255)
                     aaTGCCAAACTTGTCCCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGAAAGCGTCGGCAACATCACGCAGTCGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCa                                                                                                                              <  3:21532/149‑1 (MQ=255)
                     aaTGACAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGGAGGTGCGTGGTAACCAGGGCCACGCCCAGGCCGCGGGTGGCCAGCATAATTTcc                                                                                                                               >  4:21532/1‑148 (MQ=255)
                         caaaaCTTGTCGCGCATGTATGCATGCAGGATACGGGTTACACGTTCGACGATCTCTTCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACtt                                                                                                                           <  4:192336/146‑1 (MQ=255)
                           aaaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTtgtg                                                                                                                        >  1:227958/1‑149 (MQ=255)
                           aaaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTtgtg                                                                                                                        <  2:398903/149‑1 (MQ=255)
                           aaaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTtgtg                                                                                                                        <  5:113455/149‑1 (MQ=255)
                           aaaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTtgta                                                                                                                        <  5:113458/149‑2 (MQ=255)
                               ttGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGc                                                                                                                                        <  7:140959/129‑1 (MQ=255)
                               ttGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGc                                                                                                                                        >  8:140959/1‑129 (MQ=255)
                                  tCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTgtgcggtgc                                                                                                                  <  7:286836/148‑1 (MQ=255)
                                    gcgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCggg                                                                                                               >  6:25786/1‑149 (MQ=255)
                                    gcgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCggg                                                                                                               >  2:85893/1‑149 (MQ=255)
                                    gcgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCggg                                                                                                               >  6:31700/1‑149 (MQ=255)
                                      gcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAg                                                                                                                                         >  5:112690/1‑121 (MQ=255)
                                      gcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAg                                                                                                                                         <  6:112690/121‑1 (MQ=255)
                                           tCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTc                                                                                                        <  6:143542/149‑1 (MQ=255)
                                           tCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTc                                                                                                        <  4:248642/149‑1 (MQ=255)
                                               cATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTtgtg                                                                                                                        >  2:10845/1‑129 (MQ=255)
                                               cATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCCGGTGGCCAGCATCATTACCACTtgtg                                                                                                                        <  1:10845/129‑1 (MQ=255)
                                                 tGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAg                                                                                                  <  8:106052/149‑1 (MQ=255)
                                                   cAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCACTGGCAGGTGCGTGGTCACCAGGGCCACGCCCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAgg                                                                                                 >  4:81119/1‑148 (MQ=255)
                                                    aGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGaa                                                                                               <  1:85893/149‑1 (MQ=255)
                                                       aTACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCCTCATTTCCACTTGTGCGGGGCGGGGCAGTTCGGCAAGGAAGTc                                                                                            >  2:280438/1‑149 (MQ=255)
                                                           gggTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCATCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTa                                                                                        >  2:86213/1‑149 (MQ=255)
                                                           gggTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGt                                                                                                              >  2:25744/1‑127 (MQ=255)
                                                           gggTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGt                                                                                                              <  1:25744/127‑1 (MQ=255)
                                                           gggTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTa                                                                                        <  2:227958/149‑1 (MQ=255)
                                                           gggTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTa                                                                                        >  8:321528/1‑149 (MQ=255)
                                                                    gTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCa                                                                                                                              <  1:344799/102‑1 (MQ=255)
                                                                    gTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCa                                                                                                                              >  2:344799/1‑102 (MQ=255)
                                                                       cGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGaaaa                                                                             <  8:97428/148‑1 (MQ=255)
                                                                        gACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAgg                                                                           <  5:59057/149‑1 (MQ=255)
                                                                        gACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAg                                                                            <  2:53072/148‑1 (MQ=255)
                                                                             gcTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCAAAGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTATGCAGTGCGGGTCAGTTCGGCAAGGAAATCGGAAAGCCCGGAAAAGGCGATa                                                                      >  3:121013/1‑149 (MQ=255)
                                                                                 tGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGc                                                                  >  5:102668/1‑149 (MQ=255)
                                                                                 tGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGc                                                                  >  1:381791/1‑149 (MQ=255)
                                                                                   cAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCt                                                                 <  5:25786/148‑1 (MQ=255)
                                                                                   cAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCt                                                                 <  5:31700/148‑1 (MQ=255)
                                                                                     gTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTc                                                              <  8:246905/149‑1 (MQ=255)
                                                                                      tGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCg                                                             >  6:181362/1‑149 (MQ=255)
                                                                                        aTAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGtt                                                           <  3:316141/149‑1 (MQ=255)
                                                                                         tAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGtt                                                           <  2:381791/148‑1 (MQ=255)
                                                                                              tCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCTGGAAAAGGCGATACCGCTTTCGTTGATca                                                      >  7:68568/1‑148 (MQ=255)
                                                                                                ggCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCacac                                                   <  7:321528/149‑1 (MQ=255)
                                                                                                 gCAACATCATGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCTTTGATCACAcc                                                  <  1:307575/149‑1 (MQ=255)
                                                                                                 gCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTACGTTGATCACAcc                                                  >  5:241810/1‑149 (MQ=255)
                                                                                                     cATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTtg                                              >  5:170666/1‑149 (MQ=255)
                                                                                                        cACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCGGAAAAGGCGATACCGCTTTCATTGATCACACCCt                                                >  4:233575/1‑143 (MQ=255)
                                                                                                        cACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCGGAAAAGGCGATACCGCTTTCATTGATCACACCCt                                                <  3:233575/143‑1 (MQ=255)
                                                                                                         aCGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTACGTTGATCACACCCTTGTGc                                           <  6:241810/148‑1 (MQ=255)
                                                                                                          cGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCAc                                         >  7:26121/1‑149 (MQ=255)
                                                                                                           gCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTtgt                                             <  5:217872/144‑1 (MQ=255)
                                                                                                           gCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTtgt                                             >  6:217872/1‑144 (MQ=255)
                                                                                                           gCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCAcc                                        >  7:100753/1‑149 (MQ=255)
                                                                                                               tGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGt                                                            <  6:357531/125‑1 (MQ=255)
                                                                                                               tGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGt                                                            >  5:357531/1‑125 (MQ=255)
                                                                                                               tGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCgggg                                    >  4:21810/1‑149 (MQ=255)
                                                                                                                  cAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCg                                 >  6:343018/1‑149 (MQ=255)
                                                                                                                   aGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGt                                >  8:234338/1‑149 (MQ=255)
                                                                                                                   aGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGt                                <  8:26121/149‑1 (MQ=255)
                                                                                                                   aGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGt                                <  7:309338/149‑1 (MQ=255)
                                                                                                                     gTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGa                              <  1:86213/149‑1 (MQ=255)
                                                                                                                         gTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGatcat                          >  8:374224/1‑149 (MQ=255)
                                                                                                                          tGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGatcatc                         >  2:398009/1‑149 (MQ=255)
                                                                                                                           ggTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATcc                        >  5:158787/1‑149 (MQ=255)
                                                                                                                            gTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATccc                       >  7:267729/1‑149 (MQ=255)
                                                                                                                              cACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCgg                     <  5:343018/149‑1 (MQ=255)
                                                                                                                               aCCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGccc                                                                                   <  2:124768/86‑1 (MQ=255)
                                                                                                                               aCCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGc                    >  2:64445/1‑149 (MQ=255)
                                                                                                                                 cAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCg                   >  4:393226/1‑148 (MQ=255)
                                                                                                                                   tggCCACCCGCGGCCTGCGCGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAg                >  8:131423/2‑149 (MQ=255)
                                                                                                                                   tggCCACCCGCGGCCTGCGCGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAg                >  7:295206/2‑149 (MQ=255)
                                                                                                                                     gCCACCCGCGGCCTGCGCGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCg                                 <  2:39141/130‑1 (MQ=255)
                                                                                                                                     gCCACCCGCGGCCTGCGCGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCg                                 >  1:39141/1‑130 (MQ=255)
                                                                                                                                        aCGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCTGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGccc           <  8:68568/149‑1 (MQ=255)
                                                                                                                                        aCGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGccc           >  6:215655/1‑149 (MQ=255)
                                                                                                                                        aCGCGCAGGCCGCGGGTGGCCAGCATCATTACCAATTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGccc           >  2:184037/1‑149 (MQ=255)
                                                                                                                                         cgcgCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCg          <  3:250367/149‑1 (MQ=255)
                                                                                                                                           cgcAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTc        >  2:65804/1‑149 (MQ=255)
                                                                                                                                              attccGCGGGTGGACAGCATCATTACCACTTGTGCGGTTCGGGTCAGTTCGGCAAGTAATTCGGTATGCCCGGAAAAGTCGATACCGCTTTCGTTTATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAg     <  8:267729/146‑1 (MQ=255)
                                                                                                                                              aGGCCGCGGGTGGCCAGCATCATTACCACTTGTGTGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAg     >  6:316798/1‑149 (MQ=255)
                                                                                                                                              aGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAg     <  3:21810/149‑1 (MQ=255)
                                                                                                                                                 ccGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGTGGCTGTTCGGCTAGGATTTGTGTATGCCGGGTAAAGGCGATACCGGTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTTATCATCCGTGCGAAGTGCCCGTCCAGGCa  <  4:121013/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGCA  >  NC_002947/487687‑487980

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: