Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,937,535 C→G 7.4% 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,937,5350CG 7.4% 422.4 / 17.6 132S193W (TCG→TGG) ‡fau5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base C (56/66);  new base G (5/5);  total (61/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 139 total aligned reads displayed.

GGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAG  >  NC_002947/5937386‑5937684
                                                                                                                                                     |                                                                                                                                                     
gggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                       <  1:1486775/150‑1 (MQ=255)
 ggCGGTGGCTTCTATGGTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCg                                                                                                                                                      <  1:4268526/150‑1 (MQ=255)
 ggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAGCCCTTGCTGCTGGGGCTGGCGCGTGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCg                                                                                                                                                      <  2:967029/150‑1 (MQ=255)
 ggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAGCCCTTGCTGCTGGGGCTGGCGCGTGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCg                                                                                                                                                      <  2:4117903/150‑1 (MQ=255)
 ggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCg                                                                                                                                                      <  2:3901925/150‑1 (MQ=255)
         cTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgc                                                                                                                                              <  1:1672392/150‑2 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                             >  1:1046268/1‑150 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                             >  1:4200029/1‑150 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                             >  1:4984487/1‑150 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                             <  2:4100325/150‑1 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                             <  2:3399714/150‑1 (MQ=255)
          ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                             >  1:1821004/1‑148 (MQ=255)
          gtCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAACCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                             <  2:2700351/149‑3 (MQ=255)
           tCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCt                                                                                                                                            <  2:983508/150‑1 (MQ=255)
           tCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCt                                                                                                                                            <  2:983493/150‑1 (MQ=255)
             tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                          >  2:2427607/1‑150 (MQ=255)
             tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                          >  1:3691324/1‑150 (MQ=255)
               tGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTa                                                                                                                                        >  2:4694737/1‑150 (MQ=255)
                gATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAc                                                                                                                                       <  2:3888307/150‑1 (MQ=255)
                 aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                      >  1:998523/1‑150 (MQ=255)
                 aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                      <  1:4503442/150‑1 (MQ=255)
                 aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                      >  1:2624359/1‑150 (MQ=255)
                  tCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCACAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCt                                                                                                                                     >  1:432307/1‑150 (MQ=255)
                     cAGCCTGGCCTACGAGGCGCGTCGGCAAACCTGGAAAAATCCCTTGCTGCTGGGGCTGGCGCATGTTTGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                  <  2:2120382/150‑1 (MQ=255)
                     cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                  <  2:1791719/150‑1 (MQ=255)
                     cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                  <  1:324050/150‑1 (MQ=255)
                     cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                  <  2:4288262/150‑1 (MQ=255)
                      aGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgcg                                                                                                                                 >  2:196081/1‑150 (MQ=255)
                         cTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCa                                                                                                                              <  1:3416306/150‑1 (MQ=255)
                          tGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCAc                                                                                                                             >  2:2843072/1‑150 (MQ=255)
                           ggCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACg                                                                                                                            >  2:1569883/1‑150 (MQ=255)
                              cTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTg                                                                                                                         >  1:4623283/1‑150 (MQ=255)
                               tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                        <  1:2751969/150‑1 (MQ=255)
                               tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                        >  2:3467020/1‑150 (MQ=255)
                                  cAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgcg                                                                                                                     <  1:341657/150‑1 (MQ=255)
                                      cgcgTCGGCAAACCTGGAAAAATCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGgcgc                                                                                                                 <  2:3219463/150‑1 (MQ=255)
                                             gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                          <  2:975368/150‑1 (MQ=255)
                                             gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                          >  1:881917/1‑150 (MQ=255)
                                              cAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                         <  1:2723766/150‑1 (MQ=255)
                                                   cTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTAt                                                                                                    >  2:928476/1‑150 (MQ=255)
                                                    tGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                   <  2:5023748/150‑1 (MQ=255)
                                                    tGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGCCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                   >  2:1776887/1‑150 (MQ=255)
                                                    tGGAAAAAACCCTAGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                   >  1:599314/1‑150 (MQ=255)
                                                     ggAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCTAGCGGGGGGTATCa                                                                                                  >  2:2155907/1‑150 (MQ=255)
                                                      gAAAAAACCCGTGCGGCGGGGCCTGGCGCATGAATGCCAGAAGGTCGACCGCCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCCCGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCaa                                                                                                 <  2:1370858/150‑1 (MQ=255)
                                                          aaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGtt                                                                                             <  2:4661872/150‑1 (MQ=255)
                                                          aaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGGGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGtt                                                                                             <  2:4661866/150‑1 (MQ=255)
                                                            aCCCTTGCTGCTGGGGCTGGCGCATGAATGCCTGATGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTtgt                                                                                           <  2:4489147/150‑1 (MQ=255)
                                                            aCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGACGGTCGAGCGGCTGGCGCAGGCCAGTGGGGATGTGCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTtgt                                                                                           <  2:4488920/150‑1 (MQ=255)
                                                            aCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTtgt                                                                                           <  1:4694737/150‑1 (MQ=255)
                                                              ccTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGAACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtgg                                                                                         >  2:4157361/1‑150 (MQ=255)
                                                                ttGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                         >  2:73869/1‑132 (MQ=255)
                                                                    tgctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggct                                                                                   <  2:3882265/150‑1 (MQ=255)
                                                                    tgctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggct                                                                                   >  1:109995/1‑150 (MQ=255)
                                                                      ctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctgg                                                                                 >  1:720668/1‑150 (MQ=255)
                                                                       gggggCTGGCGCTTTACTGCCAGAAGGTCGCGGGGCTGGCGCCGGCCTGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACGGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctggc                                                                                <  1:4501451/149‑1 (MQ=255)
                                                                             tGGCTCATGAATGCCAGAAGTTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGTTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGAc                                                                          <  2:956966/150‑1 (MQ=255)
                                                                                 gcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGc                                                                      >  1:596339/1‑150 (MQ=255)
                                                                                 gcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGc                                                                      <  2:2431320/150‑1 (MQ=255)
                                                                                     gAATGCCAGAAGGTCGAGCGGCTGGCGCTGGCCTTTTGGGTTGTCCCGCTGCAGGGCCCGGTCTCGGACGGTGGGTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCt                                                                  <  1:1145477/150‑1 (MQ=255)
                                                                                     gAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCt                                                                  <  2:4623283/150‑1 (MQ=255)
                                                                                       aTGCCAGACGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCGGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                                <  2:77277/150‑1 (MQ=255)
                                                                                       aTGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                                >  2:1410928/1‑150 (MQ=255)
                                                                                         gCCATAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATa                                                              <  1:1184527/150‑1 (MQ=255)
                                                                                         gCCAGACGGTCGAGCGGGTGGCGCAGGCCTGTTGGGATGTTCCGCTGCTGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACTTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATa                                                              <  2:50323/150‑1 (MQ=255)
                                                                                         gCCAGAAGGTCGAGCGGCGGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATa                                                              <  2:3861561/150‑1 (MQ=255)
                                                                                           cAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAgg                                                            >  1:3682172/1‑150 (MQ=255)
                                                                                           cAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAgg                                                            >  1:2567335/1‑150 (MQ=255)
                                                                                            aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                           >  2:1577477/1‑150 (MQ=255)
                                                                                            aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                           <  2:2694748/150‑1 (MQ=255)
                                                                                             gAAGGTCGAGCGGCTGGCGCAGGCCTGTTGGGATGGTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCg                                                          <  2:2979501/150‑1 (MQ=255)
                                                                                              aaGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGt                                                         <  2:1046268/150‑1 (MQ=255)
                                                                                              aaGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGt                                                         >  2:4495405/1‑150 (MQ=255)
                                                                                                  tCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTc                                                     <  1:5090083/150‑1 (MQ=255)
                                                                                                          cTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgc                                             >  2:2398112/1‑150 (MQ=255)
                                                                                                                   gCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGa                                    >  1:3174744/1‑150 (MQ=255)
                                                                                                                    ccAGTTGGTATGTTCCGCTGCTGGGCAGGGTCTGGGACCGTGGCTGGTACCGGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGaa                                   <  1:81197/150‑1 (MQ=255)
                                                                                                                     cAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                                <  2:1108043/120‑1 (MQ=255)
                                                                                                                     cAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                                >  1:1108043/1‑120 (MQ=255)
                                                                                                                     cAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAg                                  >  1:4001965/1‑150 (MQ=255)
                                                                                                                      aGTTGGGAGGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGc                                 <  1:4567577/150‑1 (MQ=255)
                                                                                                                           gtttgtgCCGGTGCTGGGCACGGTCTCGGACGGTGGGTGGTACGTGGCGCCAGGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                            <  1:200707/147‑1 (MQ=255)
                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                            <  2:2631221/150‑1 (MQ=255)
                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                            <  2:599314/150‑1 (MQ=255)
                                                                                                                           ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGGTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                            >  1:1668688/1‑150 (MQ=255)
                                                                                                                                ttCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATc                                                                   <  2:3347443/106‑1 (MQ=255)
                                                                                                                                ttCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATc                                                                   >  1:3347443/1‑106 (MQ=255)
                                                                                                                                     cTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCgg                  >  1:3691786/1‑150 (MQ=255)
                                                                                                                                      taccagccACGGTCCGAGACCGTGCCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGt                 >  2:93064/8‑150 (MQ=255)
                                                                                                                                       gCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa                <  2:4614950/150‑1 (MQ=255)
                                                                                                                                       gCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTa                <  2:1511684/150‑1 (MQ=255)
                                                                                                                                        ccagccACGGTCCGAGACCGTGCCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAg               >  1:5025089/6‑150 (MQ=255)
                                                                                                                                         cagccACGGTCCGAGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGt              >  1:4064280/5‑150 (MQ=255)
                                                                                                                                           ggCCCGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGa            <  1:2795070/150‑1 (MQ=255)
                                                                                                                                            gCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc           <  2:1938491/150‑1 (MQ=255)
                                                                                                                                            gCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGAc           <  2:1938478/150‑1 (MQ=255)
                                                                                                                                                 gTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCacac      >  1:4897636/1‑150 (MQ=255)
                                                                                                                                                    tcGGACGGTGGCTGGTACCTGGCGCCTCGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCa   <  1:4812457/150‑1 (MQ=255)
                                                                                                                                                     cGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAg  <  2:2586094/150‑1 (MQ=255)
                                                                                                                                                     cGGACCGGGGCTGGTACCGGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAg  <  1:196081/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAG  >  NC_002947/5937386‑5937684

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