Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,225,323 A→C 6.0% E499A (GAG→GCG) ‡ metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,3230AC 6.0% 533.3 / 21.9 149E499A (GAG→GCG) ‡metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (67/73);  new base C (4/5);  total (71/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 161 total aligned reads displayed.

GGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGC  >  NC_000913/4225175‑4225461
                                                                                                                                                    |                                                                                                                                          
ggCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCga                                                                                                                                            <  4:848180/149‑1 (MQ=255)
   aaaGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAt                                                                                                                                         >  4:710349/1‑149 (MQ=255)
   aaaGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAt                                                                                                                                         >  3:208690/1‑149 (MQ=255)
        cATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGcc                                                                                                                                    >  3:811235/1‑149 (MQ=255)
           tGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTc                                                                                                                                 <  1:440970/149‑1 (MQ=255)
            gTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCg                                                                                                                                >  3:1468128/1‑149 (MQ=255)
            gTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCAATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCgg                                                                                                                               <  2:1441270/149‑1 (MQ=255)
               aaCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                                                             <  1:692061/149‑1 (MQ=255)
               aaCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                                                             <  1:996191/149‑1 (MQ=255)
               aaCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                                                             <  4:208690/149‑1 (MQ=255)
                aCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCg                                                                                                                            >  2:630075/1‑149 (MQ=255)
                  tctATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGc                                                                                                                                     >  2:88847/1‑138 (MQ=255)
                  tctATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGc                                                                                                                                     <  1:88847/138‑1 (MQ=255)
                  cttATCTTGATTAAAGGGGGCGGCCATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTa                                                                                                                          <  2:159875/147‑1 (MQ=255)
                       ctcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACaaaa                                                                                                                     <  3:968839/149‑1 (MQ=255)
                        tcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtgcg                                                                                                                               >  2:474522/1‑135 (MQ=255)
                        tcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtgcg                                                                                                                               >  2:253507/1‑135 (MQ=255)
                        tcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtgcg                                                                                                                               <  1:474522/138‑4 (MQ=255)
                        tcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtgcg                                                                                                                               <  1:253507/138‑4 (MQ=255)
                         cGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACaaaa                                                                                                                     >  3:280807/1‑147 (MQ=255)
                         cGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACaaaa                                                                                                                     <  4:280807/147‑1 (MQ=255)
                         cGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATc                                                                                                                   <  1:630075/149‑1 (MQ=255)
                          gATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATcc                                                                                                                  >  1:1307974/1‑149 (MQ=255)
                            tGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTc                                                                                                                >  1:699299/1‑149 (MQ=255)
                            tGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTc                                                                                                                >  3:707798/1‑149 (MQ=255)
                             gAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCa                                                                                                               >  4:248698/1‑149 (MQ=255)
                               aaGAGGGGGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACCAAATCCTCAcc                                                                                                             >  2:253934/1‑149 (MQ=255)
                               aaGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgag                                                                                                                                           <  4:585483/119‑1 (MQ=255)
                                 gagGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGt                                                                                                                                  <  2:1357821/126‑1 (MQ=255)
                                  agGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGaa                                                                                                          <  3:409371/149‑1 (MQ=255)
                                        tCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCa                                                                                                               >  2:816105/1‑138 (MQ=255)
                                          gATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGt                                                                                                                                  >  4:612560/1‑117 (MQ=255)
                                           aTGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGc                                                                                                 >  4:424477/1‑149 (MQ=255)
                                            tGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCt                                                                                                >  4:1077007/1‑149 (MQ=255)
                                             gCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGTACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCa                                                                                                               >  4:428332/1‑133 (MQ=255)
                                             gCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCa                                                                                                               <  3:428332/133‑1 (MQ=255)
                                               cTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                             >  2:1271709/1‑149 (MQ=255)
                                               cTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                             >  1:834420/1‑149 (MQ=255)
                                               cTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                             >  4:839796/1‑149 (MQ=255)
                                               cTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                             <  4:1362936/149‑1 (MQ=255)
                                               cTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                             >  1:1434317/1‑149 (MQ=255)
                                                  tATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCga                                                                                                                                            >  3:518601/1‑99 (MQ=255)
                                                  tATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGc                                                                                          <  3:1001699/149‑1 (MQ=255)
                                                   aTCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGcc                                                                                         <  3:1077007/149‑1 (MQ=255)
                                                       aTCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaa                                                                                     <  4:707798/149‑1 (MQ=255)
                                                        tCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaag                                                                                    >  2:1034601/1‑149 (MQ=255)
                                                         cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgag                                                                                                                                           >  1:902228/1‑93 (MQ=255)
                                                         cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgag                                                                                                                                           <  2:902228/93‑1 (MQ=255)
                                                          aCGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGtt                                                                                                    <  1:1029291/131‑1 (MQ=255)
                                                          aCGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGtt                                                                                                    >  2:1029291/1‑131 (MQ=255)
                                                          aCGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGt                                                                                                     <  3:643283/130‑1 (MQ=255)
                                                          aCGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGat                                                                                  <  2:1434317/149‑1 (MQ=255)
                                                           cgcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgaga                                                                                                                                          >  2:153969/1‑92 (MQ=255)
                                                           cgcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgaga                                                                                                                                          <  1:153969/92‑1 (MQ=255)
                                                            gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGAGAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCggg                                                                                                                              <  2:230591/103‑1 (MQ=255)
                                                            gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGt                                                                                                     >  3:607750/1‑128 (MQ=255)
                                                            gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGCGGTAATGGCCTTTGACGAACAGGGAGAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCggg                                                                                                                              >  1:230591/1‑103 (MQ=255)
                                                             cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                               <  3:1420081/149‑1 (MQ=255)
                                                             cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                               <  3:248698/149‑1 (MQ=255)
                                                             cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                               >  3:477628/1‑149 (MQ=255)
                                                              gAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatca                                                                              >  1:1234066/1‑149 (MQ=255)
                                                               aaaTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCt                                                                                                                 <  4:1017300/113‑1 (MQ=255)
                                                                 aTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCt                                                                           <  4:497410/149‑1 (MQ=255)
                                                                    gttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCg                                                                        >  3:738581/1‑149 (MQ=255)
                                                                     ttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACaaaa                                                                                                                     <  2:955274/103‑1 (MQ=255)
                                                                     ttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACaaaa                                                                                                                     >  1:955274/1‑103 (MQ=255)
                                                                     ttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGa                                                                       <  1:1432047/149‑1 (MQ=255)
                                                                     ttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGa                                                                       <  4:942612/149‑1 (MQ=255)
                                                                         gTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCa                                                                                                               <  4:739719/105‑1 (MQ=255)
                                                                                  gTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgag                                                                                                                                           <  1:629276/68‑1 (MQ=255)
                                                                                  gTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgag                                                                                                                                           >  2:629276/1‑68 (MQ=255)
                                                                                    gCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTcgc                                                        <  2:834420/149‑1 (MQ=255)
                                                                                     cGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGa                                                                       <  1:1033037/133‑1 (MQ=255)
                                                                                        cAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACTTCTTCGCGGTc                                                    <  3:839796/149‑1 (MQ=255)
                                                                                        cAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTc                                                    <  1:155469/149‑1 (MQ=255)
                                                                                                 tAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGc                                           >  4:473398/1‑149 (MQ=255)
                                                                                                  aaTGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCa                                          <  3:466469/149‑1 (MQ=255)
                                                                                                   aTGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAt                                         <  3:704771/149‑1 (MQ=255)
                                                                                                   aTGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAt                                         <  3:1064038/149‑1 (MQ=255)
                                                                                                   aTGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAt                                         >  3:554748/1‑149 (MQ=255)
                                                                                                    tGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAtt                                        <  2:1307974/149‑1 (MQ=255)
                                                                                                     ggCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGGCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTg                                       >  4:67127/1‑149 (MQ=255)
                                                                                                       ccTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGaa                                     <  2:699299/149‑1 (MQ=255)
                                                                                                        cTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGtt                                                                                                    >  4:385806/1‑85 (MQ=255)
                                                                                                        cTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGtt                                                                                                    <  3:385806/85‑1 (MQ=255)
                                                                                                        cTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAAc                                                                <  4:671516/121‑1 (MQ=255)
                                                                                                         tttGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTg                                       <  3:1009613/145‑1 (MQ=255)
                                                                                                         tttGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAaga                                   <  4:1468204/149‑1 (MQ=255)
                                                                                                         tttGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAaga                                   <  1:597651/149‑1 (MQ=255)
                                                                                                         tttGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAaga                                   <  4:554748/149‑1 (MQ=255)
                                                                                                          ttGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGGCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAagag                                  <  3:67127/149‑1 (MQ=255)
                                                                                                             aCGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAaga                                                                                                        <  4:1291956/76‑1 (MQ=255)
                                                                                                             aCGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAaga                                                                                                        >  3:1291956/1‑76 (MQ=255)
                                                                                                               gAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGca                                <  2:1066809/146‑1 (MQ=255)
                                                                                                               gAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGca                                >  1:1066809/1‑146 (MQ=255)
                                                                                                                         aGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTAcgcg                   >  4:483136/1‑149 (MQ=255)
                                                                                                                            ccGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAg                >  4:364933/1‑149 (MQ=255)
                                                                                                                                    cgcgCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTAt        <  1:253934/149‑1 (MQ=255)
                                                                                                                                      cgcACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCg      <  3:990838/149‑1 (MQ=255)
                                                                                                                                           ggAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAgg               >  4:1239029/1‑135 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGC  >  NC_000913/4225175‑4225461

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