Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,964 T→A 8.0% A13A (GCT→GCA 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,9640TA 8.0% 262.4 / 10.1 75A13A (GCT→GCAyafEputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (44/25);  new base A (2/4);  total (46/29)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTT  >  NC_000913/231821‑232111
                                                                                                                                               |                                                                                                                                                   
tACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCa                                                                                                                                                >  3:678165/1‑149 (MQ=255)
tACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCa                                                                                                                                                >  3:101187/1‑149 (MQ=255)
  ccGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCcaa                                                                                                                                              >  4:477730/1‑149 (MQ=255)
      ggTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGt                                                                                                                                          <  2:666309/149‑1 (MQ=255)
        tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                                                                                                                        <  1:629978/149‑3 (MQ=255)
        tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                                                                                                                        <  3:477730/149‑3 (MQ=255)
        tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGa                                                                                                                                        >  4:48538/1‑149 (MQ=255)
        tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGa                                                                                                                                        <  4:158394/149‑1 (MQ=255)
            aaCGTCAGTGAAGCTTCTGTACTGGTTACACGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAgcgc                                                                                                                                    >  4:24205/1‑149 (MQ=255)
                tCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtg                                                                                                                                >  4:323497/1‑149 (MQ=255)
                 cAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                                                                                                               >  3:168244/1‑149 (MQ=255)
                      aaGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCg                                                                                                                          >  1:27634/1‑149 (MQ=255)
                        gCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTa                                                                                                                        >  1:590653/1‑149 (MQ=255)
                         cTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTAt                                                                                                                       <  4:168244/149‑1 (MQ=255)
                          ttCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATg                                                                                                                      >  1:76764/1‑149 (MQ=255)
                               tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGt                                                                                                                                          <  2:669201/124‑1 (MQ=255)
                               tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGt                                                                                                                                          >  1:669201/1‑124 (MQ=255)
                                  tGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCt                                                                                                              >  3:545427/1‑149 (MQ=255)
                                    gTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATc                                                                                                               >  4:403328/1‑146 (MQ=255)
                                    gTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATc                                                                                                               <  3:403328/146‑1 (MQ=255)
                                     ttACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGcc                                                                                                           <  3:48538/149‑1 (MQ=255)
                                     ttACACGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGcc                                                                                                           <  3:24205/149‑1 (MQ=255)
                                          cgcgcTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACAt                                                                                                      <  2:27634/149‑1 (MQ=255)
                                           gcgcTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATg                                                                                                     >  1:259152/1‑149 (MQ=255)
                                                 ccGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGAt                                                                                               >  3:501911/1‑149 (MQ=255)
                                                       acaATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTg                                                                                         <  4:545427/149‑1 (MQ=255)
                                                              gCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAg                                                                                  >  3:466619/1‑149 (MQ=255)
                                                                     gTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGa                                                                           >  2:51494/1‑149 (MQ=255)
                                                                     gTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGa                                                                           >  3:33905/1‑149 (MQ=255)
                                                                      ttGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGaa                                                                          <  2:76764/149‑1 (MQ=255)
                                                                        gAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTg                                                                                         >  2:95570/1‑132 (MQ=255)
                                                                        gAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGAGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTg                                                                                         <  1:95570/132‑1 (MQ=255)
                                                                         aaTTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGcc                                                                       >  2:51483/1‑149 (MQ=255)
                                                                           ttCAGTCCCGGCAAACCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCg                                                                     >  3:48885/1‑149 (MQ=255)
                                                                               gTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGa                                                                                                >  3:645700/1‑118 (MQ=255)
                                                                               gTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGa                                                                                                <  4:645700/118‑1 (MQ=255)
                                                                               gTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGAt                                                                                               >  2:303581/1‑119 (MQ=255)
                                                                               gTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGAt                                                                                               <  1:303581/119‑1 (MQ=255)
                                                                               gTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGt                                                                                          >  3:66788/1‑124 (MQ=255)
                                                                               gTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGt                                                                                          <  4:66788/124‑1 (MQ=255)
                                                                                tCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAAc                                                                >  4:229401/1‑149 (MQ=255)
                                                                                    ggCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaa                                                            >  4:5416/1‑149 (MQ=255)
                                                                                     gCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaa                                                           <  2:590653/149‑1 (MQ=255)
                                                                                       aaGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaaa                                                         >  4:97375/1‑149 (MQ=255)
                                                                                          ccTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAatat                                                      <  4:48885/149‑1 (MQ=255)
                                                                                          ccTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAatat                                                      <  2:11452/149‑1 (MQ=255)
                                                                                            tGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATca                                                    >  1:374197/1‑149 (MQ=255)
                                                                                            tGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATca                                                    >  4:20225/1‑149 (MQ=255)
                                                                                                     aGGTATGTCAGTTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGTAACTGAAAAATATCACCACCCGcc                                           >  3:455556/1‑149 (MQ=255)
                                                                                                       gTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAg                                         >  3:365897/1‑149 (MQ=255)
                                                                                                       gTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAg                                         <  4:33905/149‑1 (MQ=255)
                                                                                                            tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaaa                                                         >  1:302859/1‑128 (MQ=255)
                                                                                                            tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGaaaaa                                                         <  2:302859/128‑1 (MQ=255)
                                                                                                            tCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGata                                    >  2:633640/1‑149 (MQ=255)
                                                                                                              aGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATg                                  >  4:24898/1‑149 (MQ=255)
                                                                                                                 tCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACaa                                                                                   <  2:483680/97‑1 (MQ=255)
                                                                                                                 tCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACaa                                                                                   >  1:483680/1‑97 (MQ=255)
                                                                                                                   tGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGaa                             >  3:68295/1‑149 (MQ=255)
                                                                                                                   tGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGaa                             >  2:663960/1‑149 (MQ=255)
                                                                                                                     ccACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGccc                                                                      <  2:590616/106‑1 (MQ=255)
                                                                                                                     ccACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGccc                                                                      >  1:590616/1‑106 (MQ=255)
                                                                                                                          gggCAGACCAACGTTTGGCGCTGCGCAAAATGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATTTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCATCACCCGCCAGGGATATGCCGAAAGTCTGc                      <  3:412922/149‑1 (MQ=255)
                                                                                                                          gggCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGc                      <  3:593613/149‑1 (MQ=255)
                                                                                                                           ggCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGcc                     >  1:675117/1‑149 (MQ=255)
                                                                                                                               gACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCAttt                 <  4:101187/149‑1 (MQ=255)
                                                                                                                               gACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCAttt                 <  4:678165/149‑1 (MQ=255)
                                                                                                                                 ccAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGc               >  3:622616/1‑149 (MQ=255)
                                                                                                                                  cAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGcc              >  3:178451/1‑149 (MQ=255)
                                                                                                                                  cAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGcc              <  3:323497/149‑1 (MQ=255)
                                                                                                                                   aaCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCg             >  1:540020/1‑149 (MQ=255)
                                                                                                                                   aaCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCg             >  1:474542/1‑149 (MQ=255)
                                                                                                                                   aaCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCg             >  1:236660/1‑149 (MQ=255)
                                                                                                                                    aCGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGa            <  3:20225/149‑1 (MQ=255)
                                                                                                                                          ggCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGc      >  1:3026/1‑149 (MQ=255)
                                                                                                                                              cTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCAttt  >  4:318895/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTT  >  NC_000913/231821‑232111

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