Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,940,561 T→C 5.5% S193P (TCG→CCG) ‡ fau → 5‑formyltetrahydrofolate cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,940,5610TC 5.5% 869.2 / 23.3 272S193P (TCG→CCG) ‡fau5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base T (136/117);  new base C (4/11);  total (140/132)
Fisher's exact test for biased strand distribution p-value = 6.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 283 total aligned reads displayed.

TGGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCA  >  NC_002947/5940412‑5940710
                                                                                                                                                     |                                                                                                                                                     
 gggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                      <  1:330139/150‑1 (MQ=255)
   gCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCgg                                                                                                                                                    >  1:355456/1‑150 (MQ=255)
    cGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGa                                                                                                                                                   <  2:2518149/150‑1 (MQ=255)
    cGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGa                                                                                                                                                   <  1:7329498/150‑1 (MQ=255)
    cGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGa                                                                                                                                                   <  1:1948649/150‑1 (MQ=255)
         gCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTg                                                                                                                                              <  1:2116725/150‑1 (MQ=255)
         gCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTg                                                                                                                                              <  2:5471412/150‑1 (MQ=255)
          cTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTgg                                                                                                                                             <  1:5749631/150‑1 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                            >  1:3704423/1‑150 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                            >  2:7380239/1‑150 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                            <  1:3593919/150‑3 (MQ=255)
            tCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCt                                                                                                                                           <  1:2779112/150‑1 (MQ=255)
             cTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTg                                                                                                                                          <  1:1752691/150‑1 (MQ=255)
              tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGTAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                         >  2:1114055/1‑150 (MQ=255)
              tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                         >  2:5980805/1‑150 (MQ=255)
              tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                         >  2:1709292/1‑150 (MQ=255)
                tGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTa                                                                                                                                       >  1:4784865/1‑150 (MQ=255)
                 gATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTg                                                                                                                                              >  1:2554755/1‑142 (MQ=255)
                 gATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTg                                                                                                                                              <  2:2554755/142‑1 (MQ=255)
                     gCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTgg                                                                                                                                  >  1:5504368/1‑150 (MQ=255)
                       aGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgcg                                                                                                                                <  1:5012980/150‑1 (MQ=255)
                           tGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCAc                                                                                                                            >  1:2478971/1‑150 (MQ=255)
                              ccTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGtt                                                                                                                         <  2:3492721/150‑1 (MQ=255)
                              ccTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGtt                                                                                                                         >  1:6187199/1‑150 (MQ=255)
                                tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                       >  2:6488162/1‑150 (MQ=255)
                                tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                       <  2:6973341/150‑1 (MQ=255)
                                tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                       >  2:531687/1‑150 (MQ=255)
                                tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                       <  2:1857967/150‑1 (MQ=255)
                                  ccAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgc                                                                                                                     >  2:6244852/1‑150 (MQ=255)
                                     ggCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCGGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgcggc                                                                                                                  <  1:6488661/150‑1 (MQ=255)
                                             ggCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgg                                                                                                          <  2:2819754/150‑1 (MQ=255)
                                              gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                         <  1:947478/150‑1 (MQ=255)
                                              gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                         >  2:6115660/1‑150 (MQ=255)
                                              gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                         >  1:615243/1‑150 (MQ=255)
                                               ctaaCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                        <  1:3943018/148‑1 (MQ=255)
                                                aaaCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggggg                                                                                                       >  2:5100946/1‑150 (MQ=255)
                                                aaaCCTGGAAAAAACCCTTGCTGCGGGGGCTGGCGCATGAAGGCCAGAAGGGCGAGCGGCTGGCGCAGGCCAGTGGGGAGGGTCCGCGGCAGGGCCCGGTCTCGGACCGTGGCGGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggggg                                                                                                       <  2:1259925/150‑1 (MQ=255)
                                                  aCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAAGCCCAGACGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGt                                                                                                     <  2:7050983/150‑1 (MQ=255)
                                                   ccTGGAAAAAACCCTTGCTGATGGGGCTGGCGAATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGATGGCATGTACCTCTGCAGGGCACGGTCTCTGACAGTCGCGGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTa                                                                                                    >  1:5617849/1‑150 (MQ=255)
                                                    cTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGtct                                                                                                   >  1:5724371/1‑148 (MQ=255)
                                                     tGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                  >  1:7167758/1‑150 (MQ=255)
                                                     tGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                  >  1:2230378/1‑150 (MQ=255)
                                                     tGGAAAAAACCCTTGCTGCTGGAGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCt                                                                                                                                    <  2:42740/116‑1 (MQ=255)
                                                        aaaaaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAAc                                                                                               >  2:834903/1‑150 (MQ=255)
                                                         aaaaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACg                                                                                              >  1:5950460/1‑150 (MQ=255)
                                                         aaaaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACg                                                                                              >  2:6629775/1‑150 (MQ=255)
                                                              cccTTGCTGCTGGGGCTGGCGCATGAATGCCATAAGGTCGAGCGGCTGGCGCAGGCCAGTTGTTATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGCCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTtgtt                                                                                         <  2:3083169/150‑2 (MQ=255)
                                                               ccTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtgg                                                                                        >  2:4864423/1‑150 (MQ=255)
                                                                 ttGCTGCTGGGGCTGGCGCATGACGGCCAGAAGGTCGCGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggct                                                                                      <  2:6368132/150‑1 (MQ=255)
                                                                 ttGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggct                                                                                      >  1:6469701/1‑150 (MQ=255)
                                                                 ttGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggct                                                                                      >  1:3680654/1‑150 (MQ=255)
                                                                 ttGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggct                                                                                      >  1:7107107/1‑150 (MQ=255)
                                                                    ctgctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggc                                                                                   >  2:5654517/1‑150 (MQ=255)
                                                                      gctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctg                                                                                 >  1:6979325/1‑150 (MQ=255)
                                                                        tgGGGCTGGCGCTTGAAGGCCAGATGGTCGCGCGGCTGGCGCTGGCCTGTTGGGTTGTCCCGCTGCCGGGCCCGGTCTGGGACGGTGGGTGGTACCGGGCGCCAGGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctggc                                                                               <  1:4363515/150‑1 (MQ=255)
                                                                        tgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctggc                                                                               >  1:4127918/1‑150 (MQ=255)
                                                                         ggggCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCg                                                                              >  2:3621422/1‑150 (MQ=255)
                                                                          gggCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGt                                                                             >  2:1816493/1‑150 (MQ=255)
                                                                           ggCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCGGGCGCCCCGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTc                                                                            <  2:3773040/150‑1 (MQ=255)
                                                                             cTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGa                                                                          <  2:4050831/150‑1 (MQ=255)
                                                                              gggCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGAc                                                                         <  2:731141/149‑1 (MQ=255)
                                                                                  gcATGAATGCCAGATGGTCGCGCGGCTGGCGCTGGCCTGTTGGTTTGTTCCGCTGCCGGGCACGGTCTCGGACGGTGGCTGGTACCTGGCGCCACGTTGCGGGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGc                                                                     <  1:1014802/150‑1 (MQ=255)
                                                                                   cATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCa                                                                    <  2:4238804/150‑1 (MQ=255)
                                                                                     tGAATGCCCGACGGTCGAGCGGCTGGCGCAGGCCTGTTGGGATGTTCCGCTGCTGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACTTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATc                                                                  <  2:3464640/150‑1 (MQ=255)
                                                                                           ccAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAg                                                            >  2:317255/1‑150 (MQ=255)
                                                                                            cAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAgg                                                           >  1:4543857/1‑150 (MQ=255)
                                                                                              gAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGtt                                                                                                                         <  1:2620786/86‑1 (MQ=255)
                                                                                              gAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGtt                                                                                                                         >  2:2620786/1‑86 (MQ=255)
                                                                                               aaGGTCGCGCGGCGGGCGCCGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGt                                                        <  2:95009/150‑1 (MQ=255)
                                                                                                aGGTCGAGCGGCTGGCGCAGGCCTGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGtt                                                       <  1:6630741/150‑1 (MQ=255)
                                                                                                aGGTCGAGCGGCTGGCGCAGGCCATTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGTTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGtt                                                       <  1:6753735/150‑1 (MQ=255)
                                                                                                 ggTCGAGCGGCTGGCGCTGGCCTGTTGGGATGTTCGGCTGCTGGGCAGGTTCTCGGACCGTGGGTGGTACCTGGCGCCACGTTGAGGGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTg                                                      <  1:1131124/150‑1 (MQ=255)
                                                                                                  gTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGt                                                     >  1:3996609/1‑150 (MQ=255)
                                                                                                  gTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGt                                                     <  1:4255288/150‑1 (MQ=255)
                                                                                                   tCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTc                                                    <  2:5182297/150‑1 (MQ=255)
                                                                                                     gAGCGGCTGGCGCAGGCCTGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCa                                                  <  1:2972654/150‑1 (MQ=255)
                                                                                                       gCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGc                                                >  1:1938672/1‑150 (MQ=255)
                                                                                                        cGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCt                                               <  2:3204647/150‑1 (MQ=255)
                                                                                                          ggtGGCGCCGGCCTGTTGGGATGTTCCGCTGCTGGGCAGGGTCTGGGACGGGGGCTGTTCCCTGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTg                                             <  1:1357262/148‑1 (MQ=255)
                                                                                                           cTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgc                                            >  2:2329573/1‑150 (MQ=255)
                                                                                                             ggCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCGGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgcgc                                          <  2:1472933/150‑1 (MQ=255)
                                                                                                              gcgcAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCt                                         >  1:3198330/1‑150 (MQ=255)
                                                                                                                 ctggCCTGTTGGGATGTTCCGCTGCTGGGCACGGTCTCGGCCCGTGGCTGGTACCGGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTccc                                      <  1:5232573/148‑1 (MQ=255)
                                                                                                                 ccggCCTGTTGGTTTGTTCCGCGGCCGGGCCCGGTCTGGGACGGTGGGTGGTACCGGGCGCAAGTTTGCGGGGCCCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTccc                                      <  1:3101441/148‑1 (MQ=255)
                                                                                                                  aGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCa                                     >  1:1413371/1‑150 (MQ=255)
                                                                                                                    gCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTCGCGCTCCCAGa                                   >  2:3918923/1‑150 (MQ=255)
                                                                                                                     cctgTTGGTATGTTCCGGTGCAGGGCACGGTATCGGACGGTGGCTGGTACCGGGCGCCACTTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGaa                                  <  1:3799736/147‑1 (MQ=255)
                                                                                                                            gtaTGTTCCGCTGCAGGGCACGGTCTCGGACGGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           <  1:5811977/148‑1 (MQ=255)
                                                                                                                            ggATGTTCCGCTGCAGGGCACGGTCTCGTACCGTGGCTGTTACCGGGCGCCACTTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           <  1:3459556/150‑1 (MQ=255)
                                                                                                                            ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           <  2:4383869/150‑1 (MQ=255)
                                                                                                                               tGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCa                        <  1:7400698/150‑1 (MQ=255)
                                                                                                                               tGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCa                        >  2:1294080/1‑150 (MQ=255)
                                                                                                                                       tGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGt                >  1:62990/1‑150 (MQ=255)
                                                                                                                                       tGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGt                >  1:2683695/1‑150 (MQ=255)
                                                                                                                                        gctgggCAGGTTCTCGGACGGGGGCTGGTACCTGGCGCCAGGTTGCGGGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa               <  1:3280362/147‑1 (MQ=255)
                                                                                                                                            ggCAGGGTCTGGGACGGTGGCTGGTACGGGGCGCCAGGTTGTGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGa           <  1:7051984/150‑1 (MQ=255)
                                                                                                                                            ggCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGGCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGa           <  2:4454840/150‑1 (MQ=255)
                                                                                                                                            gccACGGTCCGAGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGa           >  1:6844398/3‑150 (MQ=255)
                                                                                                                                                   tctcGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACAcc    <  2:2257204/150‑1 (MQ=255)
                                                                                                                                                     tcGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCa  >  2:2148669/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCA  >  NC_002947/5940412‑5940710

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: