Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,258 G→C 6.6% E641Q (GAG→CAG)  fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,2580GC 6.6% 527.0 / 23.4 151E641Q (GAG→CAG) fhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (58/83);  new base C (7/3);  total (65/86)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 151 total aligned reads displayed.

ATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGAT  >  NC_000913/2856111‑2856398
                                                                                                                                                   |                                                                                                                                            
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)
          cGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGAt                                                                                                                                   >  3:47182/1‑149 (MQ=255)
           gggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATg                                                                                                                                  >  3:1048690/1‑149 (MQ=255)
            ggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGa                                                                                                                                 <  1:209949/149‑1 (MQ=255)
              tAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAat                                                                                                                               <  4:888355/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)
                             ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgc                                                                                                                <  2:1656728/149‑1 (MQ=255)
                                   cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGc                                                                                                          <  2:535779/149‑1 (MQ=255)
                                       aTTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACTTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaa                                                                                                      <  4:274221/149‑1 (MQ=255)
                                        ttGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                     <  1:145578/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)
                                           aGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaa                                                                                                  <  3:77583/149‑1 (MQ=255)
                                            gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATGTTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                 <  1:1157048/149‑1 (MQ=255)
                                                gcgGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCaa                                                                                             <  3:668739/149‑1 (MQ=255)
                                                gcgGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCaa                                                                                             >  4:556607/1‑149 (MQ=255)
                                                   gTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACgg                                                                                          >  3:44340/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)
                                                          tAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGc                                                                                                          >  2:1295719/1‑126 (MQ=255)
                                                          tAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAAc                                                                                            <  2:1383005/140‑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)
                                                                 cgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGc                                                                            <  3:283185/149‑1 (MQ=255)
                                                                 cgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGc                                                                            <  3:1163941/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                                                                    >  3:38090/1‑149 (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                                                          <  4:44340/149‑1 (MQ=255)
                                                                                    tCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACg                                                         >  3:893560/1‑149 (MQ=255)
                                                                                     cATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGt                                                        <  4:1005743/149‑1 (MQ=255)
                                                                                       ttGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTAGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCt                                                      >  3:473544/1‑149 (MQ=255)
                                                                                         gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                    <  1:1159851/149‑1 (MQ=255)
                                                                                         gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                    <  2:790815/149‑1 (MQ=255)
                                                                                         gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                    <  3:1560680/149‑1 (MQ=255)
                                                                                         gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                    <  4:139773/149‑1 (MQ=255)
                                                                                          ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAAc                                                                                            >  1:959960/1‑108 (MQ=255)
                                                                                          ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAAc                                                                                            <  2:959960/108‑1 (MQ=255)
                                                                                          ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                   >  1:396572/1‑149 (MQ=255)
                                                                                           cAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgggg                                                  <  2:396572/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                                                                                                     <  4:777624/93‑1 (MQ=255)
                                                                                                atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                     >  3:777624/1‑93 (MQ=255)
                                                                                                atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTga                                                                                                       <  1:700799/91‑1 (MQ=255)
                                                                                                 ttGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGaaa                                            <  3:1457581/149‑1 (MQ=255)
                                                                                                  tGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAAc                                           <  3:529475/149‑1 (MQ=255)
                                                                                                  tGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAAc                                           <  2:1370029/149‑1 (MQ=255)
                                                                                                    ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGc                                                                                 >  4:1006143/1‑109 (MQ=255)
                                                                                                    ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGc                                                                                 <  3:1006143/109‑1 (MQ=255)
                                                                                                    ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                         <  4:1048690/149‑1 (MQ=255)
                                                                                                    ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                         >  4:113550/1‑149 (MQ=255)
                                                                                                    ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                         >  1:330629/1‑149 (MQ=255)
                                                                                                    ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                         >  1:423154/1‑149 (MQ=255)
                                                                                                     tttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCa                                        <  3:556607/149‑1 (MQ=255)
                                                                                                         ccGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                    >  3:919283/1‑149 (MQ=255)
                                                                                                               ccTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgct                              <  4:621116/149‑1 (MQ=255)
                                                                                                                 tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                            >  1:1010870/1‑149 (MQ=255)
                                                                                                                     aCGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACg                        >  4:1193260/1‑149 (MQ=255)
                                                                                                                       gccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCa                                        >  2:678342/1‑131 (MQ=255)
                                                                                                                       gccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCa                                        <  1:678342/131‑1 (MQ=255)
                                                                                                                        ccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCCAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGAt                     <  3:673877/149‑1 (MQ=255)
                                                                                                                         cgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATg                    <  3:923299/149‑1 (MQ=255)
                                                                                                                             tGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGc                >  2:1125243/1‑149 (MQ=255)
                                                                                                                                 gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAtt                                                                                                                    >  3:747440/1‑45 (MQ=255)
                                                                                                                                 gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAtt                                                                                                                    <  4:747440/45‑1 (MQ=255)
                                                                                                                                      ggTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa       >  1:1631547/1‑149 (MQ=255)
                                                                                                                                       gTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAt      <  3:1268396/149‑1 (MQ=255)
                                                                                                                                       atcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAt      >  3:1447514/5‑149 (MQ=255)
                                                                                                                                        tcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt     >  1:1239213/4‑149 (MQ=255)
                                                                                                                                         cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAg                                                                      <  2:1508615/81‑1 (MQ=255)
                                                                                                                                         cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAg                                                                      >  1:1508615/3‑83 (MQ=255)
                                                                                                                                         cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTg    >  1:736693/3‑149 (MQ=255)
                                                                                                                                          ttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTgg          <  2:230159/141‑1 (MQ=255)
                                                                                                                                          ttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa   >  4:1225444/2‑149 (MQ=255)
                                                                                                                                          gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGcc                                                                                <  4:1387608/72‑1 (MQ=255)
                                                                                                                                          gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa   <  2:1013954/149‑1 (MQ=255)
                                                                                                                                             gCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                    >  3:117854/1‑97 (MQ=255)
                                                                                                                                             gCCCTGGAGGGCGAAGATGAATATCAGTTCATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                    <  4:117854/97‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
ATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGAT  >  NC_000913/2856111‑2856398

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: