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,5270AG98.9% 334.2 / ‑4.5 93G354G (GGT→GGCfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (48/44);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAA  >  NC_000913/1685381‑1685672
                                                                                                                                                  |                                                                                                                                                 
ttGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGcc                                                                                                                                                 <  2:59868/149‑1 (MQ=255)
  gAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCgg                                                                                                                                               >  4:657214/1‑149 (MQ=255)
    gTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGaa                                                                                                                                             <  2:669904/149‑1 (MQ=255)
     tAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAg                                                                                                                                            >  4:602319/1‑149 (MQ=255)
             tgattCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGcccccc                                                                                                                                    <  1:210188/149‑1 (MQ=255)
                ttCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTTCAGGAAATTGTGGATCCCCATTGGACGGAAGACGTTCAGTTCCCAGTTGCCGGAAGCGCCCCCCATg                                                                                                                                 >  4:45748/1‑149 (MQ=255)
                 tCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCATCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGt                                                                                                                                <  1:782658/149‑1 (MQ=255)
                 tCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGt                                                                                                                                <  1:327909/149‑1 (MQ=255)
                     tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGAt                                                                                                                            <  1:544383/149‑1 (MQ=255)
                     tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGAt                                                                                                                            <  1:749679/149‑1 (MQ=255)
                      cACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATc                                                                                                                           <  2:452076/149‑1 (MQ=255)
                      cACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATc                                                                                                                           <  1:606151/149‑1 (MQ=255)
                       aCGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCg                                                                                                                          >  1:740443/1‑149 (MQ=255)
                         gATTGGGTCCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATGGTGGATCGCGATTGGACGGAAGACGCTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGCTCGcc                                                                                                                        >  4:27772/1‑149 (MQ=255)
                             cGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgt                                                                                                                    >  1:667580/1‑149 (MQ=255)
                              ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtc                                                                                                                   <  2:790283/149‑1 (MQ=255)
                               gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                                  >  1:60567/1‑149 (MQ=255)
                               gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                                  >  3:581293/1‑149 (MQ=255)
                                 tCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCGCCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGtt                                                                                                                <  2:722006/149‑1 (MQ=255)
                                 tCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGtt                                                                                                                <  4:339151/149‑1 (MQ=255)
                                  cAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTc                                                                                                               <  2:325770/149‑1 (MQ=255)
                                   aaTACCCGCTGCGCAGTGTTTGTTAAAGCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTcc                                                                                                              >  2:384472/1‑149 (MQ=255)
                                   aaTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTcc                                                                                                              <  4:321186/149‑1 (MQ=255)
                                    aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                             >  3:379577/1‑149 (MQ=255)
                                     tACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                            <  2:712327/149‑1 (MQ=255)
                                        ccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATc                                                                                                         >  3:287781/1‑149 (MQ=255)
                                            tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                                                     <  4:7916/149‑1 (MQ=255)
                                                 aGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                                <  1:633677/149‑1 (MQ=255)
                                                  gtgtTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                                                               >  2:289588/1‑149 (MQ=255)
                                                     tttgttAAAACTTTCCATGCCATCTGCCCGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGAAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCCTGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                                                                            >  3:116670/1‑149 (MQ=255)
                                                     tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                                                                            <  1:201983/149‑1 (MQ=255)
                                                     tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                                                                            >  1:285430/1‑149 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTc                                                                                                               >  1:473014/1‑126 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                        >  4:209728/1‑149 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAAGTCGTTc                                                                                                               <  2:473014/126‑1 (MQ=255)
                                                          tAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                                                                       <  4:94568/149‑1 (MQ=255)
                                                             aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGtt                                                                                    >  2:865296/1‑149 (MQ=255)
                                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGAGTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                                  >  4:514589/1‑149 (MQ=255)
                                                                tttCCATGCCGTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                                                                 >  1:745289/1‑149 (MQ=255)
                                                                tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                                                                 <  2:641206/149‑1 (MQ=255)
                                                                  tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                               <  3:334155/149‑1 (MQ=255)
                                                                  tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                               <  3:150863/149‑1 (MQ=255)
                                                                     aTGCCATCTGCCAGCAAGCGCACCAATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                            >  1:745575/1‑149 (MQ=255)
                                                                      tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTca                                                                           <  1:76253/149‑1 (MQ=255)
                                                                         cATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                        >  2:797692/1‑133 (MQ=255)
                                                                         cATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                        <  1:797692/133‑1 (MQ=255)
                                                                          aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTCGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                                                                       >  2:468795/1‑149 (MQ=255)
                                                                           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtg                                                                      >  3:40726/1‑149 (MQ=255)
                                                                             tGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtgtg                                                                    >  4:847618/1‑149 (MQ=255)
                                                                               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                                  >  3:512512/1‑149 (MQ=255)
                                                                               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                                  >  2:310122/1‑149 (MQ=255)
                                                                               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                                  <  2:610855/149‑1 (MQ=255)
                                                                                 aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                                <  1:289588/149‑1 (MQ=255)
                                                                                 aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                                >  1:405196/1‑149 (MQ=255)
                                                                                 aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                                >  1:837571/1‑149 (MQ=255)
                                                                                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                         >  3:165024/1‑123 (MQ=255)
                                                                                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                         <  4:165024/123‑1 (MQ=255)
                                                                                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTggg                                                               <  4:777482/149‑1 (MQ=255)
                                                                                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTggg                                                               >  1:713059/1‑149 (MQ=255)
                                                                                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc                                                          <  4:40471/149‑1 (MQ=255)
                                                                                          aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACttt                                                       >  4:573445/1‑149 (MQ=255)
                                                                                             gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTccc                                                    <  3:657214/149‑1 (MQ=255)
                                                                                              aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                                                                      <  4:328477/114‑1 (MQ=255)
                                                                                              aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                                                                      >  3:328477/1‑114 (MQ=255)
                                                                                                              gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGc                                   >  4:129833/1‑149 (MQ=255)
                                                                                                               aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                                  >  1:760461/1‑149 (MQ=255)
                                                                                                               aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                                  >  4:382524/1‑149 (MQ=255)
                                                                                                                tCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGccc                                 <  4:845233/149‑1 (MQ=255)
                                                                                                                  aCCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCgg                               <  2:456606/149‑1 (MQ=255)
                                                                                                                     aTTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                            >  2:678056/1‑149 (MQ=255)
                                                                                                                        ggACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCAtt                         >  3:849959/1‑149 (MQ=255)
                                                                                                                         gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGCTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCAttt                        >  2:849969/1‑149 (MQ=255)
                                                                                                                          aCGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCAtttt                       >  1:471635/1‑149 (MQ=255)
                                                                                                                           cGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                      <  2:776048/149‑1 (MQ=255)
                                                                                                                             gAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                    <  1:384472/149‑1 (MQ=255)
                                                                                                                                gACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGa                 <  2:149157/149‑1 (MQ=255)
                                                                                                                                 aCGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAt                <  2:713059/149‑1 (MQ=255)
                                                                                                                                 aCGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAt                <  2:285430/149‑1 (MQ=255)
                                                                                                                                     tcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTgag            <  2:740443/149‑1 (MQ=255)
                                                                                                                                       agttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtgt                                                                     <  3:594054/90‑1 (MQ=255)
                                                                                                                                       agttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtgt                                                                     >  4:594054/1‑90 (MQ=255)
                                                                                                                                        gttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAtt         >  2:822523/1‑149 (MQ=255)
                                                                                                                                         ttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAttt        >  4:649209/1‑149 (MQ=255)
                                                                                                                                         ttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAttt        >  3:565337/1‑149 (MQ=255)
                                                                                                                                         ttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt          <  4:360651/147‑1 (MQ=255)
                                                                                                                                         ttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt          >  3:360651/1‑147 (MQ=255)
                                                                                                                                         ttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCAATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAttt        >  4:803656/1‑149 (MQ=255)
                                                                                                                                          tcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTc       >  3:543254/1‑149 (MQ=255)
                                                                                                                                          tcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATACCTCACACTGTGTTGGGTTCACTTTCCCCGGCATAATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGCTTTc       >  2:692519/1‑149 (MQ=255)
                                                                                                                                            aaaGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCAc     <  4:550384/149‑1 (MQ=255)
                                                                                                                                            aaaGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCAc     <  1:678056/149‑1 (MQ=255)
                                                                                                                                            aaaGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCAc     >  2:856735/1‑149 (MQ=255)
                                                                                                                                               gTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCaa  >  3:95806/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAA  >  NC_000913/1685381‑1685672

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