Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,940,562 C→G 5.6% S193W (TCG→TGG)  fau → 5‑formyltetrahydrofolate cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,940,5620CG 5.6% 1055.5 / 30.1 329S193W (TCG→TGG) fau5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base C (149/157);  new base G (8/12);  total (158/171)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

GGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAG  >  NC_002947/5940413‑5940711
                                                                                                                                                     |                                                                                                                                                     
 ggCGGTGGCCTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCg                                                                                                                                                      <  1:5850949/150‑1 (MQ=255)
         cTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTgg                                                                                                                                              <  1:3594616/150‑1 (MQ=255)
         cTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgc                                                                                                                                              <  2:8770684/150‑2 (MQ=255)
         cTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgc                                                                                                                                              <  1:6151435/150‑2 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                             >  1:7393017/1‑150 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                             <  2:1595549/150‑3 (MQ=255)
           tCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgccc                                                                                                                                            <  1:7240971/150‑4 (MQ=255)
            cTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTg                                                                                                                                           <  1:4657636/150‑1 (MQ=255)
            cTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTg                                                                                                                                           <  2:439411/150‑1 (MQ=255)
             tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                          >  1:7787167/1‑150 (MQ=255)
             tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                          >  1:4717061/1‑150 (MQ=255)
                 aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                      >  1:7323522/1‑150 (MQ=255)
                  tCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCt                                                                                                                                     >  1:6224683/1‑150 (MQ=255)
                  tCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCt                                                                                                                                     >  2:6858353/1‑150 (MQ=255)
                  tCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCt                                                                                                                                     >  1:3825463/1‑150 (MQ=255)
                    gCAGCCTGGCCTACCAGGCGCGTCTGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTgg                                                                                                                                   >  2:2411126/1‑150 (MQ=255)
                    gCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTgg                                                                                                                                   <  1:951783/150‑1 (MQ=255)
                       gCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgcgc                                                                                                                                <  1:8773700/150‑1 (MQ=255)
                        ccTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGcc                                                                                                                               >  2:1279288/1‑150 (MQ=255)
                        ccTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGcc                                                                                                                               >  2:6424270/1‑150 (MQ=255)
                              cTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTg                                                                                                                         <  2:2843826/150‑1 (MQ=255)
                              cTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTg                                                                                                                         <  1:2005992/150‑1 (MQ=255)
                               tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                        >  1:3082074/1‑150 (MQ=255)
                               tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                        >  2:7459/1‑150 (MQ=255)
                               tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                        >  1:265230/1‑150 (MQ=255)
                                 ccAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgc                                                                                                                      >  1:7158581/1‑150 (MQ=255)
                                 ccAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgc                                                                                                                      <  1:3428947/150‑1 (MQ=255)
                                      cgcgTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGgcgc                                                                                                                 <  1:5443585/150‑1 (MQ=255)
                                            ggCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgg                                                                                                           >  2:5477556/1‑150 (MQ=255)
                                            ggCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgg                                                                                                           <  1:8429086/150‑1 (MQ=255)
                                             gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                          >  1:7637108/1‑150 (MQ=255)
                                              cAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                         <  2:7734980/150‑1 (MQ=255)
                                               aaaCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggggg                                                                                                        <  1:4432328/150‑1 (MQ=255)
                                               aaaCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggggg                                                                                                        <  2:1903196/150‑1 (MQ=255)
                                               aaaCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTAGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggggg                                                                                                        <  2:5670636/150‑1 (MQ=255)
                                                    tGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                   >  2:4753350/1‑150 (MQ=255)
                                                     ggAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCa                                                                                                  >  2:823378/1‑150 (MQ=255)
                                                     ggAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGGTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCCGCGGGGGGTCTCa                                                                                                  >  2:7421051/1‑150 (MQ=255)
                                                       aaaaaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAAc                                                                                                >  2:1752936/1‑150 (MQ=255)
                                                         aaaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGt                                                                                              >  2:3387279/1‑150 (MQ=255)
                                                         aaaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGAAGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGt                                                                                              >  2:8018692/1‑150 (MQ=255)
                                                             cccTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTtgtg                                                                                          >  1:1457138/1‑150 (MQ=255)
                                                              ccgtgctgcGGGGGGTGGCGCATGACTGCCAGACGGTCGAGCGGCGGGCGCCGGCCAGTTGGGATGTTCCGCTGCAGGGCCCGGTCTCGGACCGGGGCTGGTACCGGGCGCCACGTGGAGCGGCGCCAGCGGGGGGTATCAACGTTGtgg                                                                                         <  2:6168940/147‑1 (MQ=255)
                                                                      ctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTACAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctgg                                                                                 >  1:2548730/1‑150 (MQ=255)
                                                                         gggCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGt                                                                              >  1:6469551/1‑150 (MQ=255)
                                                                              ggCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGTCTGGCGTCGAcc                                                                         <  2:2123716/150‑1 (MQ=255)
                                                                              ggCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGAcc                                                                         <  2:3267001/150‑1 (MQ=255)
                                                                               gcgcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCg                                                                        >  1:6847692/1‑150 (MQ=255)
                                                                                 gcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGc                                                                      >  1:4071343/1‑150 (MQ=255)
                                                                                 gcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGc                                                                      <  2:1503913/150‑1 (MQ=255)
                                                                                   aTGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCAt                                                                    >  1:2050680/1‑150 (MQ=255)
                                                                                    tGAATGCCAGACGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATc                                                                   <  2:388384/150‑1 (MQ=255)
                                                                                    tGAATGCCAGAAGGTCGCGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCATGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATc                                                                   <  2:2967885/150‑1 (MQ=255)
                                                                                     gAATGCCAGACGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCt                                                                  <  2:1182617/150‑1 (MQ=255)
                                                                                      aaTGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTg                                                                 >  2:6683466/1‑150 (MQ=255)
                                                                                        ggCCAGATGGTCGCGGGGGTGGCGCTGGCCTGTTGGGTTGTTCCGCTGCTGGGCACGGTCTCGGACGGTGGCTGGTACCGGGCGCCAGGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGAt                                                               <  1:8267267/149‑1 (MQ=255)
                                                                                           cAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCt                                                                                                                                            >  1:5973562/1‑70 (MQ=255)
                                                                                           cAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCt                                                                                                                                            <  2:5973562/70‑1 (MQ=255)
                                                                                            aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                           >  2:6981817/1‑150 (MQ=255)
                                                                                            aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                           >  2:6270441/1‑150 (MQ=255)
                                                                                            aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                           <  2:6309875/150‑1 (MQ=255)
                                                                                               aGGTCGTGCGGGTGGCGCAGGCCTGTGGGGATGTGCCGCTGCAGGGCACGGTCTCGGACGGTGGGTGGTACCGGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGtt                                                        <  1:7036124/150‑1 (MQ=255)
                                                                                                ggTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTg                                                       >  1:3526135/1‑150 (MQ=255)
                                                                                                 gTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGt                                                      >  1:107415/1‑150 (MQ=255)
                                                                                                   cGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTcc                                                    <  2:3675516/150‑1 (MQ=255)
                                                                                                    gAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCa                                                   >  1:5231000/1‑150 (MQ=255)
                                                                                                     aGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAg                                                  >  1:5722996/1‑150 (MQ=255)
                                                                                                         gCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTg                                              <  2:5351780/150‑1 (MQ=255)
                                                                                                          cTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgc                                             >  1:4164530/1‑150 (MQ=255)
                                                                                                          cTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgc                                             <  1:5000447/150‑1 (MQ=255)
                                                                                                               gcAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTcc                                        <  2:7787167/150‑1 (MQ=255)
                                                                                                                   gCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGa                                    >  2:5095798/1‑150 (MQ=255)
                                                                                                                      aGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGc                                 >  2:7260322/1‑150 (MQ=255)
                                                                                                                        ttGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCt                               >  2:2532169/1‑150 (MQ=255)
                                                                                                                        ttGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCt                               <  1:8635640/150‑1 (MQ=255)
                                                                                                                         tGGGATGTTCCGCTGCTGGGCAGGGTCTGGGCCGGTGGGTGGTACGGGGCGCCTGGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCtt                              <  1:1161039/150‑1 (MQ=255)
                                                                                                                         tGGGATGTTCCGCTGCAGGGCAGGGTCTCGTACCGTGGCTGGTACCTGGCGCCACTTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCtt                              <  1:395162/150‑1 (MQ=255)
                                                                                                                         tGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCtt                              <  1:4342258/150‑1 (MQ=255)
                                                                                                                          gggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTg                             <  2:3524425/150‑1 (MQ=255)
                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                            <  2:424741/150‑1 (MQ=255)
                                                                                                                                ttCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATa                       >  2:5602381/1‑150 (MQ=255)
                                                                                                                                   cGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAAcc                    >  1:3235871/1‑150 (MQ=255)
                                                                                                                                    gCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCg                   >  1:2743367/1‑150 (MQ=255)
                                                                                                                                      taccagccACGGTCCGAGACCGTGCCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGt                 >  2:421106/8‑150 (MQ=255)
                                                                                                                                       gCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa                <  2:4405200/150‑1 (MQ=255)
                                                                                                                                       gCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa                <  2:8552679/150‑1 (MQ=255)
                                                                                                                                       gCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa                <  2:6492078/150‑1 (MQ=255)
                                                                                                                                         cagccACGGTCCGAGACCGTGGCTGGTACCTGGCGCCACGTCGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGt              >  1:5520991/5‑150 (MQ=255)
                                                                                                                                           ggCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGa            <  1:1242920/150‑1 (MQ=255)
                                                                                                                                            gCACGGTCTCGGACGGTGGCTGGTACCTGGCGCCTCGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc           <  1:5593898/150‑1 (MQ=255)
                                                                                                                                            gCACGGTCTCGGACGGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc           <  1:4920925/150‑1 (MQ=255)
                                                                                                                                            gCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc           <  1:6424175/150‑1 (MQ=255)
                                                                                                                                            gCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc           <  2:3790679/150‑1 (MQ=255)
                                                                                                                                            ccACGGTCCGAGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc           >  1:6564884/2‑150 (MQ=255)
                                                                                                                                             cACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAGCCGGTAGTGAcc          >  2:6996600/1‑150 (MQ=255)
                                                                                                                                             cACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAcc          >  1:1647258/1‑150 (MQ=255)
                                                                                                                                               cGGTCTCGGACCGTGGAGTGTACGTGGCTCCACTTGGAGTGGCGCCAGCGGGGGGTATCAACGTTGTGGCTTGCTGGCGTCGACCGGCATATGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACcac        <  1:643877/150‑1 (MQ=255)
                                                                                                                                                ggTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACcaca       >  2:721509/1‑150 (MQ=255)
                                                                                                                                                    tcGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCa   >  1:2298743/1‑150 (MQ=255)
                                                                                                                                                     cGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAg  <  2:4634697/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAG  >  NC_002947/5940413‑5940711

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