Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,685,527 A→G 100% G354G (GGT→GGC fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,5270AG100.0% 318.0 / NA 85G354G (GGT→GGCfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (41/44);  total (41/44)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685381‑1685668
                                                                                                                                                  |                                                                                                                                             
ttGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGcc                                                                                                                                             >  4:337320/1‑149 (MQ=255)
    gTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGaa                                                                                                                                         <  1:396108/149‑1 (MQ=255)
    gTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGaa                                                                                                                                         <  1:484149/149‑1 (MQ=255)
    gTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGaa                                                                                                                                         <  1:625230/149‑1 (MQ=255)
      aaTTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAgc                                                                                                                                       >  2:486456/1‑149 (MQ=255)
      aaTTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAgc                                                                                                                                       >  2:269170/1‑149 (MQ=255)
        ttgattgaTTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGTATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAgcgc                                                                                                                                     >  2:74828/1‑149 (MQ=255)
           attgattCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGa                                                                                                                                          >  3:34434/1‑141 (MQ=255)
           attgattCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGa                                                                                                                                          <  4:34434/141‑1 (MQ=255)
               attCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCAt                                                                                                                              <  1:633567/149‑1 (MQ=255)
                ttCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCAt                                                                                                                              <  3:166962/148‑1 (MQ=255)
                ttCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCAt                                                                                                                              >  4:166962/1‑148 (MQ=255)
                  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGtt                                                                                                                           <  3:351922/149‑1 (MQ=255)
                  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGtt                                                                                                                           <  2:353411/149‑1 (MQ=255)
                    cTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGa                                                                                                                         >  4:429934/1‑149 (MQ=255)
                     tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAGAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCACCCATGTTGAt                                                                                                                        <  3:664739/149‑1 (MQ=255)
                      cACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATGTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATc                                                                                                                       <  1:554266/149‑1 (MQ=255)
                               gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                              >  3:684681/1‑149 (MQ=255)
                               gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                              >  1:329475/1‑149 (MQ=255)
                                ttCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcgt                                                                                                             >  3:596365/1‑149 (MQ=255)
                                ttCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcgt                                                                                                             >  4:191474/1‑149 (MQ=255)
                                   aaTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAGTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTc                                                                                                           >  3:280074/1‑148 (MQ=255)
                                   aaTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTc                                                                                                           <  4:280074/148‑1 (MQ=255)
                                      aCCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                       <  3:51538/149‑1 (MQ=255)
                                      aCCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                       <  4:460802/149‑1 (MQ=255)
                                      aCCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                       >  4:51538/1‑149 (MQ=255)
                                       cccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                      >  2:6981/1‑149 (MQ=255)
                                           cTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGcccccc                                                                                                                                <  3:348864/119‑1 (MQ=255)
                                           cTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGcccccc                                                                                                                                >  4:348864/1‑119 (MQ=255)
                                               gcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                                                              <  1:681343/149‑1 (MQ=255)
                                                cAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAca                                                                                             <  1:348828/149‑1 (MQ=255)
                                                     tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCAGCTGACAGCaga                                                                                        <  3:429934/149‑1 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATc                                                                                                                       <  4:268253/114‑1 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATc                                                                                                                       >  3:268253/1‑114 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                    >  2:592713/1‑149 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                    >  2:116144/1‑149 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                    >  1:532398/1‑149 (MQ=255)
                                                           aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCa                                                                                                                               >  4:6090/1‑104 (MQ=255)
                                                           aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCa                                                                                                                               <  3:6090/104‑1 (MQ=255)
                                                           aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                                                                  <  2:33103/149‑1 (MQ=255)
                                                           aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                                                                  >  1:352766/1‑149 (MQ=255)
                                                             aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGtt                                                                                <  4:74279/149‑1 (MQ=255)
                                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                              <  3:191474/149‑1 (MQ=255)
                                                                tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                                                             <  2:461316/149‑1 (MQ=255)
                                                                tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                                                             <  4:416753/149‑1 (MQ=255)
                                                                 ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATg                                                                            >  3:660395/1‑149 (MQ=255)
                                                                  tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                           <  1:626876/149‑1 (MQ=255)
                                                                    cATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                         <  2:722696/149‑1 (MQ=255)
                                                                     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                        <  1:428391/149‑1 (MQ=255)
                                                                           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCAATGGCCGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtg                                                                  >  4:91569/1‑149 (MQ=255)
                                                                               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                              <  1:180687/149‑1 (MQ=255)
                                                                                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTg                                                             <  4:391084/149‑1 (MQ=255)
                                                                                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTggg                                                           >  3:592090/1‑149 (MQ=255)
                                                                                      gcgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCa                                                       <  1:486456/149‑1 (MQ=255)
                                                                                      gcgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCa                                                       <  1:269170/149‑1 (MQ=255)
                                                                                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc                                                      <  4:592090/149‑1 (MQ=255)
                                                                                          aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGAt                                                                                                                        <  3:218866/80‑1 (MQ=255)
                                                                                          aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGAt                                                                                                                        >  4:218866/1‑80 (MQ=255)
                                                                                          aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACttt                                                   >  1:59584/1‑149 (MQ=255)
                                                                                                tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                             <  4:650703/149‑1 (MQ=255)
                                                                                                   aGGAAATTGTGGATCACCATTGGCCGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGCTCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCAt                                          >  2:545990/1‑149 (MQ=255)
                                                                                                   aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                           <  1:604767/116‑1 (MQ=255)
                                                                                                   aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                           >  2:604767/1‑116 (MQ=255)
                                                                                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        <  4:596365/149‑1 (MQ=255)
                                                                                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        >  1:288822/1‑149 (MQ=255)
                                                                                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        >  3:272991/1‑149 (MQ=255)
                                                                                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        >  4:154736/1‑149 (MQ=255)
                                                                                                       aaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGAtt                                      <  3:72202/149‑1 (MQ=255)
                                                                                                             ggATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTg                                >  2:654127/1‑149 (MQ=255)
                                                                                                             ggATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTg                                >  4:149010/1‑149 (MQ=255)
                                                                                                              gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGc                               >  4:251975/1‑149 (MQ=255)
                                                                                                               aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                              >  4:368161/1‑149 (MQ=255)
                                                                                                                tCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGATGGGTTCACTTTCCCCGGCATGATTGAGCTGccc                             >  3:53971/1‑149 (MQ=255)
                                                                                                                    cATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                         >  4:509442/1‑101 (MQ=255)
                                                                                                                    cATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                         <  3:509442/101‑1 (MQ=255)
                                                                                                                            ggAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTcc                 <  1:654127/149‑1 (MQ=255)
                                                                                                                            ggAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTcc                 >  3:198470/1‑149 (MQ=255)
                                                                                                                              aaGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGATCATTTTCCgg               <  1:573563/149‑1 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga         >  2:106602/1‑149 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga         >  4:431842/1‑149 (MQ=255)
                                                                                                                                      cagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTgaga       <  3:149010/149‑1 (MQ=255)
                                                                                                                                      cagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTgaga       <  3:305289/149‑1 (MQ=255)
                                                                                                                                         ttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAttt    >  1:493832/1‑149 (MQ=255)
                                                                                                                                          tcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTc   <  3:258506/149‑1 (MQ=255)
                                                                                                                                           caAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCa  <  1:592713/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685381‑1685668

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