Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,937,536 G→A 6.8% S193S (TCG→TCA) ‡ fau → 5‑formyltetrahydrofolate cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,937,5360GA 6.8% 402.7 / 12.3 132S193S (TCG→TCA) ‡fau5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base G (57/66);  new base A (5/4);  total (62/70)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

GGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAGG  >  NC_002947/5937387‑5937685
                                                                                                                                                     |                                                                                                                                                     
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:4524455/1‑150 (MQ=255)
         ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                              <  2:4100325/150‑1 (MQ=255)
         ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                              <  2:2700292/150‑3 (MQ=255)
         gtCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAACCCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                              <  2:2700351/149‑3 (MQ=255)
          tCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCt                                                                                                                                             <  2:983508/150‑1 (MQ=255)
          tCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCt                                                                                                                                             <  2:983493/150‑1 (MQ=255)
            tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                           >  1:3691324/1‑150 (MQ=255)
            tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                           >  2:2427607/1‑150 (MQ=255)
            tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                           >  1:4198452/1‑150 (MQ=255)
              tGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTa                                                                                                                                         >  2:4694737/1‑150 (MQ=255)
              tGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTa                                                                                                                                         >  2:4520753/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)
                aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                       >  1:3570290/1‑150 (MQ=255)
                 tCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCACAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCt                                                                                                                                      >  1:432307/1‑150 (MQ=255)
                    cAGCCTGGCCTACGAGGCGCGTCGGCAAACCTGGAAAAATCCCTTGCTGCTGGGGCTGGCGCATGTTTGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                   <  2:2120382/150‑1 (MQ=255)
                    cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                   <  1:324050/150‑1 (MQ=255)
                    cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                   <  2:1791719/150‑1 (MQ=255)
                    cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                   <  2:5086427/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                                                                                                                         >  2:3467020/1‑150 (MQ=255)
                              tACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGa                                                                                                                         <  1:2751969/150‑1 (MQ=255)
                                 cAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgcg                                                                                                                      <  1:341657/150‑1 (MQ=255)
                                     cgcgTCGGCAAACCTGGAAAAATCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGgcgc                                                                                                                  <  2:3219463/150‑1 (MQ=255)
                                            gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                           <  2:975368/150‑1 (MQ=255)
                                             cAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                          <  2:2443312/150‑1 (MQ=255)
                                             cAAACCTGGAAAAAACCCTAGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                          <  2:1674546/150‑1 (MQ=255)
                                                  cTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCTGGGCACGGTCTCGGACCGTGGCGGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTAt                                                                                                     >  2:2110188/1‑150 (MQ=255)
                                                  cTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTAt                                                                                                     >  2:928476/1‑150 (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)
                                                         aaaCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGGGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGtt                                                                                              <  2:4661866/150‑1 (MQ=255)
                                                           aCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGACGGTCGAGCGGCTGGCGCAGGCCAGTGGGGATGTGCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTtgt                                                                                            <  2:4488920/150‑1 (MQ=255)
                                                           aCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTtgt                                                                                            <  1:4694737/150‑1 (MQ=255)
                                                           aCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCACCAGCGGGGGGTATCAACGTtgt                                                                                            <  2:3680182/150‑1 (MQ=255)
                                                             ccTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGAACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtgg                                                                                          >  2:4157361/1‑150 (MQ=255)
                                                               ttGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                          <  1:73869/132‑1 (MQ=255)
                                                               ttGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                          >  2:73869/1‑132 (MQ=255)
                                                                   tgctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggct                                                                                    >  1:109995/1‑150 (MQ=255)
                                                                   tgctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggct                                                                                    <  2:3882265/150‑1 (MQ=255)
                                                                      gggggCTGGCGCTTTACTGCCAGAAGGTCGCGGGGCTGGCGCCGGCCTGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACGGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctggc                                                                                 <  1:4501451/149‑1 (MQ=255)
                                                                                gcATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGc                                                                       >  1:596339/1‑150 (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:2567335/1‑150 (MQ=255)
                                                                                          cAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAgg                                                             >  1:2694748/1‑150 (MQ=255)
                                                                                           aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                            <  2:2694748/150‑1 (MQ=255)
                                                                                           aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                            >  2:285683/1‑150 (MQ=255)
                                                                                             aaGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGt                                                          <  2:1046268/150‑1 (MQ=255)
                                                                                             aaGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGt                                                          >  2:4495405/1‑150 (MQ=255)
                                                                                                 tCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTc                                                      <  1:5090083/150‑1 (MQ=255)
                                                                                                       ggCTGGCGCAGGCCTGTTGGGTTGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCtt                                                <  2:434695/150‑1 (MQ=255)
                                                                                                         cTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgc                                              >  2:2398112/1‑150 (MQ=255)
                                                                                                                aGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCa                                       >  2:4759896/1‑150 (MQ=255)
                                                                                                                  gCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGa                                     >  1:3174744/1‑150 (MQ=255)
                                                                                                                   ccAGTTGGTATGTTCCGCTGCTGGGCAGGGTCTGGGACCGTGGCTGGTACCGGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGaa                                    <  1:81197/150‑1 (MQ=255)
                                                                                                                    cAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                                 >  1:1108043/1‑120 (MQ=255)
                                                                                                                    cAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAg                                   >  1:4001965/1‑150 (MQ=255)
                                                                                                                       ttGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCt                                >  2:4646338/1‑150 (MQ=255)
                                                                                                                         gggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGAGAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTg                              >  2:3397271/1‑150 (MQ=255)
                                                                                                                          gtttgtgCCGGTGCTGGGCACGGTCTCGGACGGTGGGTGGTACGTGGCGCCAGGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                             <  1:200707/147‑1 (MQ=255)
                                                                                                                          ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                             <  2:1210742/150‑1 (MQ=255)
                                                                                                                          ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGGTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                             >  1:1668688/1‑150 (MQ=255)
                                                                                                                               ttCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATc                                                                    >  1:3347443/1‑106 (MQ=255)
                                                                                                                               ttCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATc                                                                    <  2:3347443/106‑1 (MQ=255)
                                                                                                                                    cTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCgg                   >  1:3691786/1‑150 (MQ=255)
                                                                                                                                     taccagccACGGTCCGAGACCGTGCCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGt                  >  2:93064/8‑150 (MQ=255)
                                                                                                                                     taccagccACGGTCCGAGACCGTGCCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGt                  >  2:93303/8‑150 (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)
                                                                                                                                          ggCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGa             >  1:4960512/1‑150 (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)
                                                                                                                                                   tcGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCa    >  1:1774440/1‑150 (MQ=255)
                                                                                                                                                    cGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAg   <  2:2586094/150‑1 (MQ=255)
                                                                                                                                                    cGGACCGGGGCTGGTACCGGGCGCCACGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAg   <  1:196081/150‑1 (MQ=255)
                                                                                                                                                     ggACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTc                                                      <  2:548900/98‑1 (MQ=255)
                                                                                                                                                     ggACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTc                                                      >  1:548900/1‑98 (MQ=255)
                                                                                                                                                     ggACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAgg  <  2:4198452/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCAGG  >  NC_002947/5937387‑5937685

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: