Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,225,322 G→T 5.4% E499* (GAG→TAG) ‡ metH → cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,3220GT 5.4% 336.6 / 5.8 99E499* (GAG→TAG) ‡metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (51/42);  new base T (2/4);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 102 total aligned reads displayed.

GCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGAC  >  NC_000913/4225176‑4225471
                                                                                                                                                  |                                                                                                                                                     
gCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgaga                                                                                                                                                    <  2:86991/150‑1 (MQ=255)
gCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgaga                                                                                                                                                    <  2:340264/150‑1 (MQ=255)
 cAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAt                                                                                                                                                   >  2:406931/1‑150 (MQ=255)
 cAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAt                                                                                                                                                   >  1:1458064/1‑150 (MQ=255)
   aaGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAttt                                                                                                                                                 >  2:1410799/1‑150 (MQ=255)
     ggCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGc                                                                                                                                               <  2:936373/150‑1 (MQ=255)
     ggCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGc                                                                                                                                               >  1:1468511/1‑150 (MQ=255)
         ttGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtg                                                                                                                                           <  1:1461819/150‑2 (MQ=255)
         ttGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTc                                                                                                                                           >  1:677123/1‑150 (MQ=255)
          tGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtg                                                                                                                                           >  2:863895/1‑148 (MQ=255)
          tGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtg                                                                                                                                           <  1:863895/149‑2 (MQ=255)
             tAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                                                                       >  1:998743/1‑150 (MQ=255)
              aaCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAACTTGTTGCGTCGCTAAGGTGCGGCAGTGGTGGGAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCg                                                                                                                                      <  2:1403988/150‑1 (MQ=255)
              aaCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCg                                                                                                                                      >  1:737580/1‑150 (MQ=255)
              aaCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCg                                                                                                                                      <  1:958387/150‑1 (MQ=255)
               aCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGCAAATCTCGATTTTCCGtgcgcgag                                                                                                                                     <  1:1525460/150‑8 (MQ=255)
                 acaaTCTAGATGAAAAAGAGAGTCCAATCCTTTATCAATCACTCGAAATAGTAGCATAACTACAATACGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTAc                                                                                                                                   >  2:618916/4‑150 (MQ=255)
                    aTCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTTGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACaaa                                                                                                                                <  1:1288362/150‑1 (MQ=255)
                    aTCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACaaa                                                                                                                                <  1:705519/150‑1 (MQ=255)
                             aaaGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCAcc                                                                                                                       <  1:952728/150‑1 (MQ=255)
                             aaaGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCAcc                                                                                                                       >  2:1055497/1‑150 (MQ=255)
                                     cGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGt                                                                                                               >  1:1522024/1‑150 (MQ=255)
                                      gTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGtt                                                                                                              <  2:1458064/150‑1 (MQ=255)
                                       tCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTg                                                                                                             >  1:453067/1‑150 (MQ=255)
                                         gATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGc                                                                                                           <  1:766703/150‑1 (MQ=255)
                                          aTGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCt                                                                                                          <  1:192178/150‑1 (MQ=255)
                                            gCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTc                                                                                                        <  2:1322798/150‑1 (MQ=255)
                                             ccTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                                       >  1:454783/1‑150 (MQ=255)
                                             ccTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                                       >  2:510877/1‑150 (MQ=255)
                                               tttATCCATCACGCGAAATTGTTGCGTCGCTACGTTGCGGCAGTGTTGGTAATGGCCTTTGACGAAGAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCg                                                                                                     <  1:1291147/150‑1 (MQ=255)
                                               tttATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCg                                                                                                     <  2:525782/150‑1 (MQ=255)
                                                 tATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGcc                                                                                                   <  1:406931/150‑1 (MQ=255)
                                                     cATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaa                                                                                               >  2:1120471/1‑150 (MQ=255)
                                                     cATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaa                                                                                               >  2:1150048/1‑150 (MQ=255)
                                                     cATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaa                                                                                               >  2:658902/1‑150 (MQ=255)
                                                      aTCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaag                                                                                              <  2:327008/150‑1 (MQ=255)
                                                      aTCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaag                                                                                              <  1:736498/150‑1 (MQ=255)
                                                      aTCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaag                                                                                              <  2:489515/150‑1 (MQ=255)
                                                      aTCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaag                                                                                              <  2:273957/150‑1 (MQ=255)
                                                       tCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAACGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaaga                                                                                             <  2:3332/150‑1 (MQ=255)
                                                        cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGat                                                                                            >  1:1264791/1‑150 (MQ=255)
                                                         aCGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGata                                                                                           >  2:39223/1‑150 (MQ=255)
                                                          cgcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGatat                                                                                          >  1:1465695/1‑150 (MQ=255)
                                                           gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                                         <  2:1522024/150‑1 (MQ=255)
                                                           gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                                         <  1:360632/150‑1 (MQ=255)
                                                            cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatca                                                                                        <  2:881151/150‑1 (MQ=255)
                                                            cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatca                                                                                        >  1:301269/1‑150 (MQ=255)
                                                            cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatca                                                                                        >  1:319654/1‑150 (MQ=255)
                                                            cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatca                                                                                        >  1:319656/1‑150 (MQ=255)
                                                            cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatca                                                                                        >  2:1300109/1‑150 (MQ=255)
                                                               aaTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTTGCTTCCCGCCAGAAGATATCATCt                                                                                     >  1:796190/1‑150 (MQ=255)
                                                                    ttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCTATATTTTTATTTGTGAGTACTAAATCCTAACTGAAGAGGTTTTATTACTTTTAGAATATATAATCTTCGAc                                                                                <  2:1453993/150‑1 (MQ=255)
                                                                    ttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGAc                                                                                <  1:1267363/150‑1 (MQ=255)
                                                                    ttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGAc                                                                                >  2:1136173/1‑150 (MQ=255)
                                                                     tgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGAcc                                                                               >  1:791285/1‑150 (MQ=255)
                                                                     tgCGTCGCTACGGTGCGGCAGTGGGGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGAcc                                                                               <  1:533344/150‑1 (MQ=255)
                                                                         tCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCaaa                                                                           <  1:1136173/150‑1 (MQ=255)
                                                                         tCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCaaa                                                                           <  1:1336889/150‑1 (MQ=255)
                                                                         tCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCaaa                                                                           <  1:220272/150‑1 (MQ=255)
                                                                         tCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCaaa                                                                           <  1:1410799/150‑1 (MQ=255)
                                                                          cGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAAc                                                                          >  1:898260/1‑150 (MQ=255)
                                                                                      gCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTc                                                              <  1:618916/150‑1 (MQ=255)
                                                                                       cAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCg                                                             <  1:604698/150‑1 (MQ=255)
                                                                                          tggtggTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCaa                                                          >  2:1032560/1‑150 (MQ=255)
                                                                                           ggtggtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAAc                                                         >  2:676002/1‑150 (MQ=255)
                                                                                            gtggtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACt                                                        <  2:1468511/150‑1 (MQ=255)
                                                                                              ggtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTgg                                                      >  1:1354861/1‑150 (MQ=255)
                                                                                              ggtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTgg                                                      <  2:737580/150‑1 (MQ=255)
                                                                                               gtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGc                                                     <  2:1136687/150‑1 (MQ=255)
                                                                                               gtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGc                                                     <  2:1465695/150‑1 (MQ=255)
                                                                                                 aaTGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAt                                                   >  1:951500/1‑150 (MQ=255)
                                                                                                 aaTGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAt                                                   >  2:135845/1‑150 (MQ=255)
                                                                                                      ccTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAag                                              >  1:1492196/1‑150 (MQ=255)
                                                                                                        tttGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAagag                                            <  1:658902/150‑1 (MQ=255)
                                                                                                           taCGATCAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGcac                                         <  2:1499899/149‑1 (MQ=255)
                                                                                                           gACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGcac                                         <  1:1300109/150‑1 (MQ=255)
                                                                                                            aCGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGcaca                                        >  2:1108773/1‑150 (MQ=255)
                                                                                                                   gggACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTa                                 >  1:434810/1‑150 (MQ=255)
                                                                                                                       cAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTAcgcg                             >  2:772659/1‑150 (MQ=255)
                                                                                                                       cAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTAcgcg                             >  1:1136782/1‑150 (MQ=255)
                                                                                                                         ggCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCa                           >  2:549456/1‑150 (MQ=255)
                                                                                                                         ggCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCa                           >  2:1178657/1‑150 (MQ=255)
                                                                                                                         ggCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCa                           >  1:580236/1‑150 (MQ=255)
                                                                                                                           ccGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAgg                         >  2:109416/1‑150 (MQ=255)
                                                                                                                           ccGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAgg                         >  1:690937/1‑150 (MQ=255)
                                                                                                                             gATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGAc                       <  1:1055497/150‑1 (MQ=255)
                                                                                                                               tACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACtt                     >  1:615274/1‑150 (MQ=255)
                                                                                                                                 cTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTa                   >  1:706209/1‑150 (MQ=255)
                                                                                                                                 cTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTa                   >  2:1391208/1‑150 (MQ=255)
                                                                                                                                  cgcccgaCGGCAAATCTCGATTTTCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTAt                  >  2:533598/7‑150 (MQ=255)
                                                                                                                                   cgcgCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATc                 >  1:859806/1‑150 (MQ=255)
                                                                                                                                         cGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGcgcg           <  2:1492196/150‑1 (MQ=255)
                                                                                                                                           gAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCg                >  1:766810/1‑143 (MQ=255)
                                                                                                                                           gAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCg                <  2:766810/143‑1 (MQ=255)
                                                                                                                                           gAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCgtg         <  2:453067/150‑1 (MQ=255)
                                                                                                                                           gAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCgtg         <  2:301269/150‑1 (MQ=255)
                                                                                                                                           gAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCgtg         >  1:845964/1‑150 (MQ=255)
                                                                                                                                           gAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCgtg         >  1:1160344/1‑150 (MQ=255)
                                                                                                                                               aTCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGaa     >  1:813034/1‑150 (MQ=255)
                                                                                                                                                  gagaTTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGAc  >  2:65029/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGAC  >  NC_000913/4225176‑4225471

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: