Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 487,837 G→C 5.4% T155T (ACC→ACG)  pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947487,8370GC 5.4% 296.9 / 6.2 92T155T (ACC→ACG) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base C (3/2);  total (48/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 = 2.36e-01

GATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGCAGCCTT  >  NC_002947/487693‑487985
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                             <  5:113455/149‑1 (MQ=255)
                     aaaCTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTtgtg                                                                                                                             <  2:398903/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                                                                                                                    >  6:31700/1‑149 (MQ=255)
                              gcgcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCggg                                                                                                                    >  2:85893/1‑149 (MQ=255)
                                gcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAg                                                                                                                                              <  6:112690/121‑1 (MQ=255)
                                gcATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAg                                                                                                                                              >  5:112690/1‑121 (MQ=255)
                                     tCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTc                                                                                                             <  4:248642/149‑1 (MQ=255)
                                     tCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTc                                                                                                             <  6:143542/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                                                                       >  1:381791/1‑149 (MQ=255)
                                                                           tGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGc                                                                       >  5:102668/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)
                                                                                                  cACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATG‑CCGGAAAAGGCGATACCGCTTTCATTGATCACACCCt                                                     <  3:233575/143‑1 (MQ=255)
                                                                                                  cACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATG‑CCGGAAAAGGCGATACCGCTTTCATTGATCACACCCt                                                     >  4:233575/1‑143 (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                                                                 >  5:357531/1‑125 (MQ=255)
                                                                                                         tGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGt                                                                 <  6:357531/125‑1 (MQ=255)
                                                                                                         tGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCgggg                                         >  4:21810/1‑149 (MQ=255)
                                                                                                            cAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCg                                      >  6:343018/1‑149 (MQ=255)
                                                                                                             aGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGt                                     <  8:26121/149‑1 (MQ=255)
                                                                                                             aGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGt                                     >  8:234338/1‑149 (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)
                                                                                                                                              cgGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGCAGcc    <  5:181362/149‑1 (MQ=255)
                                                                                                                                                ggTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGCAGCCtt  <  1:64445/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCACCAGGGCCACGCGCAGGCCGCGGGTGGCCAGCATCATTACCACTTGTGCGGTGCGGGTCAGTTCGGCAAGGAATTCGGTATGCCCGGAAAAGGCGATACCGCTTTCGTTGATCACACCCTTGTGCACCGGGGCCGTGATCATCCCGGCGAAGTGCCCGTCCAGGCAGCCTT  >  NC_002947/487693‑487985

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: ‑