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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,937,5340TC 6.9% 392.1 / 12.9 132S193P (TCG→CCG) ‡fau5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base T (56/65);  new base C (5/4);  total (61/71)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 138 total aligned reads displayed.

TGGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCA  >  NC_002947/5937385‑5937683
                                                                                                                                                     |                                                                                                                                                     
 gggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                      <  1:3332735/150‑1 (MQ=255)
 gggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                      <  1:1486775/150‑1 (MQ=255)
 gggCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTCGGGATGTTCCGCTGCAGGGCACGGtctc                                                                                                                                                      <  1:986260/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:4524455/1‑150 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                            <  2:3399714/150‑1 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                            >  1:1046268/1‑150 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                            <  2:4100325/150‑1 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGc                                                                                                                                            >  1:4200029/1‑150 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                            <  2:1821004/150‑3 (MQ=255)
           ttCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCCGAGACCGTgcc                                                                                                                                            >  1:1821004/1‑148 (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)
             cTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTg                                                                                                                                          <  1:2852999/150‑1 (MQ=255)
              tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                         >  1:2454446/1‑150 (MQ=255)
              tATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTgg                                                                                                                                         >  1:4198452/1‑150 (MQ=255)
                tGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTa                                                                                                                                       >  2:4520753/1‑150 (MQ=255)
                  aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                     >  2:3895516/1‑150 (MQ=255)
                  aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                     >  1:3570290/1‑150 (MQ=255)
                  aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                     >  1:2624359/1‑150 (MQ=255)
                  aTCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTAcc                                                                                                                                     <  1:4503442/150‑1 (MQ=255)
                   tCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCACAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCt                                                                                                                                    >  1:432307/1‑150 (MQ=255)
                      cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                 <  1:324050/150‑1 (MQ=255)
                      cAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgc                                                                                                                                 <  2:1791719/150‑1 (MQ=255)
                       aGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGgcg                                                                                                                                >  2:196081/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)
                                 aCCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAg                                                                                                                      >  1:4804611/1‑150 (MQ=255)
                                   cAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAgcg                                                                                                                    <  1:341657/150‑1 (MQ=255)
                                       cgcgTCGGCAAACCTGGAAAAATCCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGgcgc                                                                                                                <  2:3219463/150‑1 (MQ=255)
                                           tCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAgg                                                                                                            <  1:4497583/150‑2 (MQ=255)
                                              gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                         <  2:975368/150‑1 (MQ=255)
                                              gCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCggg                                                                                                         >  1:881917/1‑150 (MQ=255)
                                               cAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                        <  2:2443312/150‑1 (MQ=255)
                                               cAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCgggg                                                                                                        <  1:2723766/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)
                                                     tGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                  <  2:5023748/150‑1 (MQ=255)
                                                     tGGAAAAAACCCTAGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATc                                                                                                  >  1:599314/1‑150 (MQ=255)
                                                      ggAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCTAGCGGGGGGTATCa                                                                                                 >  2:2155907/1‑150 (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)
                                                               ccTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtgg                                                                                        >  1:540664/1‑150 (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)
                                                                       ctgGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGtggctggctgg                                                                                >  1:720668/1‑150 (MQ=255)
                                                                              tGGCTCATGAATGCCAGAAGTTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGTTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGAc                                                                         <  2:956966/150‑1 (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:2694748/1‑150 (MQ=255)
                                                                                             aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                          >  2:285683/1‑150 (MQ=255)
                                                                                             aGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGc                                                          >  2:1577477/1‑150 (MQ=255)
                                                                                              gAAGGTCGAGCGGCTGGCGCAGGCCTGTTGGGATGGTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCg                                                         <  2:2979501/150‑1 (MQ=255)
                                                                                               aaGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGt                                                        <  2:1046268/150‑1 (MQ=255)
                                                                                                   tCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTc                                                    <  1:5090083/150‑1 (MQ=255)
                                                                                                     gAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCa                                                  >  2:1501751/1‑150 (MQ=255)
                                                                                                         ggCTGGCGCAGGCCTGTTGGGTTGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCtt                                              <  2:434695/150‑1 (MQ=255)
                                                                                                           cTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGTGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgc                                            >  2:2398057/1‑150 (MQ=255)
                                                                                                           cTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTgc                                            >  2:2398112/1‑150 (MQ=255)
                                                                                                                  aGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCa                                     >  2:4759896/1‑150 (MQ=255)
                                                                                                                      cAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                               >  1:1108043/1‑120 (MQ=255)
                                                                                                                      cAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGa                                                               <  2:1108043/120‑1 (MQ=255)
                                                                                                                       aGTTGGGAGGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGc                                <  1:4567577/150‑1 (MQ=255)
                                                                                                                         ttGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCt                              <  1:3731306/150‑1 (MQ=255)
                                                                                                                           gggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGAGAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTg                            >  2:3397271/1‑150 (MQ=255)
                                                                                                                            gtttgtgCCGGTGCTGGGCACGGTCTCGGACGGTGGGTGGTACGTGGCGCCAGGTTGCGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           <  1:200707/147‑1 (MQ=255)
                                                                                                                            ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           <  2:599314/150‑1 (MQ=255)
                                                                                                                            ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           <  2:2631221/150‑1 (MQ=255)
                                                                                                                            ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           <  2:1210742/150‑1 (MQ=255)
                                                                                                                            ggATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGGTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTgc                           >  1:1668688/1‑150 (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)
                                                                                                                                         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:1938478/150‑1 (MQ=255)
                                                                                                                                                  gTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCacac     >  1:4897636/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGCGGTGGCTTCTATGATCGCAGCCTGGCCTACCAGGCGCGTCGGCAAACCTGGAAAAAACCCTTGCTGCTGGGGCTGGCGCATGAATGCCAGAAGGTCGAGCGGCTGGCGCAGGCCAGTTGGGATGTTCCGCTGCAGGGCACGGTCTCGGACCGTGGCTGGTACCTGGCGCCACGTTGAGCGGCGCCAGCGGGGGGTATCAACGTTGTGGCTGGCTGGCGTCGACCGGCATCTGATAGGCGTTGTCCAGCTTGCGCTCCCAGAAGCCTTGCGCATAACCGGTAGTGACCACACCCA  >  NC_002947/5937385‑5937683

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: