Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,257 G→C 6.6% L640L (CTG→CTC fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,2570GC 6.6% 521.1 / 22.7 151L640L (CTG→CTCfhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (57/83);  new base C (7/3);  total (65/86)
Fisher's exact test for biased strand distribution p-value = 9.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 152 total aligned reads displayed.

ACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAA  >  NC_000913/2856109‑2856400
                                                                                                                                                    |                                                                                                                                               
aCATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTg                                                                                                                                                 >  2:209949/1‑149 (MQ=255)
  aTGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGa                                                                                                                                               <  1:788346/149‑1 (MQ=255)
  aTGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGa                                                                                                                                               <  4:1486713/149‑1 (MQ=255)
  aTGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGa                                                                                                                                               >  1:30145/1‑149 (MQ=255)
         ggCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGaa                                                                                                                                        <  2:696124/149‑1 (MQ=255)
            cGGGTAACGTACGCGATCTGGAAAACGTCATTGAGCGCGCGGTATTGCTATCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAATCGCCGCCTGGCGCAACGGTTGTCGCCCTGGAGGGCGAAGAt                                                                                                                                     <  2:284000/149‑1 (MQ=255)
             gggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATg                                                                                                                                    >  3:1048690/1‑149 (MQ=255)
              ggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGa                                                                                                                                   <  1:209949/149‑1 (MQ=255)
               gTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                  <  2:30145/149‑1 (MQ=255)
                tAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAat                                                                                                                                 <  2:1482519/149‑1 (MQ=255)
                tAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAat                                                                                                                                 <  4:888355/149‑1 (MQ=255)
                     ttcacgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAg                                                                                                                            <  3:1690069/145‑1 (MQ=255)
                       cgcgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGtt                                                                                                                          <  3:348513/149‑1 (MQ=255)
                        gcgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTg                                                                                                                         >  4:942998/1‑149 (MQ=255)
                         cgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGa                                                                                                                        <  3:942998/149‑1 (MQ=255)
                          gAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAt                                                                                                                       <  3:750417/149‑1 (MQ=255)
                              tGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATtgtg                                                                                                                   <  4:1642898/149‑1 (MQ=255)
                              tGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATtgtg                                                                                                                   >  1:994028/1‑149 (MQ=255)
                                   aaCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGt                                                                                                              <  3:983639/149‑1 (MQ=255)
                                      gTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCt                                                                                                           >  3:505811/1‑149 (MQ=255)
                                         aTTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACTTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaa                                                                                                        <  4:274221/149‑1 (MQ=255)
                                          ttGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                       <  1:145578/149‑1 (MQ=255)
                                          ttGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                       <  2:1346309/149‑1 (MQ=255)
                                          ttGAGCGCGAGGTATTGATAACACGCGGTAACGTGCTGCAGCTGTCATTTCCAGATATTGTTTTACCGGAACCTAAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                       <  1:1547330/149‑1 (MQ=255)
                                            gAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaaga                                                                                                     >  2:1106387/1‑149 (MQ=255)
                                            gAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaaga                                                                                                     <  4:47182/149‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATGTTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                   <  1:1157048/149‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                   <  2:22743/149‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                   <  2:994028/149‑1 (MQ=255)
                                                  gcgGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCaa                                                                                               >  4:556607/1‑149 (MQ=255)
                                                     gTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACgg                                                                                            >  3:44340/1‑149 (MQ=255)
                                                     gTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACgg                                                                                            >  3:78869/1‑149 (MQ=255)
                                                         tGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTg                                                                                        <  3:2905/149‑1 (MQ=255)
                                                            tAACACGCGGTAACGTGCTGCAGCTGTCGTTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAATGTGCGCGTGCTGAAAGAAACCAAc                                                                                              >  1:1383005/1‑140 (MQ=255)
                                                            tAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATTAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGtt                                                                                     >  4:1241617/1‑149 (MQ=255)
                                                            tAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGc                                                                                                            <  1:1295719/126‑1 (MQ=255)
                                                            tAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGtt                                                                                     <  4:78869/149‑1 (MQ=255)
                                                             aaCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATCGTGCGCGCGCTGAAAGAAACCAACGGCGTGGTTg                                                                                    >  3:1380017/1‑149 (MQ=255)
                                                              acacGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGc                                                                                   <  4:269902/149‑1 (MQ=255)
                                                               cacGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGcc                                                                                  >  1:932219/1‑149 (MQ=255)
                                                                 cgcgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGCGCTGAAAGAAACCAACGGCGTGGTTGCCgg                                                                                <  4:1380017/149‑1 (MQ=255)
                                                                    ggTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                             <  3:1684200/149‑1 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                       >  4:579084/1‑149 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAAGGCGAAGTTGACTATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                       >  2:1227641/1‑149 (MQ=255)
                                                                           tgctgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGgc                                                                      <  2:665271/149‑1 (MQ=255)
                                                                             ctgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGgcgc                                                                    >  2:51824/1‑149 (MQ=255)
                                                                              tgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCt                                                                   <  3:1290593/149‑1 (MQ=255)
                                                                                cAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgc                                                                 <  4:88604/149‑1 (MQ=255)
                                                                                 aGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgcg                                                                >  4:1349451/1‑149 (MQ=255)
                                                                                   cTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAAc                                                                                                  <  1:107679/113‑1 (MQ=255)
                                                                                     gTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAAc                                                            <  3:1256611/149‑1 (MQ=255)
                                                                                     gTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAAc                                                            <  4:44340/149‑1 (MQ=255)
                                                                                         ttGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTAGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCt                                                        >  3:473544/1‑149 (MQ=255)
                                                                                           gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                      <  1:1159851/149‑1 (MQ=255)
                                                                                           gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                      <  4:139773/149‑1 (MQ=255)
                                                                                           accAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                      <  3:147107/148‑1 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAAc                                                                                              >  1:959960/1‑108 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                     >  1:396572/1‑149 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                     <  2:9439/149‑1 (MQ=255)
                                                                                               gATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGAGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCt                                                  >  2:57517/1‑149 (MQ=255)
                                                                                                 tatTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                  <  2:983672/67‑1 (MQ=255)
                                                                                                 tatTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTGGCCGGGCCTAAAGGCGCTGCGCAACGACTGGGGCTGa                                                >  2:83515/1‑149 (MQ=255)
                                                                                                  atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                       >  3:777624/1‑93 (MQ=255)
                                                                                                  atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                       <  4:777624/93‑1 (MQ=255)
                                                                                                  atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTga                                                                                                         <  1:700799/91‑1 (MQ=255)
                                                                                                  atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTga                                                                                                         >  2:700799/1‑91 (MQ=255)
                                                                                                   ttGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGaaa                                              <  3:1457581/149‑1 (MQ=255)
                                                                                                    tGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAAc                                             >  3:621116/1‑149 (MQ=255)
                                                                                                    tGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAAc                                             <  2:1370029/149‑1 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGc                                                                                   >  4:1006143/1‑109 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                           >  1:330629/1‑149 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                           >  1:423154/1‑149 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                           >  4:113550/1‑149 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                           <  4:1048690/149‑1 (MQ=255)
                                                                                                       tttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCa                                          <  3:556607/149‑1 (MQ=255)
                                                                                                           ccGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                      >  3:919283/1‑149 (MQ=255)
                                                                                                            cGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTc                                                         >  3:755013/1‑129 (MQ=255)
                                                                                                                 ccTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgct                                <  4:621116/149‑1 (MQ=255)
                                                                                                                  cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctg                               >  3:921339/1‑149 (MQ=255)
                                                                                                                   tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                              >  1:1010870/1‑149 (MQ=255)
                                                                                                                    gAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTc                             >  2:52722/1‑149 (MQ=255)
                                                                                                                       aCGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACg                          >  4:1193260/1‑149 (MQ=255)
                                                                                                                         gccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCa                                          <  1:678342/131‑1 (MQ=255)
                                                                                                                          ccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCCAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGAt                       <  3:673877/149‑1 (MQ=255)
                                                                                                                               tGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGc                  <  3:471567/149‑1 (MQ=255)
                                                                                                                                   gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAtt                                                                                                                      <  4:747440/45‑1 (MQ=255)
                                                                                                                                        ggTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa         >  1:1631547/1‑149 (MQ=255)
                                                                                                                                         gTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAt        >  3:737396/1‑149 (MQ=255)
                                                                                                                                          tcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt       >  1:1239213/4‑149 (MQ=255)
                                                                                                                                           cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAg                                                                        >  1:1508615/3‑83 (MQ=255)
                                                                                                                                           cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAg                                                                        <  2:1508615/81‑1 (MQ=255)
                                                                                                                                           cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTg      >  3:1059567/3‑149 (MQ=255)
                                                                                                                                           cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTg      >  1:736693/3‑149 (MQ=255)
                                                                                                                                            ttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa     >  4:1225444/2‑149 (MQ=255)
                                                                                                                                            gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGcc                                                                                  <  4:1387608/72‑1 (MQ=255)
                                                                                                                                            gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGcc                                                                                  >  3:1387608/1‑72 (MQ=255)
                                                                                                                                            gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa     <  2:1013954/149‑1 (MQ=255)
                                                                                                                                               gCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                      >  3:117854/1‑97 (MQ=255)
                                                                                                                                               gCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtgacagc  <  2:1170593/149‑7 (MQ=255)
                                                                                                                                               gCCCTGGAGGGCGAAGATGAATATCAGTTCATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                      <  4:117854/97‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAA  >  NC_000913/2856109‑2856400

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: