Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,849,945 A→T 7.7% 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 7.7% 494.3 / 23.0 143intergenic (‑96/+29)hycB/hycAhydrogenase 3, Fe‑S subunit/regulator of the transcriptional regulator FhlA
Reads supporting (aligned to +/- strand):  ref base A (78/53);  new base T (4/7);  total (82/61)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 149 total aligned reads displayed.

AAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAAC  >  NC_000913/2849798‑2850093
                                                                                                                                                   |                                                                                                                                                    
   tgtgGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGc                                                                                                                                                  <  2:900602/149‑1 (MQ=255)
         aGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgttt                                                                                                                                            <  4:279247/149‑4 (MQ=255)
         aGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCaa                                                                                                                                            >  2:152902/1‑149 (MQ=255)
           ccGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttgg                                                                                                                                          <  1:474407/149‑6 (MQ=255)
             gATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACg                                                                                                                                        >  4:775418/1‑149 (MQ=255)
              aTACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAATCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCAGCCCCAAACAGGCCTGACGCCTGCAAAACgg                                                                                                                                       <  4:572140/149‑1 (MQ=255)
                aCAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGc                                                                                                                                     >  4:199087/1‑149 (MQ=255)
                  agagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaa                                                                                                                                   >  4:90636/1‑149 (MQ=255)
                   gagCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCaaa                                                                                                                                  <  2:311201/149‑1 (MQ=255)
                    agCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAg                                                                                                                                 >  1:730815/1‑149 (MQ=255)
                     gCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGc                                                                                                                                <  3:121651/149‑1 (MQ=255)
                      cGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCACTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGcc                                                                                                                               >  3:261341/1‑149 (MQ=255)
                           aGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGc                                                                                                                          <  4:678694/149‑1 (MQ=255)
                             tCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTc                                                                                                                        >  3:622081/1‑149 (MQ=255)
                               aGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttg                                                                                                                                           <  4:681913/128‑5 (MQ=255)
                               aGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttg                                                                                                                                           >  3:681913/1‑124 (MQ=255)
                               aGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAt                                                                                                                      >  3:472636/1‑149 (MQ=255)
                                  aaTTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgct                                                                                                                   >  1:210817/1‑149 (MQ=255)
                                  aaTTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgct                                                                                                                   >  1:549139/1‑149 (MQ=255)
                                  aaTTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgct                                                                                                                   >  2:350999/1‑149 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttgg                                                                                                                                          >  4:55801/1‑115 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgtttgg                                                                                                                                          <  3:55801/120‑6 (MQ=255)
                                        aaaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCggg                                                                                                             >  1:193511/1‑149 (MQ=255)
                                         aaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGc                                                                                                            >  3:700842/1‑149 (MQ=255)
                                          aaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                           >  4:325978/1‑149 (MQ=255)
                                          aaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                           >  1:412568/1‑149 (MQ=255)
                                                 cACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGTGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTccc                                                                                                    >  1:695236/1‑149 (MQ=255)
                                                 cACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTccc                                                                                                    <  3:90636/149‑1 (MQ=255)
                                                  aCCAGGCATTCTTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCACGCTGCCGGGCTTTGTCCCt                                                                                                   >  3:524800/1‑149 (MQ=255)
                                                  aCCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACTCGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCt                                                                                                   >  4:604896/1‑149 (MQ=255)
                                                  aCCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCt                                                                                                   <  3:115700/149‑1 (MQ=255)
                                                   ccAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCtt                                                                                                  >  3:870135/1‑149 (MQ=255)
                                                   ccAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCtt                                                                                                  <  3:29657/149‑1 (MQ=255)
                                                   ccAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCtt                                                                                                  >  4:873118/1‑149 (MQ=255)
                                                   ccAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTATGTCCCtt                                                                                                  <  3:852567/149‑1 (MQ=255)
                                                   ccAGGCATTCCTCAGGTGATAGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCtt                                                                                                  >  2:270022/1‑149 (MQ=255)
                                                    cAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                           <  3:907912/139‑1 (MQ=255)
                                                    cAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCt                                                                                                           >  4:907912/1‑139 (MQ=255)
                                                    cAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCttt                                                                                                 <  1:639434/149‑1 (MQ=255)
                                                     aGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTa                                                                                                >  4:847089/1‑149 (MQ=255)
                                                            ccTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTg                                                                                         <  2:773051/149‑1 (MQ=255)
                                                                   tGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGTAAAACGGGCAAAGCCTCAGCTCACGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTaaa                                                                                  <  4:524800/149‑1 (MQ=255)
                                                                   tGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTaaa                                                                                  <  2:484565/149‑1 (MQ=255)
                                                                    gATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTaaa                                                                                  <  2:461124/148‑1 (MQ=255)
                                                                    gATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTaaa                                                                                  >  1:461124/1‑148 (MQ=255)
                                                                     aTTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAAtt                                                                                >  3:531256/1‑149 (MQ=255)
                                                                         tCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGc                                                                                                                                     <  3:179268/92‑1 (MQ=255)
                                                                         tCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGc                                                                                                                                     >  4:179268/1‑92 (MQ=255)
                                                                           aTTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgttt                                                                                                                                            <  4:501103/83‑4 (MQ=255)
                                                                           aTTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTTACGCCTgttt                                                                                                                                            >  3:501103/1‑80 (MQ=255)
                                                                            tttttGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGc                                                                         >  2:622686/1‑149 (MQ=255)
                                                                                 gACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACAtt                                                                    <  4:115671/149‑1 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTaa                                                                                                                                                     <  3:846751/65‑1 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTaa                                                                                                                                                     >  4:846751/1‑65 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCt                                                                                                                                                >  4:458797/1‑70 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCt                                                                                                                                                <  3:458797/70‑1 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                                 >  4:862283/1‑149 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                                 >  2:297949/1‑149 (MQ=255)
                                                                                    gAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtt                                                                 <  2:351045/149‑1 (MQ=255)
                                                                                      aaaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGtttt                                                               <  2:730815/149‑1 (MQ=255)
                                                                                          aTGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAAt                                                           >  2:17846/1‑149 (MQ=255)
                                                                                              cGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCg                                                       >  4:12139/1‑149 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                     <  4:531256/149‑1 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                     >  4:431457/1‑149 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                     <  4:870135/149‑1 (MQ=255)
                                                                                                aTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                     <  3:431457/149‑1 (MQ=255)
                                                                                                      tGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATaa                                               >  1:519023/1‑149 (MQ=255)
                                                                                                       gTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAg                                              >  3:839013/1‑149 (MQ=255)
                                                                                                           ttttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtg                                                                                                                     >  1:386064/1‑74 (MQ=255)
                                                                                                                gACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTaa                                     >  1:167243/1‑149 (MQ=255)
                                                                                                                   acTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACg                                  >  3:435284/1‑149 (MQ=255)
                                                                                                                   acTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACg                                  >  4:331686/1‑149 (MQ=255)
                                                                                                                   acTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACg                                  >  4:455398/1‑149 (MQ=255)
                                                                                                                     tCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGcc                                >  3:553208/1‑149 (MQ=255)
                                                                                                                       aTCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGc                                 <  1:320109/146‑1 (MQ=255)
                                                                                                                       aTCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAg                              >  3:800686/1‑149 (MQ=255)
                                                                                                                               gCCCATCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAtttttt                     >  2:23367/1‑149 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                    <  4:435284/149‑1 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                    <  2:695706/149‑1 (MQ=255)
                                                                                                                                 ccATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTc                    >  3:29131/1‑149 (MQ=255)
                                                                                                                                   aTCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCga                  >  3:774420/1‑149 (MQ=255)
                                                                                                                                      gttttgcAGGCGTTACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCa                               >  4:868483/7‑133 (MQ=255)
                                                                                                                                      gttttgcAGGCGTTACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCa                               <  3:868483/127‑1 (MQ=255)
                                                                                                                                      cccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCt                                                                                      <  1:424353/78‑1 (MQ=255)
                                                                                                                                      cccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCt                                                                                      >  2:424353/1‑78 (MQ=255)
                                                                                                                                      cccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGAt               >  2:553536/1‑149 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa              <  4:273856/149‑1 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa              >  2:697139/1‑149 (MQ=255)
                                                                                                                                       ccAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATa              >  2:781978/1‑149 (MQ=255)
                                                                                                                                           aCAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAgtgt          <  2:651647/149‑1 (MQ=255)
                                                                                                                                             aGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTc        <  3:604896/149‑1 (MQ=255)
                                                                                                                                             aGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTAAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTc        <  4:472636/149‑1 (MQ=255)
                                                                                                                                                 gTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGa    >  2:892563/1‑149 (MQ=255)
                                                                                                                                                  tAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGaa   >  2:850996/1‑149 (MQ=255)
                                                                                                                                                  tAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGaa   <  4:615399/149‑1 (MQ=255)
                                                                                                                                                  tAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGaa   <  2:39837/149‑1 (MQ=255)
                                                                                                                                                  tAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGaa   <  2:664524/149‑1 (MQ=255)
                                                                                                                                                  tAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATGCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGaa   <  4:768474/149‑1 (MQ=255)
                                                                                                                                                   aaCGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAAc  >  3:870681/1‑149 (MQ=255)
                                                                                                                                                   aaCGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAAc  <  1:152902/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAAC  >  NC_000913/2849798‑2850093

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