Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,849,945 A→T 6.3% 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 6.3% 328.1 / 9.1 95intergenic (‑96/+29)hycB/hycAhydrogenase 3, Fe‑S subunit/regulator of the transcriptional regulator FhlA
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base T (3/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

AAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCA  >  NC_000913/2849798‑2850089
                                                                                                                                                   |                                                                                                                                                
aaGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTaa                                                                                                                                                 <  3:125852/149‑1 (MQ=255)
     tgGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCt                                                                                                                                            <  4:95874/149‑1 (MQ=255)
     tgGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCt                                                                                                                                            <  2:294142/149‑1 (MQ=255)
     tgGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCt                                                                                                                                            <  3:469924/149‑1 (MQ=255)
       gCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGc                                                                                                                                          >  3:134642/1‑149 (MQ=255)
           ccGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttgg                                                                                                                                      <  3:504605/149‑6 (MQ=255)
             gATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCTTCAGGTGATTGTAATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACg                                                                                                                                    >  4:556333/1‑149 (MQ=255)
             gATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACg                                                                                                                                    <  1:322903/149‑1 (MQ=255)
             gATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACg                                                                                                                                    <  3:262943/149‑1 (MQ=255)
               tACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACggg                                                                                                                                  >  3:579949/1‑149 (MQ=255)
                aCAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGc                                                                                                                                 >  2:95325/1‑149 (MQ=255)
                  agagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaa                                                                                                                               >  1:612864/1‑149 (MQ=255)
                   gagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtt                                                                                                                                         >  3:538410/1‑136 (MQ=255)
                   gagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtt                                                                                                                                         <  4:538410/138‑3 (MQ=255)
                   gagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaaa                                                                                                                              <  3:560970/149‑1 (MQ=255)
                   gagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaaa                                                                                                                              >  3:473661/1‑149 (MQ=255)
                        tGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTc                                                                                                                         >  4:123085/1‑149 (MQ=255)
                         ggAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCa                                                                                                                        <  1:573679/149‑1 (MQ=255)
                         ggAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCa                                                                                                                        <  1:588258/149‑1 (MQ=255)
                             tCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTc                                                                                                                    >  2:231711/1‑149 (MQ=255)
                               aGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAt                                                                                                                  <  3:123085/149‑1 (MQ=255)
                                gCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtg                                                                                                                 >  1:3790/1‑149 (MQ=255)
                                  aaTTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgct                                                                                                               >  3:636708/1‑149 (MQ=255)
                                  aaTTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgct                                                                                                               >  3:525951/1‑149 (MQ=255)
                                    ttACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgctgc                                                                                                             >  1:696798/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                         >  1:397629/1‑149 (MQ=255)
                                         aaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGc                                                                                                        <  1:235601/149‑1 (MQ=255)
                                          aaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                       >  4:592238/1‑149 (MQ=255)
                                          aaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                       <  1:755855/149‑1 (MQ=255)
                                           aCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCtt                                                                                                      <  3:423356/149‑1 (MQ=255)
                                               ttCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTc                                                                                                  >  4:535578/1‑149 (MQ=255)
                                                tCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTcc                                                                                                 <  2:373289/149‑1 (MQ=255)
                                                 cACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTccc                                                                                                <  1:728682/149‑1 (MQ=255)
                                                 cACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTccc                                                                                                <  3:51309/149‑1 (MQ=255)
                                                     aGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTa                                                                                            >  1:644655/1‑149 (MQ=255)
                                                           tCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGtt                                                                                      <  1:602448/149‑1 (MQ=255)
                                                                   tGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACgg                                                                                                                                   <  4:703362/96‑1 (MQ=255)
                                                                   tGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACgg                                                                                                                                   >  3:703362/1‑96 (MQ=255)
                                                                   tGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTaaa                                                                              <  1:235118/149‑1 (MQ=255)
                                                                     aTTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAAtt                                                                            >  1:447504/1‑149 (MQ=255)
                                                                     aTTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAAtt                                                                            >  4:497454/1‑149 (MQ=255)
                                                                       tGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCa                                                                          <  3:35850/149‑1 (MQ=255)
                                                                         tCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACAATCATCGACACGCCAATCCCCAAACAGGCGTAACGCCTGCAAAACGGGAAAAGCCTCAGCTCATGATACCGGGCTTTGTCCCTTTACCAGCTGGCTTGAATTCAcc                                                                        >  2:228012/1‑149 (MQ=255)
                                                                           aTTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCgg                                                                      >  4:170740/1‑149 (MQ=255)
                                                                                tGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACAt                                                                 >  3:453114/1‑149 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaaa                                                                                                                              >  4:225913/1‑84 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaaa                                                                                                                              <  3:225913/84‑1 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                             >  1:246475/1‑149 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                             >  1:545547/1‑149 (MQ=255)
                                                                                      aaaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtttt                                                           >  2:523573/1‑149 (MQ=255)
                                                                                       aaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTc                                                          <  2:446257/149‑1 (MQ=255)
                                                                                           tGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATa                                                      >  3:236973/1‑149 (MQ=255)
                                                                                             ccGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGc                                                    <  4:579949/149‑1 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                 <  2:164027/149‑1 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                 <  2:696798/149‑1 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                 <  1:275565/149‑1 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGAGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                 <  4:68971/149‑1 (MQ=255)
                                                                                                      tGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATaa                                           <  1:40889/149‑1 (MQ=255)
                                                                                                       gTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAg                                          <  4:434083/149‑1 (MQ=255)
                                                                                                       gTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAg                                          <  2:214726/149‑1 (MQ=255)
                                                                                                       gTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAg                                          <  1:675637/149‑1 (MQ=255)
                                                                                                         cATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCt                                        >  4:373730/1‑149 (MQ=255)
                                                                                                          aTTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTg                                       >  2:241759/1‑149 (MQ=255)
                                                                                                            tttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTc                                     <  1:601592/149‑1 (MQ=255)
                                                                                                                 acacTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTaaa                                >  1:395415/1‑149 (MQ=255)
                                                                                                                   acTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACg                              >  3:262055/1‑149 (MQ=255)
                                                                                                                    cTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGc                             <  4:453114/149‑1 (MQ=255)
                                                                                                                     tCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGcc                            >  3:333575/1‑149 (MQ=255)
                                                                                                                      cATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCa                           <  4:333575/149‑1 (MQ=255)
                                                                                                                       aTCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTc                                                          <  4:685289/117‑1 (MQ=255)
                                                                                                                       aTCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTc                                                          >  3:685289/1‑117 (MQ=255)
                                                                                                                       aTCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAg                          >  1:39801/1‑149 (MQ=255)
                                                                                                                       aTCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAg                          >  3:298637/1‑149 (MQ=255)
                                                                                                                        tCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGc                         >  3:376751/1‑149 (MQ=255)
                                                                                                                         cGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCg                        >  1:90603/1‑149 (MQ=255)
                                                                                                                           acacGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAt                      >  4:354005/1‑149 (MQ=255)
                                                                                                                           acacGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAt                      <  2:644655/149‑1 (MQ=255)
                                                                                                                              cGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAtttt                   <  2:487620/149‑1 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCATTTGGCTTAAATTCACCGGCACATTGTTTTCAACAGCGGTGCATAAGCTGTCCTAAACCCCAGCGATTTTTTc                <  4:579857/149‑1 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                <  4:34793/149‑1 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                <  2:90603/149‑1 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                <  3:592238/149‑1 (MQ=255)
                                                                                                                                   aTCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGTTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCga              >  3:461734/1‑149 (MQ=255)
                                                                                                                                   aTCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCga              >  1:280952/1‑149 (MQ=255)
                                                                                                                                      gttttgcAGGGGTTACGCCTGCAAAACGGGCAAAGCCTCAGCGCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGAt           >  4:752640/7‑149 (MQ=255)
                                                                                                                                      cccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGAt           >  2:174900/1‑149 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa          >  2:527010/1‑149 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa          <  3:535578/149‑1 (MQ=255)
                                                                                                                                       ccAAACAGGCGGAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa          >  3:655554/1‑149 (MQ=255)
                                                                                                                                          tgcAGGCGTTACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCa                           >  1:556237/3‑129 (MQ=255)
                                                                                                                                          tgcAGGCGTTACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCa                           <  2:556237/127‑1 (MQ=255)
                                                                                                                                             aGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTc    <  2:246475/149‑1 (MQ=255)
                                                                                                                                             aGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTc    <  1:272794/149‑1 (MQ=255)
                                                                                                                                             aGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTc    >  3:580249/1‑149 (MQ=255)
                                                                                                                                               gCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCa  <  2:397629/149‑1 (MQ=255)
                                                                                                                                               gCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCa  >  3:535458/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCA  >  NC_000913/2849798‑2850089

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: