Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,959 G→T 5.2% G12C (GGC→TGC)  yafE → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,9590GT 5.2% 775.3 / 22.0 212G12C (GGC→TGC) yafEputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (108/93);  new base T (2/9);  total (110/102)
Fisher's exact test for biased strand distribution p-value = 2.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 217 total aligned reads displayed.

TTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCA  >  NC_000913/231811‑232108
                                                                                                                                                    |                                                                                                                                                     
ttttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTgg                                                                                                                                                      >  1:88751/1‑150 (MQ=255)
ttttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTgg                                                                                                                                                      >  1:3044717/1‑150 (MQ=255)
 tttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGgc                                                                                                                                                     >  2:1550087/1‑150 (MQ=255)
  ttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGgcg                                                                                                                                                    <  2:142416/150‑1 (MQ=255)
  ttGTTTAATACCGTGGTCTGAACGTCAGTTAAGCTTCAGTAATGGATACGCTCGATTCATATCACAATCCGATAATCGTTGAATTCAGAACCGGCAAGCCTGAGAAATAAGGTATGTCAGGTCTGACACAGGGCAGACCAACGTTTGtcg                                                                                                                                                    >  2:3018564/1‑150 (MQ=255)
       ttCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTgcg                                                                                                                                               <  2:2985611/150‑1 (MQ=255)
        tCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTgcgc                                                                                                                                              >  2:1214887/1‑150 (MQ=255)
          tACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaa                                                                                                                                            >  2:1084287/1‑150 (MQ=255)
           aCCGTGGTCTGAACGTCTGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaaa                                                                                                                                           <  1:141578/150‑1 (MQ=255)
            ccGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaaaa                                                                                                                                          >  1:2908979/1‑150 (MQ=255)
             cGTGGTCTGAACGTCTGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAAc                                                                                                                                         <  1:2510591/150‑1 (MQ=255)
               tGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGt                                                                                                                                       >  1:2415773/1‑150 (MQ=255)
                ggTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTg                                                                                                                                      <  2:937176/150‑1 (MQ=255)
                ggTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGAACACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACCGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTg                                                                                                                                      <  1:1274677/150‑1 (MQ=255)
                 gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                                                                                                                     <  2:3050682/150‑3 (MQ=255)
                 gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                                                                                                                     <  2:2605760/150‑3 (MQ=255)
                 gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                                                                                                                     <  2:2512184/150‑3 (MQ=255)
                 gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                                                                                                                     <  2:1100812/150‑3 (MQ=255)
                 gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGa                                                                                                                                     <  1:2744554/150‑1 (MQ=255)
                  tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAg                                                                                                                                    >  1:1180468/1‑150 (MQ=255)
                  tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAg                                                                                                                                    >  1:52999/1‑150 (MQ=255)
                       aCGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCgg                                                                                                                               <  1:1579299/150‑1 (MQ=255)
                         gTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtg                                                                                                                             <  2:2253191/150‑1 (MQ=255)
                              tGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGc                                                                                                                        <  1:1845332/150‑1 (MQ=255)
                              tGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGc                                                                                                                        >  1:1812241/1‑150 (MQ=255)
                              tGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGc                                                                                                                        >  1:28607/1‑150 (MQ=255)
                               gAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCg                                                                                                                       <  1:349879/150‑1 (MQ=255)
                                aaGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGt                                                                                                                      >  2:2282852/1‑150 (MQ=255)
                                 aGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTa                                                                                                                     <  1:2371763/150‑1 (MQ=255)
                                    ttCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTGTGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGa                                                                                                                  >  1:2608081/1‑150 (MQ=255)
                                      cTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACt                                                                                                                >  2:2166981/1‑150 (MQ=255)
                                        ggaCTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTa                                                                                                              <  1:2426861/148‑1 (MQ=255)
                                        gTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTa                                                                                                              >  1:2795089/1‑150 (MQ=255)
                                         tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                                                                                                             >  1:3078609/1‑150 (MQ=255)
                                         tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                                                                                                             <  1:921083/150‑1 (MQ=255)
                                          aCTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATc                                                                                                            >  2:543894/1‑150 (MQ=255)
                                           cTGGTTACACGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCt                                                                                                           >  2:774921/1‑150 (MQ=255)
                                            tGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTg                                                                                                          >  2:2781095/1‑150 (MQ=255)
                                             ggTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGc                                                                                                         <  1:2455696/150‑1 (MQ=255)
                                              gTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGcc                                                                                                        <  2:1415122/150‑1 (MQ=255)
                                                     gcgcGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGc                                                                                                 <  1:2535927/150‑1 (MQ=255)
                                                        cTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTgg                                                                                              >  1:1016444/1‑150 (MQ=255)
                                                         ttCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGa                                                                                             >  1:2102710/1‑150 (MQ=255)
                                                          tCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGAt                                                                                            <  2:203449/150‑1 (MQ=255)
                                                          tCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGAt                                                                                            <  2:814612/150‑1 (MQ=255)
                                                           ccGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATg                                                                                           >  2:395088/1‑150 (MQ=255)
                                                             gATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTc                                                                                         >  1:3060565/1‑150 (MQ=255)
                                                              aTCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCg                                                                                        >  1:439001/1‑150 (MQ=255)
                                                              aTCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGGCTGACTGATCTGCCCACATGCTGGATGTCg                                                                                        >  1:2670824/1‑150 (MQ=255)
                                                                cacaATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGATCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTg                                                                                      <  1:1650221/150‑1 (MQ=255)
                                                                cacaATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTg                                                                                      <  1:2002180/150‑1 (MQ=255)
                                                                        gCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGc                                                                              >  1:1162867/1‑150 (MQ=255)
                                                                        gCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGc                                                                              >  1:1254993/1‑150 (MQ=255)
                                                                            cTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGcc                                                                          >  2:1614851/1‑150 (MQ=255)
                                                                               gTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGaa                                                                       <  2:873722/150‑1 (MQ=255)
                                                                                  gAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGcc                                                                    >  2:254422/1‑150 (MQ=255)
                                                                                  gAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGcc                                                                    >  2:1753901/1‑150 (MQ=255)
                                                                                  gAATTCAGTCCCGGCAAGCATGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACAAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCACACATGCTGGATGTCGTGGCACAAGCTGCCGAAGcc                                                                    >  2:607230/1‑150 (MQ=255)
                                                                                    aTTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCg                                                                  <  1:2157577/150‑1 (MQ=255)
                                                                                      tCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGc                                                                >  2:189867/1‑150 (MQ=255)
                                                                                          tCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACt                                                            <  2:2735107/150‑1 (MQ=255)
                                                                                               gCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaa                                                       >  1:2638326/1‑150 (MQ=255)
                                                                                                cAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaaa                                                      <  2:1043547/150‑1 (MQ=255)
                                                                                                   gCCTGATAAATAAGGTATGTCAGGTCTGCCACTGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCTTATGACTTATCTGACCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTTAAAAatat                                                   <  2:2873791/150‑1 (MQ=255)
                                                                                                      tGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcac                                                <  2:2687322/150‑1 (MQ=255)
                                                                                                      tGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcac                                                >  1:970714/1‑150 (MQ=255)
                                                                                                         taaataaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcaccac                                             >  2:1691872/1‑150 (MQ=255)
                                                                                                           aataaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCAccc                                           >  1:1947000/1‑150 (MQ=255)
                                                                                                           aataaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCAccc                                           >  2:2392282/1‑150 (MQ=255)
                                                                                                             taaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGc                                         <  2:572640/150‑1 (MQ=255)
                                                                                                             taaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGc                                         <  1:2177480/150‑1 (MQ=255)
                                                                                                             taaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGc                                         <  1:5697/150‑1 (MQ=255)
                                                                                                                ggTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAg                                      <  2:2000001/150‑1 (MQ=255)
                                                                                                                 gTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAgg                                     >  1:2883446/1‑150 (MQ=255)
                                                                                                                  tATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAggg                                    <  1:730241/150‑1 (MQ=255)
                                                                                                                     gTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGata                                 <  2:668244/150‑1 (MQ=255)
                                                                                                                      tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCACCTGAAAAATATCACCACCCGCCAGGGatat                                >  2:95954/1‑150 (MQ=255)
                                                                                                                      tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGatat                                <  1:2104097/150‑1 (MQ=255)
                                                                                                                        aGGTCTGCCACAGGCCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGc                              <  2:2882709/150‑1 (MQ=255)
                                                                                                                         ggTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGcc                             <  1:1780038/150‑1 (MQ=255)
                                                                                                                          gacTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCg                            >  1:149506/3‑150 (MQ=255)
                                                                                                                          gTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCg                            >  1:1626455/1‑150 (MQ=255)
                                                                                                                           tCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGa                           >  1:2914613/1‑150 (MQ=255)
                                                                                                                               ccACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGt                       >  1:2279353/1‑150 (MQ=255)
                                                                                                                               ccACAGGGCACACCAACGGTTGGCGCAGCGCAAAACGTGATCGAGCTGGTGGCGTATGACTTATCTGCCCACATCCTGGAAGACGTGTCACAAGCTTCCGAAGCACGGCAACTCAAAAATATCACCACCCGCCAGGGATATGCCGAAAGt                       >  2:1848864/1‑150 (MQ=255)
                                                                                                                                    gggCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGcc                  <  1:395088/150‑1 (MQ=255)
                                                                                                                                      gCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCAt                <  1:745852/150‑1 (MQ=255)
                                                                                                                                      gCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCAt                <  2:1704925/150‑1 (MQ=255)
                                                                                                                                          aCCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGc            <  1:2030226/150‑1 (MQ=255)
                                                                                                                                          aCCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGc            >  2:2817505/1‑150 (MQ=255)
                                                                                                                                          aCCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCAGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGc            >  2:2763026/1‑150 (MQ=255)
                                                                                                                                           ccAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGcc           <  1:1663102/150‑1 (MQ=255)
                                                                                                                                            tcaCGTTTTGCGCAGCGCCAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCg          >  2:2414253/3‑150 (MQ=255)
                                                                                                                                            cAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCGTTTGCCg          >  2:2785489/1‑150 (MQ=255)
                                                                                                                                            cAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCg          >  2:583613/1‑150 (MQ=255)
                                                                                                                                             aaCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGa         <  2:421821/150‑1 (MQ=255)
                                                                                                                                              aCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGAt        >  1:1281764/1‑150 (MQ=255)
                                                                                                                                                gTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATaa      >  1:1879884/1‑150 (MQ=255)
                                                                                                                                                    ggCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGACAGTCTGCCATTTGCCGATAACGCa  >  2:1907903/1‑150 (MQ=255)
                                                                                                                                                    ggCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCa  >  1:635161/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCA  >  NC_000913/231811‑232108

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