Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,969 A→C 5.0% Q15P (CAA→CCA)  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,9690AC 5.0% 784.3 / 22.8 218Q15P (CAA→CCA) yafEputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (107/100);  new base C (2/9);  total (109/109)
Fisher's exact test for biased strand distribution p-value = 5.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 227 total aligned reads displayed.

TACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTG  >  NC_000913/231821‑232118
                                                                                                                                                    |                                                                                                                                                     
 aCCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaaa                                                                                                                                                     >  2:2528348/1‑150 (MQ=255)
  ccGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaaaa                                                                                                                                                    >  2:735949/1‑150 (MQ=255)
  ccGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaaaa                                                                                                                                                    >  1:2908979/1‑150 (MQ=255)
     tGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGt                                                                                                                                                 >  1:2415773/1‑150 (MQ=255)
     tGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGACAGACCAACGTATGGCGCTGCGCAAAACGt                                                                                                                                                 >  2:1013594/1‑150 (MQ=255)
      ggTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTg                                                                                                                                                <  2:937176/150‑1 (MQ=255)
      ggTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGAACACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACCGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTg                                                                                                                                                <  1:1274677/150‑1 (MQ=255)
       gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                                                                                                                               <  2:1100812/150‑3 (MQ=255)
       gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGa                                                                                                                                               <  1:2744604/150‑1 (MQ=255)
        tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAg                                                                                                                                              >  1:52999/1‑150 (MQ=255)
        tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAg                                                                                                                                              >  1:1180468/1‑150 (MQ=255)
             aCGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCgg                                                                                                                                         <  1:1579299/150‑1 (MQ=255)
                tCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                                                                                                                      >  2:793066/1‑150 (MQ=255)
                tCAGTGAAGCTGCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATACATAAGGTATGTCAGGTCTGCCACAGGGCAGACCACCGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                                                                                                                      >  2:152701/1‑150 (MQ=255)
                  aGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGgtggtg                                                                                                                                    <  2:1076982/150‑1 (MQ=255)
                    tGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGc                                                                                                                                  >  1:1610271/1‑150 (MQ=255)
                    tGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGc                                                                                                                                  >  1:1812241/1‑150 (MQ=255)
                      aaGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGt                                                                                                                                >  2:2282852/1‑150 (MQ=255)
                       aGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTa                                                                                                                               <  1:2371763/150‑1 (MQ=255)
                          ttCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGa                                                                                                                            >  1:2705314/1‑150 (MQ=255)
                            cTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACt                                                                                                                          >  2:2166981/1‑150 (MQ=255)
                              gTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTa                                                                                                                        <  1:2486598/150‑1 (MQ=255)
                               tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                                                                                                                       >  2:2827417/1‑150 (MQ=255)
                               tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                                                                                                                       <  1:921030/150‑1 (MQ=255)
                               tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                                                                                                                       >  1:3078520/1‑150 (MQ=255)
                                aCTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATc                                                                                                                      >  2:543894/1‑150 (MQ=255)
                                  tGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACTAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTg                                                                                                                    >  2:2025957/1‑150 (MQ=255)
                                  tGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCGGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTg                                                                                                                    >  2:3013852/1‑150 (MQ=255)
                                              cTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTgg                                                                                                        >  1:1016444/1‑150 (MQ=255)
                                               ttCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGa                                                                                                       <  2:2763706/150‑1 (MQ=255)
                                               ttCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGa                                                                                                       <  2:2654219/150‑1 (MQ=255)
                                               ttCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGa                                                                                                       >  1:2102710/1‑150 (MQ=255)
                                                 ccGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATg                                                                                                     >  2:395088/1‑150 (MQ=255)
                                                   gATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTc                                                                                                   >  1:1419733/1‑150 (MQ=255)
                                                   gATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTc                                                                                                   >  1:3060329/1‑150 (MQ=255)
                                                   gATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTc                                                                                                   >  1:1398646/1‑150 (MQ=255)
                                                   gATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGGTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACCTGCTGGATGTc                                                                                                   >  2:2254431/1‑150 (MQ=255)
                                                    aTCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCg                                                                                                  <  1:2028758/150‑1 (MQ=255)
                                                         aagccGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGca                                                                                             <  2:1906494/147‑1 (MQ=255)
                                                              gCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGc                                                                                        >  1:1162867/1‑150 (MQ=255)
                                                                  cTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGcc                                                                                    >  2:1614851/1‑150 (MQ=255)
                                                                     gTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGaa                                                                                 <  2:873715/150‑1 (MQ=255)
                                                                      ttGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAg                                                                                >  2:3100431/1‑150 (MQ=255)
                                                                          aTTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCg                                                                            <  1:2157577/150‑1 (MQ=255)
                                                                                tCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACt                                                                      >  1:2877496/1‑150 (MQ=255)
                                                                                     gCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaa                                                                 >  1:2638333/1‑150 (MQ=255)
                                                                                     gCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaa                                                                 >  1:2638289/1‑150 (MQ=255)
                                                                                      cAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaaa                                                                <  2:1043547/150‑1 (MQ=255)
                                                                                         gCCTGATAAATAAGGTATGTCAGGTCTGCCACTGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCTTATGACTTATCTGACCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTTAAAAatat                                                             <  2:2873791/150‑1 (MQ=255)
                                                                                         gCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAatat                                                             <  2:1787862/150‑1 (MQ=255)
                                                                                            tGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcac                                                          >  1:970714/1‑150 (MQ=255)
                                                                                            tGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcac                                                          >  1:662715/1‑150 (MQ=255)
                                                                                             gATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcacc                                                         >  1:1242001/1‑150 (MQ=255)
                                                                                             gATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcacc                                                         >  1:381278/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)
                                                                                                 aataaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCAccc                                                     >  2:2718336/1‑150 (MQ=255)
                                                                                                   taaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGc                                                   <  1:2177480/150‑1 (MQ=255)
                                                                                                   taaGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGc                                                   <  2:572640/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)
                                                                                                           gTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCCCCACCCGCCAGGGata                                           >  2:2084817/1‑150 (MQ=255)
                                                                                                           gTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGata                                           <  2:668244/150‑1 (MQ=255)
                                                                                                            tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCACCTGAAAAATATCACCACCCGCCAGGGatat                                          >  2:95954/1‑150 (MQ=255)
                                                                                                            tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGatat                                          <  1:2104097/150‑1 (MQ=255)
                                                                                                            tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGatat                                          >  2:2585116/1‑150 (MQ=255)
                                                                                                            tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTGATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGatat                                          <  2:2738797/150‑1 (MQ=255)
                                                                                                             cAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATg                                         <  1:2282852/150‑1 (MQ=255)
                                                                                                              aGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGc                                        <  1:2528348/150‑1 (MQ=255)
                                                                                                              aGGTCTGCCACAGGCCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGc                                        <  2:2882709/150‑1 (MQ=255)
                                                                                                               ggTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGcc                                       <  1:1780038/150‑1 (MQ=255)
                                                                                                                 tCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGa                                     >  1:1643941/1‑150 (MQ=255)
                                                                                                                  cTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGaa                                    >  1:426600/1‑150 (MQ=255)
                                                                                                                     ccACAGGGCACACCAACGGTTGGCGCAGCGCAAAACGTGATCGAGCTGGTGGCGTATGACTTATCTGCCCACATCCTGGAAGACGTGTCACAAGCTTCCGAAGCACGGCAACTCAAAAATATCACCACCCGCCAGGGATATGCCGAAAGt                                 >  2:1848864/1‑150 (MQ=255)
                                                                                                                           ggCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCa                           >  1:2560262/1‑150 (MQ=255)
                                                                                                                           ggCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCa                           >  1:318592/1‑150 (MQ=255)
                                                                                                                            gCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCAt                          <  2:1704925/150‑1 (MQ=255)
                                                                                                                            gCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCAt                          <  1:745852/150‑1 (MQ=255)
                                                                                                                                aCCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGc                      >  2:1885664/1‑150 (MQ=255)
                                                                                                                                aCCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGc                      <  1:2030226/150‑1 (MQ=255)
                                                                                                                                 ccAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGcc                     <  1:1663100/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)
                                                                                                                                   aaCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGa                   >  2:1349391/1‑150 (MQ=255)
                                                                                                                                    aCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGAt                  >  1:1281764/1‑150 (MQ=255)
                                                                                                                                      gTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATaa                >  1:2452847/1‑150 (MQ=255)
                                                                                                                                      gTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATaa                >  1:1879884/1‑150 (MQ=255)
                                                                                                                                      gTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATaa                >  1:2828889/1‑150 (MQ=255)
                                                                                                                                          ggCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCa            >  1:635161/1‑150 (MQ=255)
                                                                                                                                              cTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTg        <  2:140301/150‑1 (MQ=255)
                                                                                                                                              cTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTg        <  2:922876/150‑1 (MQ=255)
                                                                                                                                               tGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGa       >  1:380984/1‑150 (MQ=255)
                                                                                                                                               tGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGa       <  1:1761640/150‑1 (MQ=255)
                                                                                                                                                   cAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATAtt   >  1:787408/1‑150 (MQ=255)
                                                                                                                                                   cAAAACGTGAGCGCGGTGGTGGCGTATGACGTATCTGCCCACATGCTGGATGTCGTGGCACACGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATAtt   <  2:4103/150‑1 (MQ=255)
                                                                                                                                                    aaaaCGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTg  >  1:2944592/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTG  >  NC_000913/231821‑232118

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