Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,849,945 A→T 5.0% intergenic (‑96/+29) hycB ← / ← hycA hydrogenase 3, Fe‑S subunit/regulator of the transcriptional regulator FhlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,849,9450AT 5.0% 418.8 / 8.1 118intergenic (‑96/+29)hycB/hycAhydrogenase 3, Fe‑S subunit/regulator of the transcriptional regulator FhlA
Reads supporting (aligned to +/- strand):  ref base A (60/52);  new base T (2/4);  total (62/56)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

AAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAA  >  NC_000913/2849798‑2850092
                                                                                                                                                   |                                                                                                                                                   
aaGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTaa                                                                                                                                                    <  3:234065/149‑1 (MQ=255)
aaGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTaa                                                                                                                                                    >  3:540275/1‑149 (MQ=255)
 aGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAAc                                                                                                                                                   >  3:15184/1‑149 (MQ=255)
     tgGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCt                                                                                                                                               <  4:15184/149‑1 (MQ=255)
        cAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtt                                                                                                                                            <  1:114055/149‑3 (MQ=255)
        cAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCa                                                                                                                                            >  2:376036/1‑149 (MQ=255)
           ccGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttgg                                                                                                                                         <  3:779262/149‑6 (MQ=255)
           ccGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttgg                                                                                                                                         <  3:842464/149‑6 (MQ=255)
              aTACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACgg                                                                                                                                      <  1:164248/149‑1 (MQ=255)
              aTACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACgg                                                                                                                                      <  2:438354/149‑1 (MQ=255)
              aTACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACgg                                                                                                                                      <  4:146839/149‑1 (MQ=255)
                aCAGAGCGTGGAGTCAGCAATTACAAAACGGTTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGc                                                                                                                                    >  2:315433/1‑149 (MQ=255)
                  agagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaa                                                                                                                                  <  2:2271/149‑1 (MQ=255)
                    agCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAg                                                                                                                                >  1:368028/1‑149 (MQ=255)
                         ggAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCa                                                                                                                           <  4:410576/149‑1 (MQ=255)
                             tCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTc                                                                                                                       >  2:783575/1‑149 (MQ=255)
                             tCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTc                                                                                                                       >  4:800383/1‑149 (MQ=255)
                               aGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAt                                                                                                                     >  1:382448/1‑149 (MQ=255)
                                     tACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATCCTGcc                                                                                                               <  4:540275/149‑1 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCTTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                            >  4:598618/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                            >  4:273498/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                            >  3:402907/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                            >  2:89216/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                            >  2:557249/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                            >  2:78083/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                            >  1:308462/1‑149 (MQ=255)
                                          aaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgttt                                                                                                                                           >  2:514799/1‑113 (MQ=255)
                                          aaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                          >  4:818377/1‑149 (MQ=255)
                                          aaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                          >  1:71788/1‑149 (MQ=255)
                                             gATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTg                                                                                                       >  2:657809/1‑149 (MQ=255)
                                              aTTGACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGt                                                                                                      <  2:233243/149‑1 (MQ=255)
                                              aTTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGc                                                                                                                         <  4:410914/130‑1 (MQ=255)
                                              aTTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGc                                                                                                                         >  3:410914/1‑130 (MQ=255)
                                              aTTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGt                                                                                                      <  3:5158/149‑1 (MQ=255)
                                               ttCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTc                                                                                                     <  4:519056/149‑1 (MQ=255)
                                                tCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTcc                                                                                                    <  1:873785/149‑1 (MQ=255)
                                                tCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTcc                                                                                                    <  2:689377/149‑1 (MQ=255)
                                                  aCCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCt                                                                                                  >  3:139160/1‑149 (MQ=255)
                                                  aCCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCt                                                                                                  >  2:628854/1‑149 (MQ=255)
                                                      ggCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAACCAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTAc                                                                                              <  2:553484/149‑1 (MQ=255)
                                                         aTTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAg                                                                                           >  3:857862/1‑149 (MQ=255)
                                                              tCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCAGCCCCAAACAAGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCGCTTTACCAGTTGGc                                                                                      >  2:536064/1‑149 (MQ=255)
                                                                aGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCtt                                                                                    <  4:503499/149‑1 (MQ=255)
                                                                   tGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTaaa                                                                                 <  4:136352/149‑1 (MQ=255)
                                                                   tGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTaaa                                                                                 <  3:181229/149‑1 (MQ=255)
                                                                         tCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCAcc                                                                           >  1:412106/1‑149 (MQ=255)
                                                                          cATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAGACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCg                                                                          <  2:812566/149‑1 (MQ=255)
                                                                          cATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCg                                                                          <  1:583505/149‑1 (MQ=255)
                                                                            tttttGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGc                                                                        >  2:307856/1‑149 (MQ=255)
                                                                             ttttGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGca                                                                       >  2:434417/1‑149 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                                <  2:221329/149‑1 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                                >  2:136212/1‑149 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                                >  4:255543/1‑149 (MQ=255)
                                                                                     aaaaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGttt                                                               <  1:628854/149‑1 (MQ=255)
                                                                                      aaaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtttt                                                              <  4:172395/149‑1 (MQ=255)
                                                                                      aaaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtttt                                                              <  3:654396/149‑1 (MQ=255)
                                                                                      aaaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtttt                                                              <  3:255761/149‑1 (MQ=255)
                                                                                              cGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtttt                                                              >  2:358253/1‑141 (MQ=255)
                                                                                              cGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCg                                                      <  1:691785/149‑1 (MQ=255)
                                                                                              cGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCg                                                      <  2:524431/149‑1 (MQ=255)
                                                                                               gATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCgg                                                     >  4:694989/1‑149 (MQ=255)
                                                                                                  gAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgtgt                                                  >  4:127467/1‑149 (MQ=255)
                                                                                                      tGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATaa                                              <  3:371340/149‑1 (MQ=255)
                                                                                                       gTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAg                                             >  1:528388/1‑149 (MQ=255)
                                                                                                       gTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAg                                             <  2:382448/149‑1 (MQ=255)
                                                                                                          aTTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTg                                          >  3:253410/1‑149 (MQ=255)
                                                                                                          aTTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTg                                          <  1:536064/149‑1 (MQ=255)
                                                                                                           ttttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGt                                         >  3:634106/1‑149 (MQ=255)
                                                                                                           ttttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGt                                         >  4:285738/1‑149 (MQ=255)
                                                                                                            tttCGACACTCATCGACACGCCCATCCCCAAACAGGTGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTc                                        <  4:682042/149‑1 (MQ=255)
                                                                                                            tttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTc                                        <  3:800383/149‑1 (MQ=255)
                                                                                                            tttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTc                                        <  2:80199/149‑1 (MQ=255)
                                                                                                             ttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCg                                       <  1:350181/149‑1 (MQ=255)
                                                                                                              tCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGt                                      >  3:709339/1‑149 (MQ=255)
                                                                                                                   acTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACg                                 >  2:1290/1‑149 (MQ=255)
                                                                                                                     tCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGcc                               >  2:136585/1‑149 (MQ=255)
                                                                                                                     tCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGcc                               >  2:625833/1‑149 (MQ=255)
                                                                                                                      cATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCa                              >  1:529698/1‑149 (MQ=255)
                                                                                                                          gACACGCCCACCACCAGACAGGCGTAAAGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGTTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACAATGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGa                          <  4:346161/149‑1 (MQ=255)
                                                                                                                            cacGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAtt                        <  3:273498/149‑1 (MQ=255)
                                                                                                                               gCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAttttt                     >  1:420749/1‑149 (MQ=255)
                                                                                                                                cccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAtttttt                    <  1:376036/149‑1 (MQ=255)
                                                                                                                                cccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAtttttt                    >  1:783034/1‑149 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                   <  4:139160/149‑1 (MQ=255)
                                                                                                                                     ccccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCgaga               <  1:315433/149‑1 (MQ=255)
                                                                                                                                      cccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGAt              <  3:694989/149‑1 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATaa                                              <  1:580874/116‑1 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATaa                                              >  2:580874/1‑116 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                   >  2:819772/1‑143 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                   <  1:819772/143‑1 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa             >  2:873813/1‑149 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa             >  1:619933/1‑149 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa             >  1:546952/1‑149 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa             >  1:711267/1‑149 (MQ=255)
                                                                                                                                         aaaCAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAgt           <  1:568732/149‑1 (MQ=255)
                                                                                                                                         aaaCAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAgt           <  3:818377/149‑1 (MQ=255)
                                                                                                                                            cAGGCGTTACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGtt        >  1:52810/1‑149 (MQ=255)
                                                                                                                                              ggCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTcc      <  3:367538/149‑1 (MQ=255)
                                                                                                                                                 gTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGa   >  3:304488/1‑149 (MQ=255)
                                                                                                                                                  tAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCTAGCGATTTTTTCGAGATAGTGTTCCAGaa  <  1:747396/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAA  >  NC_000913/2849798‑2850092

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