Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,144,083 G→T 9.2% S30I (AGC→ATC)  pflD → putative glycine radical domain‑containing pyruvate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,144,0830GT 9.2% 361.7 / 21.0 109S30I (AGC→ATC) pflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (47/52);  new base T (7/3);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 113 total aligned reads displayed.

ACCGCTTAAGTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTATGTCGCC  >  NC_000913/4143936‑4144230
                                                                                                                                                   |                                                                                                                                                   
aCCGCTTAAGTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCgccagc                                                                                                                                                    >  4:882335/1‑149 (MQ=255)
    cTTAAGTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATc                                                                                                                                                <  4:915462/149‑1 (MQ=255)
    cTTAAGTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATc                                                                                                                                                <  2:209909/149‑1 (MQ=255)
     ttAAGTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCg                                                                                                                                               >  1:714928/1‑149 (MQ=255)
         gTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCATCCATCGGCaa                                                                                                                                           >  4:45820/1‑149 (MQ=255)
         gTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCaa                                                                                                                                           <  3:193570/149‑1 (MQ=255)
         gTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTAGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCaa                                                                                                                                           >  1:667368/1‑149 (MQ=255)
          ttcttcTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCaaa                                                                                                                                          <  3:674026/149‑1 (MQ=255)
             ttctCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCg                                                                                                                                       >  4:512720/1‑149 (MQ=255)
               ctcCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCACTACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGaa                                                                                                                                     <  4:881208/149‑1 (MQ=255)
                tcCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAg                                                                                                                                    <  3:282958/149‑1 (MQ=255)
                tcCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAg                                                                                                                                    >  1:407776/1‑149 (MQ=255)
                  ccGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTGTCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCATCCATCGGCAAACCGAAGGc                                                                                                                                  >  2:800511/1‑149 (MQ=255)
                  ccGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGc                                                                                                                                  >  4:152729/1‑149 (MQ=255)
                       cGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAAcc                                                                                                                             <  3:60101/149‑1 (MQ=255)
                        gCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCg                                                                                                                            >  2:589089/1‑149 (MQ=255)
                          aaGGGCGGGTTCGCTTTCCCACAGGAGGTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGt                                                                                                                          >  4:722100/1‑149 (MQ=255)
                           aGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTg                                                                                                                         >  4:150327/1‑149 (MQ=255)
                            gggCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGa                                                                                                                        <  3:654321/149‑1 (MQ=255)
                            gggCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGa                                                                                                                        >  1:352253/1‑149 (MQ=255)
                            gggCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGa                                                                                                                        <  2:853850/149‑1 (MQ=255)
                             ggCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGat                                                                                                                       <  4:622003/149‑1 (MQ=255)
                              gCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGata                                                                                                                      <  4:854531/149‑1 (MQ=255)
                                   ttCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTgcg                                                                                                                 >  2:842154/1‑149 (MQ=255)
                                   ttCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTgcg                                                                                                                 >  4:679346/1‑149 (MQ=255)
                                      gCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCg                                                                                                              >  4:916212/1‑149 (MQ=255)
                                        tttCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCggg                                                                                                            >  2:749175/1‑149 (MQ=255)
                                              acaGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCGGGGCGAAAg                                                                                                      >  1:84930/1‑149 (MQ=255)
                                                     ttCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGc                                                                                               >  2:356956/1‑149 (MQ=255)
                                                     ttCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGc                                                                                               >  3:602158/1‑149 (MQ=255)
                                                        cTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGta                                                                                            <  3:160804/149‑1 (MQ=255)
                                                         tCATGACGAATCGTATCTCTCGGCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGGGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGtat                                                                                           >  4:785597/1‑149 (MQ=255)
                                                         tCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCATCCATCGGCAAACCGAAGGCGAACCGGTGCTATTGCGCCGGGCGAAAGCAACAGCGtat                                                                                           >  3:195876/1‑149 (MQ=255)
                                                         tCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGtat                                                                                           >  1:16720/1‑149 (MQ=255)
                                                               cGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCtt                                                                                     <  4:584680/149‑1 (MQ=255)
                                                               cGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGAGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCtt                                                                                     <  3:630499/149‑1 (MQ=255)
                                                                gAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGt                                                                                             <  4:244531/140‑1 (MQ=255)
                                                                gAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTg                                                                                    <  2:16720/149‑1 (MQ=255)
                                                                     gTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACAt                                                                               >  4:172910/1‑149 (MQ=255)
                                                                      tATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGCCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCATTGAACATg                                                                              >  4:688775/1‑149 (MQ=255)
                                                                       aTCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCATCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGt                                                                             <  1:452310/149‑1 (MQ=255)
                                                                       aTCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGc                                                                                                           >  1:173264/1‑119 (MQ=255)
                                                                       aTCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGc                                                                                                           <  2:173264/119‑1 (MQ=255)
                                                                        tctctcGCCGCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCACTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAAGAGCGTATATCCTTGAACATGtt                                                                            >  1:802598/1‑149 (MQ=255)
                                                                         ctctcGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATg                                                                              <  1:483594/146‑1 (MQ=255)
                                                                         ctctcGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATg                                                                              >  2:483594/1‑146 (MQ=255)
                                                                                cTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTc                                                                    <  4:602158/149‑1 (MQ=255)
                                                                                cTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTc                                                                    >  1:809805/1‑149 (MQ=255)
                                                                                 tCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCg                                                                   >  3:686421/1‑149 (MQ=255)
                                                                                 tCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCg                                                                   >  3:329582/1‑149 (MQ=255)
                                                                                  cAAAACTGCACTGTTTGCCAATACCCGTGAAAACTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAAGATGTTGAAATTTCGa                                                                  <  2:424246/149‑1 (MQ=255)
                                                                                   aaaaCTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGAt                                                                 >  1:312537/1‑149 (MQ=255)
                                                                                   aaaaCTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGAt                                                                 <  3:761275/149‑1 (MQ=255)
                                                                                     aaCTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCATCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTc                                                               >  4:841398/1‑149 (MQ=255)
                                                                                     aaCTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTc                                                               <  1:344300/149‑1 (MQ=255)
                                                                                     aaCTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTc                                                               <  3:150327/149‑1 (MQ=255)
                                                                                       cggCACTGGTTGCCGTTATGCGTGAGATCTCGCTGGAGCGGGCGCTGGTTTATTCCTCCAGCGATCGGCCAGCCGAATGCGAACCGGTGATATTGCGCCGGACGAAAGCAACAGCGTATATCCTTGAACATGTTTAAACGTCGATTTGt                                                             <  2:40138/147‑1 (MQ=255)
                                                                                       cTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGt                                                             <  1:540130/149‑1 (MQ=255)
                                                                                          cACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGtgat                                                          <  3:11706/149‑1 (MQ=255)
                                                                                          cACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGtgat                                                          >  4:112415/1‑149 (MQ=255)
                                                                                          cACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGtgat                                                          <  1:842154/149‑1 (MQ=255)
                                                                                             tGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATgaa                                                       >  2:24582/1‑149 (MQ=255)
                                                                                              gTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCATCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATgaag                                                      >  3:235153/1‑149 (MQ=255)
                                                                                                  gCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGtga                                                           >  1:799590/1‑140 (MQ=255)
                                                                                                  gCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGtga                                                           <  2:799590/140‑1 (MQ=255)
                                                                                                      aTACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCATCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGAtt                                              <  3:45820/149‑1 (MQ=255)
                                                                                                      aTACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCGGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGAtt                                              <  2:84930/149‑1 (MQ=255)
                                                                                                      aTACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGAtt                                              <  4:528899/149‑1 (MQ=255)
                                                                                                      aTACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGAtt                                              <  4:329582/149‑1 (MQ=255)
                                                                                                        aCCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGc                                            <  3:882335/149‑1 (MQ=255)
                                                                                                         cccGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGcc                                           <  3:302724/149‑1 (MQ=255)
                                                                                                          ccGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCg                                          >  4:499578/1‑149 (MQ=255)
                                                                                                           cGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCgg                                         >  3:563014/1‑149 (MQ=255)
                                                                                                           cGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCgg                                         >  3:211116/1‑149 (MQ=255)
                                                                                                                   ctcGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCa                                 >  3:331393/1‑149 (MQ=255)
                                                                                                                    tcGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCAc                                >  3:281334/1‑149 (MQ=255)
                                                                                                                    tcGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCAc                                >  4:451133/1‑149 (MQ=255)
                                                                                                                      gCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCg                              >  4:740571/1‑149 (MQ=255)
                                                                                                                           aGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTaaaa                         <  2:809805/149‑1 (MQ=255)
                                                                                                                            gCGGGCGCTGCTTTATACCGCCATCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAAc                        >  1:744254/1‑149 (MQ=255)
                                                                                                                               ggCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACcgc                     >  3:30190/1‑149 (MQ=255)
                                                                                                                               ggCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGATGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACcgc                     <  1:504717/149‑1 (MQ=255)
                                                                                                                                  gctgctTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGTGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACcgcgcg                  >  2:297952/1‑149 (MQ=255)
                                                                                                                                   ctgctTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCgcgcgc                 <  1:629204/149‑1 (MQ=255)
                                                                                                                                    tgctTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGcc                >  1:276310/1‑149 (MQ=255)
                                                                                                                                      ctTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCgg              >  4:673155/1‑149 (MQ=255)
                                                                                                                                       tttATACCGCCATCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCg                                          <  1:342180/120‑1 (MQ=255)
                                                                                                                                       tttATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCggg             >  2:509152/1‑149 (MQ=255)
                                                                                                                                        ttATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGa            <  2:352253/149‑1 (MQ=255)
                                                                                                                                         tataCCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGAt           <  4:211116/149‑1 (MQ=255)
                                                                                                                                          ataCCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGACCTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGAtt          <  3:499578/149‑1 (MQ=255)
                                                                                                                                          ataCCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGAtt          <  1:524935/149‑1 (MQ=255)
                                                                                                                                          ataCCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGAtt          >  3:837135/1‑149 (MQ=255)
                                                                                                                                          ataCCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAAACGCACCGTAAAACCGCGCGCCGGGAtt          <  1:53505/149‑1 (MQ=255)
                                                                                                                                            aCCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTAt        <  2:923142/149‑1 (MQ=255)
                                                                                                                                            aCCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTAt        >  1:479106/1‑149 (MQ=255)
                                                                                                                                             ccGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTATg       <  2:554204/149‑1 (MQ=255)
                                                                                                                                                ccagccaTCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTATGTCg    >  3:409614/1‑149 (MQ=255)
                                                                                                                                                 cagccaTCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTATGTCGc   <  1:337642/149‑1 (MQ=255)
                                                                                                                                                  agccaTCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTATGTCGcc  >  1:775193/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCGCTTAAGTTCTTCTCCCGCTCGCAAGGGCGGGTTCGCTTTCCCACAGGAGTTCCTCATGACGAATCGTATCTCTCGCCTCAAAACTGCACTGTTTGCCAATACCCGTGAAATCTCGCTGGAGCGGGCGCTGCTTTATACCGCCAGCCATCGGCAAACCGAAGGCGAACCGGTGATATTGCGCCGGGCGAAAGCAACAGCGTATATCCTTGAACATGTTGAAATTTCGATTCGTGATGAAGAACTGATTGCCGGTAACCGCACCGTAAAACCGCGCGCCGGGATTATGTCGCC  >  NC_000913/4143936‑4144230

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